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ABSTRACT 


The  first  part  of  this  work  (Chapters  2 and  3)  considers  the  problem 
of  radiation  from  sources  in  the  presence  of  smooth,  convex,  impenetrable 
objects  and  presents  a brief  survey  of  various  high-frequency  techniques. 

A generalization  of  the  geometrical  theory  of  diffraction  and  two  new 
techniques  — based  on  the  spectral  domain  approach  and  an  asymptotic 
evaluation  of  the  radiation  integral  for  the  surface  current  — also  are 
discussed.  Some  numerical  results  derived  from  the  spectral  domain  formulas 
are  presented,  and  a comparison  with  available  theoretical  and  experimental 
data  is  included. 

The  second  part  of  this  work  (Chapters  4 through  10)  describes  a 
new  technique  for  analyzing  the  radiation  from  a point  electromagnetic 
source  located  on  an  infinitely  conducting  solid  cylinder  with  a finite 
length.  The  method  presented  in  this  work  is  based  upon  a generalization 
of  STD  to  the  case  of  curved  surfaces. 
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1.  INTRODUCTION 


The  problems  of  radiation  from  sources  in  the  presence  of  impenetrable 
smooth  convex  objects  and  the  diffraction  of  a plane  wave  by  such  objects  are 
of  great  practical  interest  in  the  design  of  antennas  on  structures,  e.g., 
conformal  arrays.  Unfortunately,  the  exact  analytical  solutions  to  these 
problems,  based  on  the  methods  of  "separation  of  variables"  or  "function- 
theoretic"  procedures  (Wiener-Hopf  technique,  residue  calculus,  etc.),  exist 
only  for  a very  limited  number  of  scattering  geometries.  Furthermore, 
the  exact  solutions  are  typically  highly  complex  in  nature;  hence, the  process 
of  extracting  numerical  results  from  them  can  be  very  time-consuming  and  is 
by  no  means  trivial.  This  situation  has  motivated  many  researchers  to  explore 
approaches  to  the  problems  of  radiation  and  scattering  from  smooth  convex 
structures . 

In  the  low  and  resonant  frequency  ranges,  several  reliable  numerical 
procedures,  e.g.,  the  moment  method,  are  available  for  solving  the  radiation 
and  scattering  problems.  However,  in  the  high  frequency  domain,  numerical 
techniques  based  on  matrix  methods  become  unwieldy,  if  not  impractical,  prompting 
one  to  employ  asymptotic  techniques  suitable  for  large  k(«2-:/\),  where  \ is 
the  wavelength  of  the  illuminating  wave. 

In  this  work,  we  begin  by  presenting,  in  Chapter  2,  a survey  of  various 
high  frequency  asymptotic  techniques  for  the  problem  stated  above.  The  survey 
will  be  necessarily  brief,  and  will  cover  only  the  highlights  of  a number  of 
important  approaches  to  the  problem  at  hand,  viz.,  Fock's  theory,  the  geometrical 
theory  of  diffraction  (GTD) , and  the  direct  integral  equation  approach.  The 
reader  interested  in  further  details  may  choose  to  consult  the  works  of  3owman, 
et  al.  [1],  Uslenghi  [2],  and  Kouyoumjian  [3]. 
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In  Chapter,  3,  we  consider  the  generalization  of  GTD  and  present  some 
new  approaches  to  the  curved  surface  radiation  and  scattering  problems.  Some 
numerical  results  based  on  one  of  these  new  approaches  are  presented  in  Chapter 
3,  and  a comparison  with  other  available  methods  is  included. 

Throughout  the  rest  of  the  work,  we  apply  the  results  of  the  new 
approach  combined  with  STD  interpretation  of  wedge  diffraction  to  some  special 
case  of  great  theoretical  and  practical  interest.  The  high  frequency  radiation 
of  an  electromagnetic  point  source  on  the  surface  of  conducting  circular 
cylinder  with  finite  length  is  analyzed  by  a new  technique  combining  the  main 
features  of  different  asymptotic  theories.  Chapter  4 is  concerned  with  the 
general  formulation  of  the  problem.  Chapters  5,  6,  and  7 describe  the  methods 
used  in  calculating  different  "constituents"  of  the  total  radiated  field, 
namely,  surface  and  wedge  diffracted  rays.  Chapter  8 is  concerned  with 
analyzing  the  effect  of  the  finiteness  of  the  cylinder.  Field  evaluation 
in  paraxial  region  is  explained  in  Chapter  9.  Numerical  results  are  discussed 
in  Chapter  10,  and  finally,  a brief  conclusion  is  given  in  Chapter  11. 


3 


2.  SURVEY  OF  AVAILABLE  HIGH-FREQUENCY  ASYMPTOTIC  TECHNIQUES 
2. 1 Watson  Trans format  Ion  and  Physical  Optics 

Ona  of  tha  first  successful  attempts  to  dariva  an  aaymptotlc  expansion 
for  the  far-fleld  generated  by  a point  source  located  in  tha  proximity  of  a 
conducting  surface  was  made  by  G.  N.  Watson  in  1918  [4].  His  method, 
essentially,  consisted  of  two  steps:  1)  transforming  the  original  infinite 
series  solution  into  a contour  Integral  (by  Cauchy's  residue  theorem); 

2)  deforming  the  contour  of  integration  so  as  to  capture  a sat  of  complex 
poles  of  the  Integrand.  The  original  Integral  is  then  expressed  in  terms 
of  an  infinite  series  which  converges  very  rapidly,  provided  the  observation 
point  is  in  the  shadow  region.  The  first  few  terms  of  this  series  were  later 
interpreted  as  "creeping  waves."  The  method  was  first  applied  to  a sphere 
and  circular  cylinder,  and  later  to  some  other  geometries  as  well.  The 
mathematical  rigor  of  the  method  was  the  subject  of  further  investigations  by 
other  researchers  ([5],  [6],  and  [7]).  Although  the  Watson  transformation 
can  onlv  be  applied  to  a few  simple  geometries,  e.g.,  the  sphere,  cylinder, 
cone,  spheroid,  etc.,  it  Is  still  regarded  as  one  of  the  cornerstones  of  the 
general  high  frequency  techniques  because  of  its  mathematical  rigor.  The 
Watson  transformation  Is  especially  powerful  in  the  shadow  region  of  the 
geometric  optics  field.  In  the  lit  region,  the  above-mentioned  contour 
integral  is  evaluated  using  the  "stationary  phase"  method  and  yields  the 
reflected  field  from  the  surface.  In  this  region,  the  moat  significant  contri- 
bution to  the  total  scattered  field  tvplcallv  comes  from  the  surface  current 
Induced  on  the  smooth  convex  part  of  the  object;  the  so-called  "Physical  Optics" 
approximation  can  he  applied  ((8|,  [9|,  and  [10])  to  derive  the  reflected 
field.  The  Physical  Optics  method  is  based  upon  approximating  the  induced 
surface  current  in  the  lit  region  of  the  object  by  the  current  that  would  be 


I 


4 


Induced  on  Che  local  tangent  plane,  and  by  assuming  that  the  surface  current 
is  *ero  In  the  shadow  region.  The  far  field  is  constructed  by  substituting 
the  above  estimate  for  the  induced  surface  current  in  the  integral  representa- 
tion of  the  scattered  field,  and  evaluating  the  same  in  an  asymptotic  sense. 

The  dominant  term  of  the  asymptotic  expansion  of  this  integral  can  be  shown 
to  be  identical  to  the  first  term  of  the  Luneberg-Kline  expansion  of  the 
geometrical  optics  far  field  ([111  and  [12)).  However,  the  higher-order  terms 
derived  from  the  physical  optics  aoproach  do  not  provide  us  with  correct 
results  in  the  shadow  or  transition  regions  where  the  diffracted  field  contri- 
butes the  most. 

In  the  next  subsection,  we  discuss  Fock's  theory,  which  can  fill  the 
gap  between  the  Physical  Optics  in  the  lit  region  and  the  "creeping  wave" 
representation  in  the  shadow  region. 

2.2  Fock's  Theory 

The  region  between  the  lit  and  the  shadow  part  on  a surface  is  called 

"penumbra  region."  The  angular  width  of  this  region  is  approximately  given 
1/3 

bv  (Xr^/i)  where  X is  the  wavelength  of  the  illumination  and  rQ  is  the 
radius  of  curvature  of  the  surface  of  the  object  in  this  region  in  the  incident 
plane  (Fig.  1).  Fock's  theory  invokes  the  principle  of  y joal  character  of 
the  field  in  the  penumbra  region  [13]  and  is  based  on  the  conjecture  that 
all  bodies  with  a smoothly  varying  curvature  have  the  same  current  distribution 
in  the  penumbra  region,  provided  that  the  curvature  and  the  incident  wave  are 
the  same  near  the  point  under  consideration.  This  principle  allows  one  to 
locally  replace  the  surface  of  the  object  by  a portion  of  a paraboloid  of 
revolution.  A unique  feature  of  the  expressions  for  Foci:  currents  is  that 
they  provide  a convenient  transformation  of  the  geometric  optics  currents 


Penumbra  region 


Figur*  i:  Section  of  che  body  in  eh*  plan*  of  iacidanca. 


;ura  I : 


Plana  v*v«  incidanc  upon  i smooch  :onv«x  bodv . 
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in  the  lit  region  into  the  creeping  wave  currents  in  the  shadow  region. 

Fock  himself  deduced  the  pertinent  formulas  for  the  surface  currents  by 
treating  a convex  body  problem  [14]  described  below: 

Consider  a convex  body  and  a plane  wave  incident  in  the  direction  of  the 
x-axis.  If  the  equation  of  the  surface  is 

f(x,y,z)  - 0 (2.1) 

then  the  curve  representing  the  boundary  of  geometrical  shadow  is  given  by 

f(r)  - 0 , - 0 (2.2) 

Consider  a point  0 on  the  boundary  of  a shadow  region  where  we  set  up 
a rectangular  coordinate  system  as  shown  in  Fig.  2 (z:  normal  to  the  surface, 
x:  in  the  direction  of  propagation,  and  y is  the  tangent  to  the  boundary  of 
shadow).  In  the  vicinity  of  this  point,  the  surface  of  the  body  could  be 
locally  replaced  by  a paraboloid  of  revolution  which  is  expressed  by  the 
equation. 

z + 1/2  (ax2  + 2bxy  + cy“)  -0  (2.3) 

Each  of  the  field  components  satisfies  the  Helmholtz  equation 

(72  + k2)f  - 0 (2.4) 

The  fact  that  the  incident  wave  travels  along  the  x-axis,  suggests  that  f 
be  written  in  the  form 

I 


y . y -Jk* 


(2.5) 


where  an  exp(Jwt)  time  dependence  has  been  assumed.  Substituting  (5)  In  (4) 
gives 

3"  ^ ' 

(2.6) 


v‘y  - 21k  •—  '»  - 0 

3x 


At  this  point,  two  basic  assumptions  are  introduced  in  Fock’s  theorv. 


viz. 


i)  f’s  are  relatively  slowly  varying  function  of  coordinates 
il)  ')  varies  more  rapidly  in  the  z-direction  than  in  x and  y, 
i. e. , 


f - o(S  ?>.  f • o,\  ?>  H-o  <*,.> 

3z  m ^x  m 3v  m’ 


(2.7) 


Based  upon  (2.7),  we  can  write  (2.6)  as 


■*  y jy 

- — i - 21k  — - 0 
3z*  J 3x 


(2.8) 


and  consequently  m'  - m~  (m  is  very  large),  where  the  terms  of  relative  order  1/m" 
have  been  omitted. 

Inserting  these  estimates  and  assumptions  into  the  Maxwell’s  equations, 
we  can  find  some  simple  expressions  for  all  the  field  components  in  terms  of 


H and  H . If  we  write  H as 

v z v 


H - H°  «-JI“  . 

V V 


(:.<») 
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where  H°  is  Che  magnitude  of  the  Incident  wave  at  infinity,  then  Y must 
satisfy 


gV 

where  H ‘ ■ external  field 

y 

G(C,q)  » e'(J/3)  5 V^r.q) 

V^(i,q)  * Fock  function  defined  in  the  Appendix  A. 

i - m(ax  + by)  - reduced  distance  from  the  shadow  boundary  - i/d. 

m - (kr  /2)1,  3 , m-1-  is  the  angular  width  of  penumbra  region, 
o 

2 1/3 

d * che  width  of  penumbra  region  ■ (2r^  /k) 

i - distance  between  the  observation  point  and  the  shadow  boundary  along  the 


incident  ray  (Fig.  3). 
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q ■ -1m/» n--(J/* n)  \rr“  (*  0 for  conducting  body), 


The  other  tangential  component  of  the  magnetic  field  11^  on  the  surface 
of  the  bodv  can  be  obtained  in  a similar  manner 


H - Hax  [-  ^ e'(3/3)C  /3  f(S)  ] 


x z 


(2.13) 


where  f(£)  is  another  Fock  function  defined  in  Appendix  A.  Fock's 
formulas  not  only  give  the  surface  value  of  the  field,  but  also  can  be 
utilized  to  find  the  field  in  the  proximity  of  the  object.  For  a plane  wave 
incidence,  the  first  order,  i.e.,  0(l/m)  terms  for  the  scattered  field  within 
a certain  layer  around  the  object,  can  be  written  as 


Hx  - 0,  Hy  - e"JkX  * U,;>.  H - V?  e"3kX 


E ■ (j  /m)  H°  e"Jkx  E - H , E - -H 

X V V 2 2 V 


(2. 14) 


where 


IT  2 

; ■ 2am“[z  +d/2)<ax*‘  + 2bxy  > cy  ) ] • reduced  height  from  the  surface 
of  the  body  (see  Fig.  •*) . 


,3, 


O j ,-J5t  it 

(2.15) 


• -Kt  W’U) 

4 - f e J"c 


dt 


The  path  of  integration  for  I and  f is  shown  in  Fig.  5 . 


Complex  t- plane 


Fijure  5:  Path  o<  inctfracion  for  12.151  in  the  complex  t-plane. 


Jcapari son  becveen  the  •srtoua  definitions  o: 
taraaecar  ; for  the  case  c:  a tirtuiar  fliacer. 
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Fock's  expressions  for  the  field  components  in  the  penumbra  region 
(5  - 0)  can  be  extended  to  the  shadow  region,  by  introducing  some  modifications 
in  the  definition  of  parameter  ?.  Goodrich  [15]  has  generalized  the 
argument  used  by  Fock  in  the  penumbra  region  to  anywhere  in  the  shadow 
region  by  introducing  a new  set  of  variables,?  and  c,for  the  incremental 
distances  along  the  path  leading  into  the  shadow  region.  In  this  generalization, 
the  parameter  ? as  defined  in  (2.12)  is  replaced  by 


, fS  ,kR(sK  1/ 3 ds 
5 i 1 2 ' R(s) 


(2.16) 


where  s is  the  arc  length  along  the  geodesics  which  originate  from  the  shadow 
boundary  and  go  into  the  shadow  region  along  the  surface,  and  R(.s)  is  the 
radius  of  curvature  of  the  surface  along  the  geodesics.  For  the  case  of  a 
circular  cylinder  of  radius  a (Fig.  6),  the  expression  of  ? simplifies  to 


£ - (ka/2)1/3  9 - s/d 


(2.17) 


Fock  also  treated  the  case  where  the  point  source  was  very  close  to 
the  surface  of  the  body.  He  analyzed  the  radiation  of  electric  dipoles  near  a 
spherical  model  of  the  earth  [16]  and  derived  the  formulas  for  the 
scattered  fields  in  terms  of  functions  (attenuation  functions)  similar  to  41 
and  ¥,  which  are  valid  both  in  the  shadow  and  transition  regions  [17]. 

Fock's  assumptions  were  later  proven  in  a more  systematic  and  mathematically 
rigorous  manner  by  Cullen  [18]  and  Hong  [19]  by  using  a direct  integral 
equation  approach.  This  method  is  described  in  the  next  section. 
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2. 3 Direct  Integral  Equation  Approaches 

This  method,  which  is  closely  related  to  Fock's  theory,  can  be 
illustrated  by  analyzing  the  diffraction  of  a plane  electromagnetic  wave 
by  an  arbitrary  conducting  body  (large  compared  with  X).  Cullen  [18]  obtained 
a first-order  asymptotic  solution  to  the  integral  equation  for  the  induced 
surface  current 

J(r)  - 2n(r)  * HinC(?)  - (l/2*)nU) 


(2.18) 


/ / ds'  (J(r')  x Re“jkR} 


where  n(r)ls  the  outward  unit  normal  to  the  surface  at  r,  H (r)  is 
the  incident  magnetic  field  on  the  surface  (S)  of  the  body, and  R - r - r' 

(r'  is  a variable  point  on  the  surface). 

Fock  used  this  integral  equation  to  deduce  the  important  principle 
of  local  character  of  the  field  in  the  penumbra  region.  Cullen  derived  a 
first-order  asymptotic  solution  to  (2.18)  which  agreed  with  Fock's  results 
given  in  (2.12)  and  (2.13).  Cullen's  method  consists  of  transforming  the  two- 
dimensional  integral  equation  (2.18),  in  the  penumbra  region,  to  a one-dimensional, 
Volterra-type  equation.  This  is  accomplished  by  applying  the  stationary 
phase  technique  to  the  original  Integral  while  Integrating  with  respect  to 
one  of  the  variables.  The  resulting  one-dimensional  Volterra  equation  is  then 
solved  in  Cullen's  method  by  the  Fourier  transform  technique.  A similar 
procedure  was  used  by  Hong  [19]  to  analyze,  asymptotically , the  diffraction 
of  electromagnetic  and  acoustic  plane  waves  bv  smooth  convex  bodies.  We  will 
now  proceed  to  explain  Hong's  method  in  a little  more  detail  by  referring 
back,  once  again,  to  the  integral  equation  (2.181.  The  surface  is  parametrized 
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by  the  geodesic  coordinate  system  (o,v)  such  that  the  shadow  boundary  for  the 
Incident  plane  wave  traveling  along  the  tangent  o(0,v)  to  the  v ■ 0 curve  is  the 
c **  0 curve.  The  quantities  o(o,v),  b(o,v)  and  n(a,v)  form  a right-hand  local 
orthonormal  basis  (n  ■ o * b)  as  shown  in  Fig.  7. 

Since  the  incident  field  has  a phase  factor  e Jko(O.O)  r(.,0)^ 
we  write  the  surface  current  in  the  form 


J(r)  - [ I ^ <r)  o(r>  + I.  (r)  b(r)]  e~Jko 

O D 


(2.19) 


where  a is  Che  arc  length  along  the  geodesic.  Substituting  (2.19)  back  into 

(2.18)  and  restricting  the  resulting  equation  to  the  points  on  the  geodesic  v»0, 

we  obtain  two  coupled,  two-dimensional  integral  equations  for  I (o,0)  and 

0 

l^(j,0).  It  can  be  shown  that  these  integrals  have  saddle  points  at  v-0  (for 

the  v-integration) . Applying  the  "steepest  descent  path"  method  to  v-integration, 

2 1 1 3 

and  keeping  the  terms  up  to  the  order  1/M“,  where  ■ ike  _ (o,0))  ' , we 

obtain  the  following  decoupled  one-dimensional, Volterra-tvpe  integral  equations 

for  I (5,0)  and  l (5,0) 

D J 


I U.O)  - 2 l^nc(i,0)  - / dr  I (5,0)  K (5-r)  + 0 (>T3) 

J J 0 0 u 

— oo 

(2.20) 

Ib(5,0)  - 2 Ibnc  (5,0)  - f dr  Ib(5,0)  ^(5-t)  + 0(b-“3) 


e (<7,v)  is  the  radius  of  curvature  of  the  surface  along  geodesics  (v  » constant 
curves)  at  point  (a.vl. 

Solving  (2.20)  by  Fourier  transforms,  we  obtain  the  expression  for  the 
induced  currents  in  the  penumbra  and  shadow  regions, and  tire  first-order  solutions 


*-  ^ J 
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are  found  to  be  the  same  as  those  of  Fock  and  Keller's  GTD  [20].  One 
of  the  important  conclusions  drawn  from  Hong's  solution  is  that  the  leading 
term  in  the  asymptotic  expansion,  which  is  the  same  for  the  acoustic  and  electro- 
magnetic problems,  is  independent  of  curvature  in  the  direction  transverse  to  the 
geodesic,  provided  the  divergence  factor  is  suppressed.  However,  we  should 
bear  in  mind  that  Hong's  method  was  designed  for  the  case  of  axial  incidence  on 
symmetric  objects,  and  in  this  case,  the  geodesics  are  torsionless.  The  above 
conclusion  does  not  seem  to  be  valid  in  the  cases  where  the  rays  have  nonzero 
torsion  ([21]  and  [22]).  In  Hong's  expressions  for  the  surface  current,  the 

transverse  curvature  has  only  a second-order  effect.  It  was  also  shown  that  up 

-2/3 

to  the  terms  of  order  (kp^)  in  the  asymptotic  expansion,  the  tangential  and 
binormal  components  of  the  creeping  waves  are  not  coupled. 

Both  Fock's  theory  and  the  "direct  integral  equation  approach"  give  the 
induced  surface  current,  or  the  scattered  field  in  the  neighborhood  of  the 
surface  of  the  scatterer,  due  to  an  incident  plane  wave.  These  expressions  can 
also  be  used  to  derive  the  radiated  field  via  the  use  of  the  reciprocity  theorem. 

The  methods  which  have  been  discussed  thus  far  are  mathematically 
rigorous.  However,  they  are  limited  in  the  scope  of  their  application  to 
geometries  satisfying  some  special  smoothness  and  symmetry  criteria.  "Geometrical 
theory  of  diffraction"  (GTD),  which  we  discuss  in  the  next  section,  has  a 
broader  scope,  although  it  does  lack  the  mathematical  rigor  of  approaches 
described  until  now. 

2.4  Geometrical  Theory  of  Diffraction  (GTD) 

Geometrical  theory  of  diffraction  (GTD),  developed  by  J.  B.  Keller 
([20],  [23],  [24],  [25],  and  [26])  is  a generalization  of  geometrical  optics. 
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It  is  based  upon  the  assumption  that  fields  propagate  along  rays.  Keller's 
major  contribution  was  to  introduce  the  new  kinds  of  rays  called  the  "diffracted 
rays,"  which  together  with  the  geometrical  optics  rays, constitute  the  total 
field.  In  our  problem,  viz.,  source  radiating  in  the  proximity  of  the  smooth 
object,  the  diffracted  rays  travel  along  the  curves  on  the  surface  of  the 
scatterer.  By  applying  Fermat's  principle  to  these  surface  rays,  we  conclude 
that  the  above-mentioned  curves  should  be  geodesics  on  the  surface  of  the  body. 

In  the  GTD  procedure,  one  assigns  a value  to  the  field  along  each  of  these  sur- 
face rays. The  total  field  at  any  point  in  the  space  is  the  sum  of  the  fields 
due  to  various  rays  (incident,  reflected  and  diffracted)  passing  through  that 
point.  An  important  advantage  of  the  GTD  approach  is  that  it  can  be  applied 
to  both  scalar  (acoustic)  and  vector  (electromagnetic)  problems  and  to  smooth 
convex  objects  of  an  arbitrary  shape. 

Consider  the  problem  of  determining  the  radiated  field  of  a scalar  point 
source  located  on  the  surface  of  a smooth  convex  opaque  body.  If  the  observation 
point  is  in  the  shadow  region,  the  ray  paths  originating  at  Q and  reaching  P 
(ooservation  point)  are  comprised  of  two  sections.  One  of  the  sections  follows 
straight  line  path  P^P,  while  the  other  travels  along  a geodesic  on  the  surface 
(Fig.  8).  Let  us  consider  the  propagation  of  the  field  along  each  section,. 

a)  Rays  in  free  space:  Behavior  of  the  fields  along  these  rays  can  be 
determined  by  obtaining  a high-frequency  asymptotic  solution  to  Maxwell's 
equation  in  a source-free  homogeneous  isotropic  medium.  We  begin  with  the 
Luneberg-Kline  asymptotic  expansion  of  the  electric  field  ([11]  and  [12]): 

-Jk  S(r)  » 

E(r)  - k e ° T (jk)  e (r)  (2.21) 

in 

m-0 
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and  insert  it  into  the  Maxwell's  equations.  This  results  in  the  following 
equations  governing  the  propagation  of  electromagnetic  fields  along  the  rays. 


[VS(r)]“  - 1 (Eikonal  equation) 


(2.22) 


2(VS*V)e  + (V2S)e  - -V2  e (Transport  equation) 
m m m-i 


(2.23) 


VS’e  ■ -7*e  , (Gauss's  Law) 

m in—  i 


(2.24) 


e_1  - 0,  m - 0,  1,  2,  . . 


The  zeroth-order  solution  to  the  above  system  of  equations, which  turns  out 
to  be  in  agreement  with  what  one  would  obtain  by  geometric  optics,  may  be 
written  as 


- * -Jk0s(0) 

E(o)  ~ fiQ(0)e  U 


(pj+o) (p „+c) 


e-jV 


(2.25) 


where  o is  the  distance  traveled  along  the  ray  from  the  reference  point 

0(o«0)  on  the  ray  path  (Fig.  9).  p^  and  are  the  principal  radii  of  curvature  of 

the  wavefront  at  o» 0.  It  is  apparent  that  the  expression  fails  when  o»-o^  or 

a " -02’  at  c^e  caustic  lines  (Fig.  9).  In  the  cases  where  it  is 

convenient  to  choose  the  point  of  diffraction  on  the  surface  of  a bodv  as 

the  reference  point  0.  the  formula  (2.25)  should  be  modified  as  follows: 


a (c+o) 


(2. 26) 


t 


d 
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In  these  cases,  the  point  of  diffraction  itself  Is  a caustic,  and  o Is  the 
distance  between  this  point  and  the  second  caustic. 

b)  Surface  Rays:  These  rays  follow  the  surface  S along  the  geodesics 
into  the  shadow  region,  and  shed  off  energy  tangentially  as  they  propagate. 

In  order  to  study  the  behavior  of  the  field  along  these  rays,  we  introduce  a 
special  ray-fixed  coordinate  system, o,n,b. 

» A 

o:  Unit  vector  tangent  to  the  ray;  n:  outward  unit  normal  to  the 
surface;  and  b ■ t * n or  binormal  direction;  a vector  field  can  be  decomposed 
into  its  components  along  these  unit  vectors  as 
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jU  -ko) 

a(a)  ■ A(a)  e 


(2.28) 


where  a is  the  distance  between  an  arbitrary  point  along  the  ray  and  the  source 
Q and  <J>Q  is  the  phase  of  the  field  at  the  source  point.  Next,  invoking  the 
principle  of  conservation  of  energy  between  two  adjacent  rays,  and  using  the  fact 
that  the  surface  rays  shed  energy  off  tangentially,  we  can  arrive  at  the 
following  expression  for  a(o) 


a(a) 


K 


exp [ -j  kc 


(o')do'] 


(2.29) 


where  a (a)  is  the  "attenuation  constant,"  K is  proportional  to  the  strength  of 

1/ 1 

the  source,  and  df^,  dV,,  and  d are  shown  in  Fig.  10.  The  quantity  [df  (pdH^)  ] 
indicates  the  "spreading  of  the  surface  ray  tube"  as  it  travels  along  the  surface. 
Equations  (2.26)  and  (2.29)  describe  the  laws  of  propagation  for  the  rays  which 
originate  from  the  source  point  Q,  are  diffracted  at  P^,  and  reach  the 
observation  point  P.  To  complete  the  solution,  we  need  to  determine  the  actual 
values  of  the  fields  from  these  equations.  These  require  the  knowledge  of  Cq 
and  K,  which,  in  turn,  are  related  to  the  initial  values  of  the  rays  QP^  and 
P^P  as  well  as  the  attenuation  constant  x ( <r ) . The  initial  value  of  the  field 
at  0 is  related  to  the  strength  of  the  source  by  L(Q),  the  so-called  "launching 
coefficient,"  while  the  initial  value  of  the  field  at  P^  is  related  to  the 
actual  field  on  the  surface  at  P^  through  the  "diffraction  coefficient"  D(P^). 

If  we  now  sum  up  the  contributions  of  all  the  modes,  we  obtain  the  final  solution 


[25]  for  the  field  radiated  in  the  shadow  region  by  an  infinitesimal  magnetic 
dipole  of  strength  M located  on  a smooth  convex  conducting  body 
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eV)  - H • [MQJnfP^F  + jfQJbfP^G]  J~ 


s(c+s) 


-Jks 


(2. 30a) 


where 


F 


I Lph(Q)D^  (Pl)  exp  [-/  xJ(o’)do’] 

p-1  ° 


(2.30b) 


and  the  expression  for  G is  obtained  by  replacing  the  superscript  "h"  by 

"s"  in  (2.30b),  where  h and  s stand  for  hard  and  soft  boundary  conditions, 

viz.,  u-0  and  3u/3n  ■ 0,  respectively.  The  quantities  L^,s,  D^,S  and 
h s 

a ’ , appearing  in  (2.30b),  in  general,  depend  upon  the  local  geometry  and  the 
electromagnetic  characteristic  of  the  surface,  frequency  k,  and  the  mode  of 
propagation.  They  are  determined  by  studying  the  asymptotic  expansions  of 
the  exact  solutions  of  some  special  canonical  problems.  Keller  and  Levy 
([20]  and  [27]),  have  derived  the  first  few  terms  of  the  asymptotic  expansions 
for  D and  a by  considering  the  canonical  problems  of  scalar  diffraction  by  a 
circular  cylinder,  sphere,  elliptical  and  parabolic  cylinder.  A study  of  the 
above-mentioned  asymptotic  expansions  and  the  works  of  Franz  and  Klante  [28] 
and  Voltmer  [29],  who  have  also  investigated  the  same  problem,  as  well  as  a 
comparison  with  the  results  of  the  "direct  integral  equation  approach," 
reveals  the  following  characteristics  of  the  solution:  i)  the  first-order  terms 
in  the  asymptotic  expansion  of  D and  a are  independent  of  whether  the  problem 
under  consideration  is  scalar  or  vector;  ii)  the  first-order  approximation  of 
D and  a are  dependent  only  on  o , the  radius  of  curvature  of  the  surface  along 

the  ray;  iii)  the  second-order  terms  are  functions  not  onlv  of  c , but  also  of 

3 

, ‘‘c,  and  c (the  radius  of  the  curvature  of  the  surface  transverse 

~ir~  ~ ** 

to  the  ray).  Finally,  the  higher-order  terms  are  different  for  scalar  and 
vector  problems. 
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The  leading  terms  in  the  asvmpcocie  expansion  of  "diffraction  coefficient" 
D,  attentuation  constant"  a and  "launching  coefficient"  L are  presented  below: 
"Soft"  polarization: 


(2.31) 


(2.32) 


(2.33) 


(2.34) 


(2.35) 


(2. 3681 


(2. 36b) 
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and  Ai(-r  ) ■ 0,  Ai’(-r’)  ■ 0,  (Ai’  is  the  derivative  of  Ai  with  respect 
P P 

to  its  argument).  Higher-order  terms  in  the  expansion  of  D,  a and  L have 
been  given  in  [25]  and  [30]  and  in  some  of  the  other  works  on  GTD  mentioned 
earlier. 

The  expression  (2.30)  is  convenient  to  use  in  the  shadow  region.  How- 
ever, in  the  shadow  part  of  the  transition  region,  since  the  exponential  decay 
of  the  terms  in  (2.30)  is  weak,  the  convergence  of  the  series  representation  is 
very  slow.  Furthermore,  the  series  diverges  in  the  lit  part  of  the  transition 

region.  Consequently,  in  these  regions,  it  is  more  reasonable  to  use  an  integral 
representation  for  the  surface  ray  field,  which,  in  our  case,  can  be  expressed 

in  terms  of  Fock  functions  [25]. 

Attempts  have  been  made  to  establish  the  mathematical  validity  of  GTD 
and  to  minimize  its  "nondeductive  parts"  (parts  which  are  based  upon  physical 
intuition  or  the  study  of  the  asymptotic  solution  of  some  simple  problem 
geometrical  concepts  of  different  kinds  of  rays,  diffraction  coefficients, 
attenuation  constants,  etc.).  Kravtsov  [31]  and  Ludwig  [32]  have  analyzed 
the  field  near  the  caustic  surface  (smooth  envelope  of  a family 
of  ravs),and  have  developed  a "uniform  asymptotic  solution"  in  the  sense  that 
it  is  finite  at  the  caustic  and  reduces  to  geometrical  optics  awav  from  the 
caustic. 
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3.  GENERALIZATION  OF  GTD  AND  INVESTIGATION  OF  ALTERNATE  METHODS 
3. 1 Generalization  of  GTD  to  Arbitrary  Surfaces 

Keller's  generalization  of  GTD  for  the  analysis  of  the  field  diffracted 
from  a smooth  convex  object  is  closely  related  to  what  is  known  as  the 
"boundary  layer  technique"  in  the  theory  of  differencial  equations  [43]. 

On  the  other  hand,  the  "uniform  asymptotic  theory"  is  analogous  to  the 
method  used  by  R.  E.  Langer  [34]  and  F.  J.  Olver  [35]  to  find  the 
asymptotic  solutions  of  the  second-order  differential  equations  near  their 
"turning  points,"  which  are  counterparts  of  the  transition  regions  in  our 
case  [33],  [34]  and  [35]. 

The  second  prodecure  is  based  upon  the  generalization  of  the  G.O. 
interpretation  of  the  circular  cylinder  problem.  The  solution  obtained  by 
this  method  involves  some  functions  with  unknown  phase  and  amplitude,  similar 
to  Bessel  and  Hankel  functions.  Since  the  surface  of  a smooth  object  is 
actually  che  caustic  surface  of  diffracted  rays,  the  above-mentioned  formulacion 
is  applicable  in  this  case,  too.  Lewis  et  al.  [36]  have  modified  this 
solution  to  make  it  satisfy  the  boundary  condition  on  a convex  body.  Using 
ray  formalism,  they  have  obtained  an  asymptotic  solution  in  a complicated  form, 
which  they  call  "creeping  wave"  and  satisfies  the  boundary  condition  on  and  is 
uniformly  valid  near  and  away  from  the  surface.  It  should  be  mentioned  that 
the  method  has  been  developed  primarily  for  scalar  diffraction  problems. 

Creeping  waves  chat  travel  on  the  surface  of  the  body  generate  ocher 
kinds  of  diffracted  rays  in  the  presence  of  any  irregularities  in  the 
geometric  or  electromagnetic  characteristics  of  the  surface.  The  effects  of 
discontinuity  in  the  surface  curvature,  its  higher-order  derivatives,  or  the 
surface  impedance  have  been  studied  by  many  authors  [37],  [38],  [39]  and  [40], 
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An  exhaustive  study  of  various  diffraction  mechanisms  and  corresponding 
diffraction  coefficients,  and  constants  associated  with  the  propagation  of 
creeping  waves,  has  been  carried  out  by  Albertsen  [41], 

At  this  point,  let  us  examine  the  most  important  features  of  the  GTD 
and  its  various  modifications.  GTD  formulation  is  essentially  scalar  in  nature 
and  is  heuristic  in  some  parts.  Thus,  when  GTD  is  applied  to  a vector  problem, 
it  is  not  surprising  that  the  coupling  between  various  components  of  the  fields 
are  neglected,  and  each  one  of  them  is  treated  as  an  uncoupled  scalar  wave. 

The  other  assumptions  in  GTD  are  concerned  with  the  directions  of  these  field 
components  and  the  kind  of  boundary  conditions  they  satisfy  (see  Sec.  2.4). 

As  mentioned  earlier,  non-deductive  parts  of  GTD  are  based  on  asymptotic 
expansions  of  known  solutions  to  some  selected  "canonical'  problems.  Quite 
often  these  canonical  problems  are  not  general  enough  to  fully  and  accurately 
describe  the  local  behavior  of  the  field  for  an  arbitrary  structure.  Finally, 
most  of  the  canonical  problems  investigated  are  two-dimensional  in  nature.  The 
only  exception  to  this  is  the  sphere.  However,  insofar  as  the  geometric 
properties  of  the  surface  are  concerned,  the  sphere  is  a very  special  case 
since  its  radius  of  curvature  is  the  same  in  all  directions  and, consequently , 
the  surface  rays  are  torsionless.  Finally,  GTD  fails  when  the  observation 
point  is  located  in  the  transition  regions,  shadow  boundaries  or  in  the  neighbor- 
hood of  a caustic.  In  each  of  these  regions,  one  needs  to  carefully  modify  the 
GTD  formulas  and  often  such  a modification  is  not  too  simple.  Nevertheless, 
in  spite  of  these  difficulties,  GTD  is  recognized  to  be  a powerful  high-frequency 
technique  for  computing  the  leading  terms  of  the  asymptotic  solution.  Two  of 
the  principal  attributes  of  GTD  are  its  simplicitv  and  wide  scope  of  application. 


Li 
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3.2  Spectral  Domain  Approach 

We  now  examine  an  approach  different  from  GTD  which  uses  the  spectrum 

of  the  induced  current,  or  the  expression  for  the  radiated  field,  as  a starting 
point.  In  order  to  gain  a better  insight  into  the  curved-surface  radiation 

and  scattering  problem  and  to  verify  the  basic  assumptionsof  GTD,  it  is 

worthwhile  to  consider  such  alternative  approaches,  particularly  if  they  apply 

to  canonical  problems  which  are  more  general  in  nature  than  those  employed  to 

derive  the  GTD  results.  An  example  of  such  a study  would  be  to  consider  the 

case  of  surface  ray  propagation  with  non-zero  torsion,  a situation  that  occurs 

when  a magnetic  dipole  source  radiates  from  a location  on  the  surface  of  a 

circular  cylinder. 

The  geometry  of  the  problem  is  shown  in  Fig.  11.  The  radius  of  the 
cylinder  is  a and  the  source, which  is  an  infinitesimal  magnetic  dipole  with 
density  M,is  located  at  the  point  Q described  by  the  spherical  polar  coordinates 
(r*a,  3*90* , >«0°).  Each  point  P on  the  surface  of  the  cylinder  is  defined  by  a 
"geodetical  polar  coordinate"  system  (o,3),  where  3 is  the  arclength  of  the 
geodesic  connecting  Q to  P and  3 is  the  angle  between  $ (at  point  Q)  and 
geodesic  QP.  The  local  orthonormal  basis  vectors  (a, 8)  are  also  associated  with 
these  two  parameters.  The  observation  point  in  the  far  field  is  specified  by 
its  spherical  polar  coordinates  (r,8,$).  The  radiated  field  at  an  arbitrary 
point  can  be  expressed  in  terms  of  two  potentials , $ and  ?,  which,  in  cylindrical 
coordinates,  can  be  written  as: 
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For  the  problem  under  consideration,  we  can  express  the  spectral  weight 
coefficients  as 
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In  order  to  derive  an  asymptotic  expansion  of  (3.1)  and  (3.2),  we  proceed 
as  follows.  As  a first  step,  we  apply  Uatsonls  transformation  to  the  infinite 
summation  with  respect  to  n and  employ  appropriate  asymptotic  formulas  for 
Hankel  functions  with  large  order  and  argument  to  derive  the  following  expressions 
for  (3.1)  aii  (3.2)  under  the  conditions  that  ka  is  large  and  tp  small  compared  to  ir: 
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where 

ft  ■ k z + k fo+Ji($-n/2) ) 
z t 

m * (k^/2)1^3 
" m(4>-ir/2) 

ffl,  8q » 8^  ■ Fock's  functions  defined  in  Appendix  A. 


fl  and  M ■ components  of  M,  (M'n-O) 
P z 


Next,  applying  the  "saddle-point"  technique  to  (42)  and  (43)  and  keeping  only 

the  first-order  terms,  the  far  field  can  be  written  in  terms  of  its  components  along 

the  normal  and  tangent  to  the  surface  at  the  "stationary  point"  as 


where 
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P^:  is  Che  stationery  point  of  il  which  turns  out  to  be  the 

same  as  the  point  of  diffraction  predicted  by  GTD. 

. ,ka. 1/ 3 ...  . 1/3  . 

Cls  • (— ) (4>— w /2 ) • sin  0 

■ radius  of  curvature  of  geodesic  QP^ 
o * arc  length  QP^ 


R ■ the  distance  between  the  point  of  diffraction  P^  and 
the  observation  point 


n0  ■ 0,  x 0,;  normal  to  the  surface  at  P, 
- ^ ^ 1 


The  details  of  the  derivations  of  (3.8)  and  (3.9)  are  given  in  Appendix  B. 

Fig.  12  Illustrates  the  geometric  meaning  of  some  of  the  parameters 
appearing  in  (3.8)  and  (3.9),  for  the  observation  point  is  located  in  the 
shadow  region.  In  this  case,  , which  is  identical  to  £ given  in  (2.  lb'', 

is  the  reduced  distance  traveled  by  the  surface  ray  before  leaving  the 
surface  tangentially. 

In  the  lit  region,  the  geometric  interpretations  of  o and  £ are  shown 
in  Fig.  13.  The  rays, like  QP  T,  that  do  not  obey  the  generalized  Fermat's 
principle  are  called  "pseudo-ravs"  125].  The  ray  QP^P  appears  to  travel 
along  the  surface  up  to  the  point  P^  and  then  leaves  the  surface  at  tangen- 
tially in  the  opposite  direction,  to  reach  the  observation  point  P.  It  should  bo 
noted  that  formulas  (3.8)  and  (3.9)  give  us  the  contribution  of  the  ray  which 
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travels  along  the  shortest  path  on  the  surface,  and  thus,  suffers  the  least 
attenuation.  It  Is  not  difficult  to  see  that.  In  general,  there  are  Infinitely 
many  rays  which  contribute  to  the  total  field  at  any  observation  point.  However, 
their  contributions  are  very  small  compared  to  that  of  QP^  and  their  phases  and 
amplitudes  can  be  determined  in  a similar  manner  (Appendix  B) . Several  other 
remarks  on  the  formulas  given  in  (3.8)  and  (3.9)  are  in  order: 

a)  Numerical  results  indicate  that  good  agreement  between  (3.8), 

(3.9),  and  the  exact  modal  solution  is  obtained  for  ka  > 10. 

b)  The  zeroth-order  terms  in  the  asymptotic  expansion  of  the  normal 

component  of  the  field  E are  identical  to  those  given  bv  GTD; 

-1/3 

however,  the  k terms  derived  from  the  two  approaches  are 
different. 

c)  Tangential  component  of  the  field,  E , given  bv  (3.9),  also 

D 

is  different  from  the  corresponding  expression  based  on  GTD 
by  a multiplicative  polarization  factor.  Specifically, 

V<45>  • teSwa'TT1  <CTD>  <3-10' 

Consequently,  our  results  agree  with  GTD  only  for  the  circumferential 
ray,  i.e.,  for  8 ■ 0. In  addition,  for  an  axial  magnetic  dipole 

a » 

(M  • ^ ■ 0) , GTD  gives  a nonzero  value  for  the  field  in  the  3,  direction; 


our  solution  predicts  that  this  field  is  identically  zero,  a result 
which  is  in  complete  agreement  with  the  exact  solution  for  the 
problem. 


35 


Yet  another  important  distinguishing  feature  of  the  spectral  domain 
formulas  from  the  corresponding  ones  based  on  GTD  is  worth  noting.  As 
mentioned  earlier,  GTD  neglects  the  effect  of  torsion  on  the  diffraction 
of  surface  rays  propagating  along  smooth  objects.  In  contrast,  the  spectral 
domain  formulas  (3.8)  and  (3.9)  indeed  contain  the  torsion  effect  of  the 
surface  rays  as  may  be  seen  from  the  discussion  below. 

To  rearrange  the  formulas  (3.8)  and  (3.9)  in  a manner  such  that  the 
effect  of  the  torsion  on  the  surface  rays  is  explicitly  illustrated,  it  is 
necessary  to  examine  the  polarization  factors  in  the  above-mentioned  formulas. 
Referring  to  Fig.  12,  the  helical  path  of  a surface  ray  with  initial  angle  8 
can  be  parametrized  as: 

x(o)  » [a  cos  (q-c?3-) , a sin  (—-”*-) , o sing]  (3.11) 

The  torsion  of  the  ray,  t, is  found  from  the  relation. 
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and  n is  the  outward  unit  normal  to  the  surface  of  cylinder.  For  a helix 
on  a circular  cylinder,  the  torsion  is  constant  along  the  path  and  depends 
solely  on  the  initial  angle  of  the  curve.  The  torsion  r is  given  by  the 
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t “ 
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Combining  Che  above  relation  with  the  expression  for 
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one  obtains. 
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The  above  relationship  between  8,  T and  p^  will  be  found  useful  in  what  follows. 

Equation  (3.16)  enables  one  to  rewrite  the  source  polarization  factors 
in  (3.8)  and  (3.9)  in  the  following  form: 
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^*ie  exPlicit  manifestation  of  the  torsion  effect  can  be  seen  very 
clearly  if  we  use  (3.17)  and  (3.18)  in  (3.8)  and  (3.9).  The  final  result 
can  be  written  as 
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where  t is  the  torsion  of  the  surface  ray.  Note  that  the  effect  of  the 
torsion  on  the  surface  ravs  has  been  isolated  and  explicitly  expressed  in 
the  above  expression.  These  appear  as  correction  terms  to  GTD  and  are 

proportional  to  to  and  (to  ).  As  expected,  the  formulas  given  above  reduce 

j o 

to  the  conventional  GTD  results  for  the  limiting  case  of  no  torsion,  i.e., 
t - 0. 

In  contrast  to  GTD,  formulas  (3.8)  and  (3.9)  are  valid  irrespective 
of  the  location  of  the  observation  point,  be  it  in  the  lit,  shadow  or  transi- 
tion regions.  Although  not  valid  in  the  paraxial  region  (t$:90°),  they  can  be 
generalized  to  work  along  this  direction  also. 

The  solution  to  the  problem  of  radiation  of  a normal  electric  dipole 
can  be  found  from  the  expressions  (3.8)  and  (3.9)  for  the  radiated  far  field 
due  to  a magnetic  dipole  through  an  application  of  the  reciprocity  theorem. 

To  show  this,  we  apply  this  theorem  to  the  following  two  reciprocal  situations. 
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For  the  first  case,  let  (E,  H}  from  an  inf initesimal  electric  dipole  source 
of  moment  p,  located  far  from  the  cylinder  at  a point  P.  In  the  second 
situation,  let  {E™,  H*}  be  the  fields  radiated  by  an  infinitesimal  dipole 
of  moment  M,  located  at  a point  Q on  the  surface  of  the  cylinder.  As  the 
outer  boundary  surface  of  the  volume,  where  the  theorem  is  applied,  goes 
to  infinity,  the  "cross  flux"  through  it  vanishes.  Also  there  is  no  "cross 
flux"  through  the  inner  surface  (cylinder).  Therefore,  invoking  the  recipro- 
city relation,  one  can  derive 

- M • 2(Q)  - • E®(P)  (3.21) 

Sow  if  v«  apply  reciprocity  to  another  pair  of  cases,  in  which  the  first  one 
is  the  same  as  Che  first  situation  in  the  previous  example  and  the  second  one 
consists  of  the  radiation  {2*,  5*}  of  a normal  electric  dipole  of  moment  J, 
we  will  have 

J • E(Q)  - p • E*(P)  (3.22) 

Suppose  Che  magnitudes  of  J and  M are  unity,  then  it  follows  from 
(3.21)  and  (3.22) 


H^(Q)  - -p  • i“C  (P)  (3.23a) 

Hr(Q)  - -p  • 2“  (P)  (3.23b) 

Ep(Q)  - f • E*n  (P)  (3.23c) 


where  H (Q) , H (Q)  and  E (Q)  are  the  local  field  components  in  the  cylindrical 

<b  7.  0 

coordinate  system  coaxial  with  the  conducting  cylindrical  surface.  This  field 
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is  generated  by  p (Fig.  14).  The  following  set  of  relations  can  be  derived 

* 

from  (3.23)  by  replacing  p by  6^  and  respectively, 
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where  E^  (Q) , a - and  is  the  normal  electric  field  at  Q due  to  an  electric 

dipole  with  moment  p - 0 , located  at  P.  HJ  (Q)  and  Ha  (Q) , a " t»,,,  and  $ 

U u $ z 0 0 

are  the  tangential  components  of  the  magnetic  field  generated  by  an  electric 
dipole  of  moment  p ■ 0^,  located  at  P.  Invoking  the  Maxwell's  equations,  a 
relationship  can  be  sec  up  between  the  desired  far  field  E (P) , due  to  a 

V*o 

mc 

normal  electric  dipole  of  unit  moment,  and  the  known  solutions  E ’ (,P) 

(Equations  (3.3)  and  (3.9))  to  the  problem  of  the  radiation  of  a circumferential 

(mc)  and  axial  (ma)  magnetic  dipole  on  a cylinder.  As  a matter  of  fact,  E*  and 

0 

{H®,  H°)  c*n  be  related  by 
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where  * and  z are  cylindrical  coordinates  of  a point  on  the  surface  of  the 
cylinder.  Using  Equations  (3.25),  we  have 
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It  is  noted  that  the  derivatives  are  taken  with  respect  to  the  source  coordin- 
ates, whereas  Equations  (3. 8)  and  (3.9)  give  us  Che  field  components  generated 
by  a tangential  mangetic  dipole,  as  a function  of  the  observation  point  P 
coordinates.  However,  the  derivatives,  with  respect  to  the  source  coordinates, 
are  related  with  those  with  respect  to  the  observation  point  location  through 
a sign  reversal.  Therefore, 
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where  and  2q  are  the  observation  point  coordinates.  The  far-zone  radiated 
field  due  to  a normal  electric  dipole  of  unit  moment  located  on  the  surface 
of  a cylinder  can  be  found  by  evaluating  relations  (3.28)  using  (3.8)  and  (3.9) 
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Finally,  let  us  consider  the  possibility  of  the  generalization  of 
(3.8)  and  (3.9)  to  other  convex  surfaces  of  more  general  nature.  By  "more 
general  surfaces,"  we  mean  those  surfaces  which  are  not  substantially  different 
from  cylinders,  some  examples  being  cylinders  with  noncircular  (elliptical, 
hyperbolic,  and  parabolic)  cross  sections  and  conical  surfaces  with  small  apex 
angles.  The  key  step  in  a systematic  approach  to  generalizing  (3.8)  and  (3.9) 
is  to  use  the  generalized  definition  of  £ given  in  (2.16). 

Fig.  15  exhibits  some  initial  results  of  the  generalization  of  these 
formulas  to  the  case  of  a cone.  In  these  graphs,  the  Hughes  results  have  been 
reproduced  from  [63].  It  is  evident  that  results  obtained  from  the  present 
approach  agree  quite  well  with  the  series  solution  which  is  rather  tedious  and 
time-consuming.  We  also  observe  from  Fig.  16  that  there  is  a noticeable 
discrepancy  between  the  analytical  solution  and  the  experiment.  Thus,  within 
the  range  of  experimental  error,  our  results  agree  quite  well  with  those 
published  in  the  literature. 

Before  concluding  this  subsection,  it  is  worth  mentioning  some  basic 
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points  concerning  the  applicability  of  this  formulation  to  the  case  of  lossy 
surfaces  which  may  be  approximately  described  in  terms  of  surface- impedances. 
For  this  case,  the  simplifying  assumption  of  infinite  conductivity  for  the 
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Figure  15  Comparison  beeveen  spectral  domain  results  (UI) , 

modal  approach  v Hughes) , and  experimental  measure- 
ments .Hughes).  Che  UI  results  are  ierivec  iron  a 
generalized  version  of  (3.3)  and  (3.9)  for  a cone. 
The  Hughes  resales  nave  beer,  reproduced  from  [63] 
and  are  based  on  a modal  series  of  13  earns.  All 
results  are  for  \,'Z  radial  sloe  on  a cone  of 
half-angle  10  . ( i^OO®  cut) 
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scatterer  no  longer  holds, the  problems  become  tremendously  complicated, 
especially  when  the  surface  of  the  scatterer  Is  curved.  However,  for  special 
cases  where  the  frequency  Is  high  and  the  geometry  of  the  body  satisfies  some 
specific  conditions,  certain  approximate  procedures  are  available  for  solving 
the  problem  [64],  [651,  [66].  Some  of  these  approaches  are  based  on  the 
modification  of  the  exact  boundary  conditions  in  a manner  such  that  the  solu- 
tion to  the  imperfectly  conducting  body  scattering  problem  can  be  obtained  as 
a perturbation  of  the  existing  solutions  for  the  ideal  case  of  a perfectly 
conducting  surface.  These  approaches  are  mostly  based  on  the  use  of  the 
Leontovich  boundary  conditions  [64],  They  relate  the  tangential  components 
of  the  electric  and  magnetic  fields  (or  the  normal  components  and  their  normal 
derivatives)  via  a surface  impedance  condition  which  is  a function  onlv  of  the 
electromagnetic  properties  of  the  material.  Thev  are  valid  for  surfaces  whose 
radii  of  curvature  are  large  compared  with  the  penetration  depth,  and  also  for 
materials  which  are  not  homogeneous,  but  whose  properties  vary  slowlv  from 
point  to  point.  As  the  refractive  index  (or  conductivity)  increases  to 
infinltv,  the  conditions  go  over  unlformlv  to  those  for  perfect  conductivity. 

The  Leontovich  boundary  conditions  can  be  expressed  as 

i-(n'E)a«nZn*it  (3.30) 

A 

where  n is  the  unit  normal  to  the  surface,  1 is  Che  impedance  of  free  space  and 
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a,  u,  e ■ electromagnetic  characteristic  of  the  body 
Z ■ Ai/e  - free  space  impedance 
N - complex  refractive  index  of  the  body. 


The  conditions  for  applicability  of  (3.30)  to  the  curved  surfaced  can  be 


formulated  as 

|n|  » 1 


Im  N ko  » 1 


(3.31a) 


(3.31b) 


where  p is  the  smallest  radius  of  curvature  or  dimension  of  the  body  [64]. 

When  n is  very  small, we  can  expand  the  fields  {E,  H}  around  their 
values  {Eg,  3g}  for  the  case  of  perfectly  conducting  surfaces  (n_0),  in  the 
ascending  powers  of  n as 


E - fQ  + + n2  E2  + . . . 


H - HQ  + + r>2H2  + ... 


Inserting  (3.32)  in  (3.30)  one  finds 


(3.32a) 


(3.32b) 


Eg  - (n  • Eg)  n - 0 


(n  • E^)  n ■ nZ  n x Hg 


(n  • Ej)  n • nZ  n * H^, 


(3.33a) 


(3.33b) 


(3.33c) 
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Equation  3.33b  is  a first-order  approximation  for  the  boundary  condition  for 
the  tangential  electric  field  on  the  surface  of  the  body.  The  function  HQ 
is  the  known  solution  for  a purely  conducting  surface.  Using  boundary  condition 
(3.33b),  one  is  able  to  find  the  first-order  solution  for  the  electric  field 
in  the  presence  of  imperfectly  conducting  homogeneous  or  inhomogeneous  objects. 
The  higher-order  solutions  can  be  found  in  a similar  manner. 

3. 3 Approach  Based  on  an  Asymptotic  Evaluation  of  the  Radiation 
Integral  of  the  Surface  Current 

As  a final  topic  in  this  chapter,  we  consider  an  approach  based  on  the 
asymptotic  evaluation  of  the  radiation  integral  expressed  in  terms  of  the 
induced  surface  current  which  is  itself  derived  in  an  asymptotic  manner  for 
surfaces  with  large  radius  of  curvature. 

It  was  shown  in  Chapter  2 that  Fock's  theory  can  provide  us  with  an 
expression  for  the  scattered  field  in  the  neighborhood  of  a smooth  convex 
body  illuminated  by  a plane  wave.  Using  this  solution  in  conjunction  with 
the  reciprocity  principle,  we  can  find  the  far  field  radiated  bv  a point 
source  located  on  the  surface  of  the  body.  By  generalizing  the  definition 
of  £ in  Fock's  theory,  we  can  also  write  the  final  result  in  a GTD 
format  and  represent  it  as  a surface  ray.  The  total  field  at  a point  on 
the  surface  is  obtained  by  adding  all  the  possible  rays  which  reach  the 
observation  point  P.  Various  techniques  can  be  used  to  determine  the  field 
propagation  along  these  rays.  For  instance,  when  the  source  is  located  on 
the  surface,  and  the  surface  is  a conical  one,  the  field  at  each  point 


can  be  decomposed  into  two  parts 
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F - Fl  + F2  (3.34) 

where  F^  is  the  geometrical  optics  field  when  the  observation  point  is  directly 
illuminated  by  the  source,  and  is  the  creeping-wave  contribution  derivable 
via  an  extension  of  Fock's  theory  when  the  point  is  in  the  shadow  region.  The 
other  term,  F^,  is  the  so-called  tip  contribution,  and  can  be  obtained  by 
physical  optics  or  GTD.  Goodrich  et  al.  [42]  have  applied  this  procedure  to 
find  the  radiation  pattern  of  slot  arrays  on  cones. 

The  approximate  induced  surface  current  distribution  can  be  obtained  by 

Fock's  theory,  GTD  [13],  [14],  [16],  and  [25]  or  some  other  appropriate 
high  frequency  technique.  The  induced  surface  current  due  to  a magnetic 
dipole  on  a perfectly  conducting  circular  cylinder  and  cone  has  been  calculated 

by  Chang,  et  al.  [44],  and  Chan,  et  al.  [45],  whose  procedure  is  based  upon 
an  asymptotic  expansion  of  the  exact  modal  solution  to  the  above-mentioned 
problems.  Lee,  et  al.  [46]  and  [22]  have  treated  the  same  problem  by  a method 
based  on  Fock's  asymptotic  solution  of  the  problem  of  a sphere  [47].  These 
expressions  for  the  current  distribution  can  be  used  in  the  radiation 
integral  representation  of  the  far  field. 

The  numerical  evaluation  of  this  integral  is  a formidable  task,  especiallv 
when  the  frequency  is  very  high.  Thus,  it  is  highly  desirable  to  have  an 
analytical  and  explicit  formula  for  the  far  field  expressed  in  terms  of  the 
surface  current.  We  now  discuss  an  approach  for  accomplishing  this  task  and 
examine  the  problem  of  deriving  an  asymptotic  expansion  of  the  far  field 
radiated  due  to  a point  source  located  on  the  surface  of  a smooth,  conducting, 
and  convex  body  of  an  arbitrary  shape. 


50 


Consider  an  arbitrary  smooth  convex  surface  S shown  in  Fig.  17.  Let 
a magnetic  dipole  source  be  located  at  a point  Q on  S.  We  parametrize  the 
surface  S introducing  a "geodetical  polar  coordinate"  system  with  the  pole 
located  at  Q such  that  an  arbitrary  point  on  the  surface  is  defined  by  a 
pair  of  numbers  (a, 8),  where  o is  the  arc  length  of  the  geodesic  QP^  and  8 
is  the  angle  between  QP^  and  some  reference  direction  at  Q.  Unit  vectors 
along  the  constant  parameter  curves  o and  6 are  locally  orthogonal.  The 
unit  normal  to  the  surface,  n,  is  given  by  n ■ o « j.  An  element  of  length 
in  this  coordinate  system  may  be  written  as 

ds2  - do2  + G dd2  (3.35) 


The  radiation  integral  for  the  scattered  far  field  can  be  written  as 


JU-RK) 


exp(-jkR) 

R 


dS 


(3.36) 


where  R is  the  distance  between  any  point  on  the  surface  and  the  observation 
point.  In  Che  geodetical  polar  coordinate  system,  we  can  rewrite  a scalar 
component  of  (3.36),  say  M,  in  terms  of  a double  integral  of  the  following 
general  form 


M - / / F(o,8;P)  exPl-lk.(R+?.l1  v£  d0  dt5 

D K 


(3.37) 


where  we  have  assumed  the  following  form  for  the  surface  current: 
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J(P.)  ■ ?(P  ) exp(-jk  o)  - J o +•  J B 
l 11  a tJ 

k^  ■ k - jkj*  k - 2*/\,  k^  > 0,  k^  <<  k 


(3.38) 


where  f(P^)  is  a relatively  slowly  varying  function  when  k,  the  free -space 
wave  number,  is  large.  This  assumption  is  based  upon  a close  scrutiny  of 
different  asymptotic  formulas  given  for  the  induced  surface  current. 

’-/hen  the  observation  point  P is  located  in  the  shadow  reeion,  the 

main  contribution  to  (3.37)  comes  from  a small  neighborhood  of  the  stationary 
point  of  the  integrand,  and  the  stationary  phase  method  for  multiple  integrals 

([48]  and  [49])  is  applicable.  The  asymptotic  expansion  of  (3.36)  has  been 

-5/3 

derived  up  to  the  order  of  k . The  details  of  calculation  can  be  found 
in  Appendix  C.  The  final  result  is 


+ Vl  + 0(k~11/6)) 


(3.39) 


where 


[(AJ.  + -rr  (J,ejkj)e~Jk,:7)3  + (BJ,  + CJ  ) n] 

3 a a 0 J 


(3.40) 
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(3.41) 
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o • r ( 2/ 3)  • r(i/2)  . P 
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(3.43) 


A,  B,  and  C are  dependent  upon  geometric  properties  of  the  surface  at  the 
stationary  point  which  turns  out  to  be  exactly  the  same  as  the  "point  of 
diffraction"  of  surface  rays.  The  quantities  A,  B,  and  C are  given  by 


A 


2G 


0 i W ’ t * 

— • — (— ) + -£•[- 
2 3o  vp  ’ 2G  lo 
o o 


SB 


+ (L30)2  - (1/2^-t]  + 0 (i) 

3o“ 


B - L3a/G1/2, 


C - - I/O. 


(3.44) 


(3.45) 


where 

o • radius  of  curvature  of  the  geodesic 
o - geodetic  radius  of  curvature 

J J - 

L , L * coefficients  of  the  second  fundamental  form  of  the  surface  (S)  . 


A geometric  interpretation  of  these  parameters  has  been  illustrated 

1/2 

in  Fig.  18.  It  is  evident  from  this  figure  that  ~1  is  simply  the 

8 

divergence  factor  of  the  rays  leaving  the  surface  tangentially  at  the  point 
of  diffraction.  In  using  formula  (3.43),  we  should  bear  in  mind  that  the 
various  terms  in  Un  and  are  not  of  the  same  order.  For  example,  in  the 
deep  shadow,  J is  exponentiallv  larger  than  J_,. 
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The  formulas  given  in  (3.43)  have  been  tested  and  compared  with  other 
available  solutions.  An  important  conclusion  derived  from  this  comparison 
is  that  although  the  method  of  radiation  integral  is  based  on  less  restrictive 
assumptions,  it  is  perhaps  not  as  useful  as  the  spectral  domain  approach 
because  the  stationary  point  of  the  phase  of  the  integrand  in  (3.37)  is  of  the 
second  order,  and  hence,  the  asymptotic  expansion  of  this  integral  converges 
rather  slowly  except  when  ko ^ is  very  large  (;40  or  more). 
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4.  SOURCE  RADIATION  IN  THE  PRESENCE  OF  AM  INFINITELY  CONDUCTING  CIRCULAR 
CYLINDER  OF  FINITE  LENGTH  (GENERAL  FORMULATION  OF  THE  PROBLEM) 

An  asymptotic  analysis  of  the  problem  of  radiation  of  an  electro- 
magnetic point  source  located  on  a solid  metallic  cylinder  of  finite  length 
is  the  topic  which  is  studied  in  detail  in  the  remaining  chapters.  Because 
of  its  theoretical  and  practical  Importance  (small  antennas  mounted  on  the 
body  of  the  satellites,  conformal  arrays,  electromagnetic  coupling  through 
the  slots  on  the  surface  of  a cylinder,  and  etc.),  the  problem  of  radiation 
in  the  presence  of  or  the  diffraction  by  the  finite  cylinder  has  been  treated 
by  many  investigators.  Various  numerical  and  asymptotic  techniques  have 
limited  application  since  they  are  usually  valid  only  for  a special  range  of 
parameters,  e.g.,  the  length  or  radius  of  the  cylinder,  the  location  of  the 
source,  or  a special  type  of  incident  field.  When  the  length  and  the  radius 
of  the  cylinder  are  small,  compared  to  the  wavelength  (fraction  of  a wavelength), 
the  2-D  integral  equation  for  the  surface  magnetic  field  or  the  two  one- 
dimensional (generally  coupled)  integral  equations  for  the  transverse  and 
axial  components  of  the  induced  surface  current  are  solved  by  the  moment  method 
or  point-matching  techniques  [67],  [63],  [69],  [70],  and  [71].  Other  approaches 
have  been  presented  for  the  cases  where  the  length  of  the  cylinder  is  very 
large  and  its  radius  is  relatively  small  or  the  incident  field  has  some  special 
characteristics  [72],  [73],  [74],  [75]  and  [76]. 

When  the  frequency  is  high  or  the  length  and  the  radius  of  the  cylinder 
are  large  compared  to  the  wavelength,  several  asymptotic  procedures,  mostly 
based  on  GTD  and  its  modification,  have  been  designed  and  applied  to  this 
problem  [77],  [78],  [79],  [30],  and  [81]  with  various  degrees  of  success. 

The  asymptotic  method  described  in  this  work  is  essentiallv  a spectral 


domain  approach  and  is  very  closely  related  to  the  Spectral  Theorv  of  Diffrac- 
tion (STD)  developed  by  R.  Mittr3  et  al.  [82],  [38],  and  [39]. 
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The  present  approach  is  presumably  more  accurate  than  the  other  avail- 
able methods  because  it  avoids  most  of  the  restrictive  assumptions  underlying 
the  conventional  high  frequency  techniques.  In  this  chapter,  we  merely 
outline  the  basic  steps  to  this  approach.  The  details  are  presented  in  the 
later  chapters. 

In  the  high  frequency  range  where  the  ray  formalism  is  the  most  con- 
venient representation  for  the  Maxwell's  equations,  the  field  at  each  point 
of  the  space  can  be  expressed  as  the  sum  of  the  contributions  of  different 
rays  passing  through  that  point.  Various  diffraction  or  reflection  mechanisms 
give  rise  to  different  types  of  reflected  or  diffracted  rays  which  together 

I 

with  the  direct  rav  constitute  the  total  field.  In  the  present  problem, 

illustrated  in  Fig.  19,  the  effect  of  the  finite  solid  cylinder  on  the  radiated 

field  at  P due  to  the  source  Q (located  on  the  surface)  has  been  represented 

by  three  types  of  rays.  F is  the  surface  diffracted  ray  (or  psuedoray,  if 

the  P is  in  the  lit  region  of  Q)  corresponding  to  the  radiated  field  due  to 

wl  w2 

current  induced  on  the  surface  of  the  cylinder.  F and  F “ are  wedge  diffracted 

rays  generated  by  the  curved  wedges  at  both  ends  of  the  cylinder.  To  calculate 

Fs,  if  the  source  is  far  from  both  ends,  the  infinite  cylinder  formula  derived 

in  Chapter  3 or  in  [85]  can  be  used.  These  formulas  are  discussed  in  Chapter 

wl  w2 

5.  Wedge  diffracted  fields  F and  F can  be  analyzed  only  when  the  surface 

field  or  current  are  available.  Chapters  6 and  7 are  devoted  to  computation 
wl  w2 

of  F and  F . In  Chapter  8,  by  calculating  the  effect  of  truncation,  we 
remove  the  assumption  made  in  Chapter  5 concerning  the  length  of  the  cylinder 
and  extend  the  range  of  applicability  of  the  procedure  to  the  cases  where  the 
source  is  not  too  far  from  the  upper  or  lower  wedge  of  the  cvlinder.  Chapter 
9 is  devoted  to  the  calculation  of  the  field  in  the  paraxial  region.  A 
summary  of  numerical  results  and  conclusions  are  given  in  Chapters  10  and  11. 
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5.  SURFACE  DIFFRACTED  RAYS 

Let  as  consider  Che  case  where  an  infinitesimal  magnetic  dipole  with 
the  moment  M placed  on  the  surface  at  a location  described  by  (p  * a, 

0*0,  z * 0) . The  orientation  of  the  dipole  is  tangent  to  the  surface  and 
the  direction  of  observation  is  given  by  ^ the  ob9ervation  point  is 

not  placed  in  the  lit  region,  no  direct  ray  can  reach  it.  In  this  case 
only  the  surface  diffracted  rays  that  travel  along  the  curved  surface 
can  pass  through  the  observation  point  after  leaving  the  surface  tangentially. 
In  general,  there  are  an  Infinite  number  of  rays  which  contribute 
to  the  field  at  P.  In  most  of  the  cases,  however,  only  two  of  these  rays 
traveling  in  opposite  directions  around  the  cylinder  are  taken  into 
consideration.  If  the  radius  of  the  cylinder  is  large,  the  others  suffer 
strong  attenuation  and  thus  have  negligible  effect  on  the  total  field  at  P. 

Propagation  of  the  field  along  these  rays  is  governed  by  formulas 
(3.8)  and  (3.9): 


(M  • 8.)  • (J-rr ) • gn(SiC) 


-jka 
• e J 


' Air  7 80  IS 

kp  >1/3 

T ‘ 81(HS)  * 


jkR  (M  • 01)(01  • Sj) 
R " Aita 


(5.1) 


(M  * »1)(*1  ' V . f^al2/3  . -jka 
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(5.2) 


The  far-zone  fields  along  rays  no.  1 and  no.  2 should  be  calculated 

separately.  The  contribution  of  surface  diffracted  rays  is  the  sum  of  these 

two  fields.  For  circular  cylinder  geometry,  simple  relationships  exist 

between  (J,3)  and  (6  ,0  ) ; (Figs. 20  and  21). 

o o 

Rav  no.  1: 


5 * t/2  - 8 


(5. 3a) 
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a * a(<p  - 7r/2)/cos  3 
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(5.3b) 

(5.3c) 

(5.4a) 

(5.4b) 

(5.4c) 


The  field  corresponding  to  ray  no.  1 is  obtained  by  inserting  (5.3a,b,c) 
into  (5.1)  and  (5.2).  In  order  to  compute  the  field  along  ray  no.  2,  the 
minus  sign  between  the  two  terms  in  (5.1)  should  be  changed  to  a plus  sign  and 
the  parameters  in  (5.1)  and  (5.2)  are  then  replaced  by  their  values  given  in 
(5.4a,b,c) . 

Throughout  this  chapter  and  later  on,  the  source  Q is  chosen  to  be 
the  phase  reference  point.  Consequently,  R must  be  replaced  by 

R - r - QF>1  • r (5.5) 

or 


R = r - (a  cos“  0q  - a sin  Sq  cos  4>q)  (5.6) 

where  r is  the  distance  between  the  source  and  the  observation  point. 

A -N 

The  components  of  the  far-zone  field  along  0q  and  4^  can  be  found  in 
terms  of  its  components  along  §-  and  n_  very  easily  ; 

E..  * -E,  (Rav  no.  1 + Rav  no.  2)  (5.7) 

0o  32 

E * -E  (Ray  no.  1 + Ray  no.  2)  (5.8) 


The  case  in  which  the  point  source  is  a normal  electric  dipole  is 
treated  in  a similar  manner.  In  this  case,  the  expressions  (3.29a)  and 
(3.29b)  of  Chapter  3 should  be  used  instead  of  (3.8)  and  (3.9).  The 
remaining  steps  remain  exactly  the  same. 
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6.  SURFACE  FIELD  DUE  TO  A POINT  SOURCE  ON  A CIRCULAR  CYLINDER 
The  cylinder  problem  has  an  exact  modal  solution,  which  is  in  the  form  of 
an  infinite  series  with  each  term  containing  an  infinite  integral  [86],  [87], 
For  a large  ka,  this  solution  is  very  slowly  convergent  and  becomes  less 
useful.  Among  several  asymptotic  solutions,  those  given  by  Hwang, 

Kouyoumj ian , Wang,  and  Pathak  [88],  [89],  [90]  and  by  Chang, 

Felsen,  Hessel  and  Shmoys  [91],  [44]  are  approximately  deduced  from 
Che  exact  modal  solution  under  the  condition  ka  ■+ ».  In  this  work  we  use 
the  asymptotic  solution  obtained  by  Lee  and  Safavi-Naini  [21],  [22]  which 
gives  the  surface  magnetic  field  everywhere  from  the  source  point  to  the 
deep  shadow  more  accurately  and  is  based  on  a classical  work  by  Fock  in  1949 
on  spheres  [65].  A nondeductive  part  of  this  solution  was  later  justified 
theoretically  by  Boersma  and  Lee  [92]. 

Let  us  first  consider  the  case  of  a tangential  magnetic  source.  The 
normal  electric  source  is  discussed  later.  At  a point  Q on  the  surface  of 
the  cylinder  (Fig.  22),  there  is  a tangential  magnetic  dipole  source  described 
by  a magnetic  dipole  moment  M.  The  problem  is  to  determine  H at  another 
point  P(a,4>,  z)  on  the  surface  under  the  assumption  that  ka  is  large. 

At  this  point  we  introduce  several  paramters.  According  to  GTD  [23], 
[24],  [25],  the  dominant  contribution  of  H at  P is  the  field  on  the  surface 
ray  which  is  a geodesic  on  the  conducting  surface,  and  in  the  present  case 
is  a helical  path  (Fig.  22).  The  arclength  of  the  surface  ray  is 

a ■ *(aor"+V  <«•» 


The  surface  of  the  cylinder  is  parameterized  by  geodetical  polar  coordinates 
(a,3)  discussed  in  the  previous  chapter.  The  tangent,  normal  and  binormal 


mmmm 


66 


of  the  surface  ray  are  (o^.-n^.S^)  at  Q and  (op,-np,Sp)  at  P.  At  any  point 
on  the  surface  ray,  the  curvature  of  the  conducting  surface  is  described 
by  two  parameters: 

(3,  * the  radius  of  curvature  in  the  direction  of  S (or  that  in  Che 

u 

longitudinal  direction  of  the  surface  ray) 

Pg  * the  radius  of  curvature  in  the  direction  of  S (or  that  in  the 
transverse  direction  of  the  surface  ray) 

On  a convex  surface,  both  p and  p are  non-negative.  For  the  present  case 

o p 

of  a conducting  cylinder,  one  has 


cos* 


B sin;  g 


(6.2) 


where  8 is  measured  from  the  a<J>-axis  in  Fig.  22  and  takes  a value  between  0 
and  2tt.  The  large  parameter  for  our  asymptotic  expansion  is 


m 


(6.3) 


Thus,  our  solution  is  an  approximate  (not  rigorous)  asymptotic  solution  valid 
for  m ■+  »,  up  to  and  including  terms  of  0(m  ).  Furthermore,  let  us  introduce 

a distance  parameter 


5 


mI/3 

UpJ 


kg 

2m‘ 


(6.4) 


which  is  the  same  as  that  defined  in  the  previous  chapter.  Note  that  £ » 0 
defines  the  lit  region  (8  * t/2)  , < 1 defines  the  penumbra  region,  and 

r >>  i defines  the  deep  shadow.  Our  solution  is  uniformly  valid  for  all  £ 0. 

Due  to  the  point  source  M at  Q,  our  final  asymptocie  solution  for  the 
magnetic  field  on  the  surface  is  given  by 


H(P) 


(iQBpH3  + JQ*pHa) 


(6.5a) 
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where  Che  transverse  component  is 


K,(P)  - 
D 


1 " to 


r(5)  - [~]  u(5>  + j(^  kPo)‘:/3  • v’(0 


+ j (*2  kpa)"2/3  * (0a/0g)  * u'(5)jG(a), 


(6.5b) 


the  longitudinal  component  is 
H,(P)  - 

and  the  function  G(j)  is 


£)£«>  + (l  - * J<>*  k^)-2/Iu’«)]ow) 


(6.5c) 


G(a) 


k‘Y  -jko 
o e 


2irj  to 

I/O  _1  , 

Here  Y * (e  /u  ) * (120tt)  , v and  u are  defined  in  Appendix  A,  and  v 


(6 . 5d) 


o o o 

is  the  derivative  of  v.  It  can  be  shown  that  if  ka  tends  to  infinity  the 


exact  solution  for  the  ground  flat  plane  will  be  recovered.  When  3 — t/2 
(paraxial  region)  but  the  radius  of  the  cylinder  remains  finite,  we  have 


Vp>  - <h8<MW«  + B "t'1 


-jko 


'Tco 


(6.o) 


as  ko  approaches  infinity.  Another  limiting  case  of  special  importance 
to  us  is  the  case  where 


y oo 


(6.7) 


which  occurs  when  observation  point  P is  in  the  deep  shadow.  Making  use  of 
(A-20)  through  (A-23) , we  have  from  (6.5)  that 


H , (P)  - - 

1528  (ka) 


2 2/3  , 

k cos  a 


5t 


■ l7,  * exp* -0.385  - j fy  + 0.515  + ko 
1/3(kc)L/"  - 


(6.8a) 


H (P)  ' r*-  H,(P) 
•J  KJ  3 


(o.8b) 


Therefore,  in  the  deep  shadow,  the  field  is  a slow  wave  and  decays  exponentially 


along  the  surface  ray.  In  later  applications,  explicit  field  expressions  are 


needed  for  axial  and  circumferential  dipoles.  They  are  listed  below 


Axial  dipole;  M » z 


Circumferential  dipole:  M 


One  can  now  derive  similar  expressions  for  the  electromagnetic  surface 


fields  of  a radial  electric  current  J on  a perfectly  conducting  circular 


cylinder.  Let  us  define  J as 


This  derivation  has  not  been  given  in 


Here,  the  main  steps  are  discussed 


briefly.  It  is  noted  that  p is  the  dipole  moment  which  is  oriented  normal  to 


the  surface  at  (a,  0,0).  The  tangential  magnetic  field  H (P)  due  to  J can  be 


obtained  by  applying  the  reciprocity  theorem  to  the  fields  of  M and  p as  follows 


?n(P) 


?m(P) 


is  the  radial  component  of  the  surface  electric  field  atQdueoothe 


tangential  magnetic  dipole  M located  at  P,and  H 


is  the  surface  magneti 


field  at  P due  to  the  radial  electric  dipole  of  moment  p located  at  Q 

.ra(P)  j . _ . _ , ^m(P)  I ... 


s equations 
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where  H 


mU?) 


and  H 


m(P) 


are  given  bv  (6.9)  and  (6.10)  when  che  source  is 

‘ »Q  * 'Q 

located  at  P and  the  field  is  calculated  at  Q (Fig.  23).  The  axial  and 

circumferential  components  of  the  surface  magnetic  field  H®^^,  H®^3  due 

2 

to  a radial  electric  point  source  now  can  be  computed  by  (6.12)  and  (6.13) 
as  follows: 

-*e(0)  -*■ 

Axial  component  due  to  p at  Q: 


^(Q) 

Z 


■+  . ?m(p) 

. (~P  * P) 

3H  v J 
z 

» -p  . E 

P p 

tq.M-a  ja)£a 

3<p ' 

Q,M«2 


Circumferential  component  due  to  p at  Q: 

P 


a*(Q) 


^(Q) 


- . sm(P) 

-p  • E 


r 3H*(P) 

(P  * . _$ 


jwe 


3z' 


(6.14) 


(6.15) 


It  is  noted  that  the  right-hand  sides  of  the  above  equations  are  derivatives 
of  the  surface  field  due  to  a magnetic  dipole  M located  at  P,  whereas  the 
left-hand  sides  are  the  field  quantities  at  P due  to  the  source  p at  Q.  Surface 
derivatives  of  H^P^  can  be  written  in  terms  of  the  derivatives  of  tangential 
and  binormal  components  of  H™^P^  with  respect  to  3 and  J: 


3Hm(P) 

2 
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3H^P) 

£ 
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! 0 

36’  3HS’ 
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) ■> 

| sin"  S' 

36' 

3c  M 

3z' 

( 33' 

3z ' 3o' 

3z’ 
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32' 

36’ 

* 3 2 ' j 
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* cos“ 

S’  + (H,,  - 
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0 

36' 
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(6 

In  deriving  the  relations  (6.16)  and  (6.17),  the  Jacobian  of  transformation 


between  the  two  coordinate  systems  (o', 8')  and  (<{> 


z')  is  needed.  Below 


electric  diple  p are  obtained  by  substituting  (6.16)  and  (6.17)  back 
into  (6.14)  and  (6.15).  The  relationship  between  the  (a', 8')  coordinate 


system  with  the  pole  at  P and  the  (o, 3)  coordinate  system  with  pole  Q is 
illustrated  in  Fig.  23: 
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7.  DIFFRACTION  OF  SURFACE  RAYS  BY  CURVED  WEDGE 

As  discussed  previously,  Che  computation  of  the  radiated  field  for 
high-frequency  electromagnetic  point  sources  located  on  a finite  solid 
cylinder  it  must  include  diffracted  rays  generated  by  the  surface  ravs 
reaching  the  edges  of  the  truncated  cylinder.  Diffraction  of  surface 
rays  (creeping  waves)  has  been  investigated  by  a number  of  workers; 
however,  diffraction  coefficients  for  these  creeping  waves  are  not 
readily  available  in  a convenient  form.  Hence,  in  this  work,  we  address 
ourselves  to  the  problem  of  deriving  the  diffraction  of  a creeping  wave  by 
a wedge  and  use  it  subsequently  to  compute  the  edge-contribution  to  the 
radiated  fields. 

The  difficulty  with  the  analysis  of  the  diffraction  process  of  surface 
rays  at  wedges  stems  from  the  complexity  of  the  associated  canonical  problem. 
It  is  well-known  that  the  solution  of  a suitably  chosen  canonical  problem 
is  a first  step  toward  deriving  the  diffraction  coefficient  for  a given 
problem  with  a complex  geometry.  The  choice  of  the  canonical  problem  is 
based  on  two  principal  considerations: 

1)  The  geometry  of  the  canonical  problem  should  contain  all  the 
important  features  of  the  local  geometry  of  the  original  problem. 

2)  The  solution  of  the  canonical  problem  must  be  expressible  in  a 
convenient,  usable  form. 

As  one  example  of  such  a geometry,  Kouyoumjian  and  Burnside  [93] 
considered  the  case  of  a cylinder-tipped  half-plane.  The  case  of  ogive  was 
treated  by  Peters  and  Ryan  [94]  in  an  empirical  manner.  Ryan  [951  and  Knott 
et  al.  [96]  studied  the  finite  cylinder'.  They  assumed  a wedge  diffracted  ray 
excited  by  creeping  waves  hitting  the  wedges  of  the  cylinder  could  be 


73 


constructed  by  using  equivalent  currents,  determined  from  the  surface  field 
of  the  waves  and  the  wedge  diffraction  coefficients  for  a straight  wedge. 

This  approach  was  consistent  with  the  results  found  by  Albertsen  and 
Christiansen  [97],  [41].  A somewhat  similar  procedure  is  applied  to  find 
the  wedge  diffracted  fields. 

7 . 1 Diffraction  Points 

For  the  case  of  the  finite  cylinder  studied  in  this  work,  the  surface 

rays  emanate  from  the  source  Q (Fig.  19)  and  propagate  along  all  directions 

on  the  surface  of  the  cylinder  along  a helical  path  (geodesics) . Surface 

rays  creep  over  the  surface  until  the  curved  wedges  are  reached.  One  of 

these  rays  has  been  shown  in  Fig.  24.  After  diffraction  on  the  upper 

wedge  (wedge  no.  2)  at  diffraction  point  , a cone  of  diffracted  rays  with 

the  axis  t (unit  tangent  to  the  wedge  at  P .,)  and  the  half-angle  3 is 

w - 

formed.  3 , is  also  the  angle  between  the  surface  rav  and  circumferential 
w 2 

• < 

direction  at  Q.  a indicates  the  direction  of  incidence  of  the  surface  ray 
on  the  wedge  and  is  tangent  to  the  rav  at  P The  angle  41  , counted  from 
the  x-axis  determines  the  location  of  the  diffraction  point  on  the  wedge, 
n - is  the  radial  unit  vector  normal  to  the  wedge  lving  in  the  plane  of  the 
upper  end  cap.  o is  the  unit  vector  along  a diffracted  ray.  0 and  0 both 
lie  on  the  diffraction  cone.  3^-,  and  are  dependent  upon  the  location  of 
diffraction  point  P^ . Therefore,  each  point  of  the  wedge  is  assigned  a cone 
of  diffracted  rays  with  a specified  semiangle.  The  ensemble  of  directions 
defined  by  these  cones  covers  the  whole  space.  Conversely,  associated  with 

each  direction  (observation  direction)  in  the  space,  t.iere  is  1 finite  number 
>(  rays  originating  from  specific  points  on  the  wedge,  which  ontribut,  to 

the  far-zone  field  in  that  direction.  Hence,  the  first  step  in  computing 
the  contribution  of  wedge-diffracted  rays  to  the  total  radiated  field  in  a 
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given  direction  is  to  determine  the  location  of  the  diffraction  points  on 
the  wedge  corresponding  to  that  special  direction.  It  is  obvious  that 
this  location  is  solely  determined  by  the  value  of  Pw2'  Now  let  us  derive 
the  relationship  between  ^w>  and  the  observation  direction  defined  by  the 
observation  angles  ( 0q , 4>0 ) * 

To  derive  the  above-mentioned  relationship,  it  is  noted  that  the 
direction  of  incidence  O*  and  the  observation  direction  must  lie  on  the 


same  cone.  In  other  words. 


^ i /\  ' d 

a • t * a • t 


(7.1) 


All  the  quantities  on  both  sides  of  (7.1)  can  be  written  in  terms  of  p 

W *- 

and  as  follows  (Fig.  24): 


o « cos  S ,t  + sin  3 ..z 
w 2 w2 


<3  = sin  3 cos  P x + sin  6 sin  p v + cos  0 z 

o o o o'  o 


t - - sin  <J>  ,x  + cos  p • y 

wz 

cos  3 ,,  * aP  ,/[(a$  ,)‘  + “ ; 

w2  w2  w 2 2 

— TT  < 4>  _ < TT 

— w2  — 

Substituting  (7.2a,b,c,d)  back  into  (7.1),  the  following  transcendental 
Equation  for  4>  0 is  obtained: 


(7.2a) 

(7.2b) 
(7.2c) 
(7. 2d) 


wJ 


*w2 


" (t-)/ a) 1 -t-  p' 

a.  W - 


sin  0 • sin  (^  - ,) 

o o w2 


(7.3) 


Equation  (7.3)  is  solved  ror  3 when  4 and  5 (observation  ancles')  are 

w2  oo 

given.  For  each  fixed  value  of  parameters  4 and  ? . Equation  ('.3)  has  several 

o o 

simple  or  multiple  roots.  Now  let  us  study  some  interesting  limiting  cases. 
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When  0^  is  equal  to  0 or  u (axial  direction),  the  right-hand  side  of 
(7.3)  vanishes.  The  single  root  of  (7.3)  in  this  case  is 

*w2  “ 0 (7-4) 

which  indicates  that  only  the  ray  which  travels  along  the  generatrix  of  the 

cylinder  passing  through  the  source  point  contributes  to  the  field  in  this 

direction.  In  the  case  where  9 « tt/2  and  4 - 0 (direction  normal  to  the 

o o 

surface  of  the  cylinder  at  the  source  point),  again  we  have  the  single  root 
given  in  (7.4). 

It  should  be  noted  that  the  range  of  variation  of  $ , has  bee’  chosen 
to  be  [— ft , tt ] . This  special  choice  excludes  the  rays  which  travel  around 
the  cylinder  more  than  it.  These  rays,  as  mentioned  previously,  have  negligible 
contribution. 

Equation  (7.3)  locates  the  diffraction  points  on  the  upper  wedge. 

A similar  equation  can  be  written  for  the  diffraction  points  on  the  lower 
wedge  (Fig.  25): 


*wl 


- Ki 


sin  9 sin  (i> 
o o 


♦wl> 


(7.5) 


Once  Equations  (7.4)  and  (7.5)  are  solved,  a set  of  values  for  ,,  and 

w «. 

is  obtained,  with  each  of  them  a specific  diffracted  ray  Is  associated. 
Construction  of  these  diffracted  rays  is  the  subject  matter  of  the  next 
subsection. 


' • - Diffracted  Ravs 

Surface  ravs  reaching  P . or  P ,,  are  diffracted  at  these  points  and 

wl  w2 

generate  cones  of  diffracted  rays  with  their  apex  placed  at  these  points.  To 
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analyze  the  propagation  of  the  field  along  the  diffracted  rays,  proper  dif- 
fraction coefficients  and  divergence  factors  for  the  curved  wedge  are  needed. 
7.2.1  Diffraction  bv  edge 

First  let  us  consider  a typical  edge  diffraction  problem  sketched 
in  Fig.  26.  T is  the  edge  of  conducting  surface  £.  According  to  Keller's 
GTD  [98],  [99],  the  total  field  £c  can  be  constructed  as 

EC(?)  - E*(r)  + E^r)  (7.6) 

where  E (r)  is  the  geometrical-optical  field  given  by 

Eg(?)  - 0(-ei)Ei(r)  + 0(-er)Er(r)  . (7.7a) 

-*r 

E and  E are  the  incident  and  reflected  fields,  and 
f 1 x > 0 

9(x)  - ^ . (7.7b) 

[0  x < 0 

i r 

£ and  E.  are  the  shadow  indicators  for  the  incident  and  reflected  fields 
defined  as 


i -*■ 

e (r) 


+1  , r is  in  the  shadow  region  of  incident  field 

-1  , r is  in  the  lit  region  of  incident  field 


(7.7c) 


£ (r)  - 


+1  , r is  in  the  shadow  region  of  reflected  field 

-1  , r is  in  the  lit  region  of  reflected  field 


(7.7d) 


^ 1 / O 

E (r)  is  the  diffracted  field  which  is  of  order  k “ and  is  asymptotically 
smaller  than  E8,  which  is  of  order  kc . As  E8,  Ed  consists  of  two  symmetrical 
parts:  Edi  associated  with  the  incident  field  and  Edr  with  the  reflected  field. 


Ed(rl  - Edi(r)  + Edr(r) 
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where  Sd(r)  (phase  function)  and  (edi  , edr}  are  obtained  by  asymptotic 

m m 


solution  of  Maxwell's  equations  [98].  Referring  to  Fig.  26,  it  is  clear 


that  all  diffracted  rays  emanate  from  the  edge  T of  the  surface  E.  As 


discussed  earlier,  when  an  incident  ray  in  the  direction  a meets  the 


edge  T at  some  point  0 (point  of  diffraction)  and  makes  an  angle  S with 


t (tangent  to  the  edge  at  0) , a cone  of  diffracted  rays  with  semiangle  8 


A 

and  vertex  at  0 and  axis  t is  generated.  According  to  the  law  of  edge 


diffraction"  (counterpart  of  Snell's  law  for  reflection), 


8 = 8 


(7.9) 


Now,  in  order  to  construct  SQ(r)  along  a particular  diffracted  ray  in 


the  direction  a and  passing  through  an  observation  point  r (Fig.  26),  we 


use  the  spherical  coordinates  (o,3,a)  with  pole  at  0 and 


a = distance  from  0 to  r, 


8 = polar  angle  measured  from  t in  the  range  (0,tt), 


a = azimuthal  angle  in  the  plane  perpendicular  to  t. 


Then  it  can  be  shown  [98]  that 


Sd(r)  = Si(0)  + ad 


(7.10) 


where  Sd(0)  is  the  value  of  Sd(r)  at  0. 


The  variation  of  amplitude  functions  e (r)  along  a given  rav,  in  general, 

m 


depends  upon  the  principal  radii  of  curvature  and  R,.  But  in  this  case. 


because  of  the  fact  that  the  edge  T itself  is  a caustic  line  of  the  diffracted 


field,  one  of  the  principal  radii  of  curvature  is  zero. 


(7.11) 


and  the  final  solution  for  e (r)  can  be  written  as 

m 


81 


-*d  -*di 

e (r)  - e (r) 
m m 


+ e (r) 
ra 


0'1/2(0'  + Rx)1/2 


• v2Vi(?,)  * da’ + {1  - r} 


(7.12) 


where  r and  r'  have  coordinates  (a, 6, a)  and  (o'.S.Qt),  respectively,  and 

(i  -►  r}  means  that  the  whole  expression  in  the  right-hand  side  of  the 

equation  is  repeated  for  the  reflected  field.  The  determination  of  the 

initial  values  , in  general,  is  very  involved.  In  GTD  only  are 

m o 

determined  by  comparison  of  (7.12)  with  the  asymptotic  expansion  of  the 
known  solution  to  the  Sommerfeld's  half-plane  problem.  The  zero-th  order 
term  (dominant  term)  in  the  asymptotic  expansion  of  diffracted  field 
can  be  written  as 


E°(r)  - g(ko) 


'l  + ^/Ri  sin  8 


V2ksl(0)rot 

L ° J 


+ {i  - r}  + 0(k“3/2) 


(7.13) 


where  g(kcr  ) is  the  cylindrical  wave  factor, 


g(kod) 


• exp  [-j  (kad  + tt/4)  ] 


2/2nkad 


(7.14) 


1 

n + Jd/RL 


i 

X 


divergence  factor  * DF 


diffraction  coefficient 


(7.15) 


17.16) 


(7.17) 


rot  = rotation  operator 


1 
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DF  (Divergence  Factor)  indicates  the  spreading  of  the  tube  of  rays  emanating 
from  0 as  it  travels  into  space.  In  general,  the  divergence  factor  DF 
depends  upon  the  principal  radii  of  curvature  of  the  wavefront  and  R9, 
but  in  our  case,  as  discussed  earlier,  one  of  these  radii  is  zero.  The 
other  radii  of  curvature  R^  can  be  computed  in  several  ways  [24],  [100], 
[101]  and  [102].  R^  is  the  radius  of  curvature  of  the  diffracted  wavefront 


in  the  direction  defined  by 
'd 


°1  = °2  * 0 


2 sin  g 


1 ~d  ~ 

a * t 


(7.18a) 

(7.18b) 


In  other  words,  R^  is  the  radius  of  curvature  of  a curve  formed  by  the 

* * J 

diffracted  wavefront  intersected  by  a plane  defined  by  t and  o at  the  point 
0.  The  method  described  by  Lewis  and  Boersma  [100]  to  calculate  R^  involves 
the  computation  of  the  Jacobian  of  the  transformation  from  (n,a,a),  where  n. 
is  the  arclength  along  the  edge  measured  positively  in  the  direction  t (Fig. 
26)  to  the  rectangular  coordinate  system.  The  resulting  expression  for  R^  is 


sin"  8 


<Pa  • nr  + (dg/dh)  sin  S 


(7.19) 


where  <p  and  n^  are  the  curvature  and  the  unit  normal  of  f at  0 (iCp  >,  0 and 
n„  points  toward  the  center  of  curvature) . 

Another  procedure  to  determine  R^  has  been  given  by  Deschamps  [102]. 
This  procedure  is  based  upon  matching  the  phase  of  the  incident  and  the 
diffracted  rays  on  T.  The  resulting  formula  for  R^  is 


Rj.  - |-T  +— r-  (51  - 3d) 

I R sin"  S 
o 


1-1 


(7.20) 


* Wave  number  k.1  of  the  incident  surface  ray  has  been  assumed  to  be  the  same 
as  that  of  a plane  wave  in  free  space. 
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where  R is  the  radius  of  curvature  of  the  incident  wavefront  at  0 in 
o 

A A< 

the  plane  through  t and  0 ; in  terms  of  the  principal  radii  of  curvature 

of  the  incident  wavefronts  (R^,R^)  and  their  corresponding  principal 
*i  i 

directions  (x^.x.,)  at  0,  Rq  is  given  by 


2 i 2 i 

1 m cos'  sin  ft 

R1  R1  R* 

o 1 2 


(7.21) 


i ^ i ( a 

where  £2  is  the  angle  between  x^  and  the  projection  of  t on  the  plane  through 

x^  and  x^.  It  can  be  shown  that  (7.19)  and  (7.20)  are  equivalent. 

V1  (Diffraction  Coefficient):  For  a conducting  surface,  the  diffraction 
i r 

coefficients  x ’ are  given  by 


,.i,r 


cosec 


(7.22) 


i r r 

To  define  ip  ’ we  refer  to  Fig.  27,  which  is  actually  the  projection  of  Fig.  26 

on  the  plane  through  0 perpendicular  to  t.  T is  the  tangent  plane  to  I at  0. 

As  illustrated  in  Fig.  27,  the  magnitude  of  namely,  \'l>  ,r|,  is  defined  as 

the  angle  of  rotation  around  t which  brings  a * to  j without  crossing  E, 
i.  r 

and  the  sign  of  ip  ' is 


i,r  i,r,-*\ 

sgn  iji  » e (r) 


(7.23) 


i it 

where  e * (shadow  indicators)  are  defined  (7.7). 

Rotation  Operator  (rot):  The  operator  rot  in  (7.13)  denotes  the  rotation 

y\  i 

about  t that  brings  a onto  a » 


. •'i  -'d 

rot  a m a 


(7.24) 


If  we  decompose  the  amplitude  e^  into  two  transverse  components: 
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-i 

e 

o 


(?) 


ai  i -► 
S e08(r) 


~i  i 

a eoa  r) 


then  one  would  have 


(7.25) 


rot  e*(r)  = 8 e\.(r)  + ctV  (r)  (7.26) 

o op  oa 

^ i '"'i  A ^ 

where  (8  ,ot  ) and  (8  ,ot  ) are  (8, a)  unit  vectors  corresponding  to  incident  and 
diffracted  directions,  respectively  (Fig.  26). 

7.2.2  Diffraction  bv  wedge 

The  previous  results  for  the  edge  diffraction  problem  can  be  applied 

to  the  analysis  of  diffraction  by  the  wedge  after  a simple  modification  of 

i it 

diffraction  coefficient  x ’ • Before  explaining  this  modification,  let  us 
examine  the  geometry  of  the  problem  shown  in  Fig.  28.  The  wedge  is  composed 
of  two  conducting  surfaces,  E^  and  Z^,  intersecting  along  the  curved  edge  E. 

The  half-planes  T^  and  T?,  tangent  at  0 to  E and  Z0,  respectively,  form  a 
planar  wedge  which  approximates  the  curved  wedge  in  the  vicinity  of  0 and 
has  an  exterior  angle  mir  with  1 < m < 2.  When  the  wedge  is  illuminated  by 
the  incident  field  E1  (Fig.  26),  the  total  field  solution  is  again  composed 
of  the  geometrical-optical  field  Eg  and  diffracted  field  E^.  For  the  case  of  a 
wedge,  in  general,  we  have  two  reflected  fields:  E^  from  E^  and  E^  from  E0 
(Fig.  29).  Associated  with  each  of  these  reflected  fields  exists  a proper 
shadow  indicator  (e^  and  c^) , illustrated  in  Fig.  29.  Therefore,  in  the 
case  of  a wedge,  the  geometrical-optical  field  is  given  by 

E®(?)  = 6(-ei)Ei(r)  + 0(-e^)E^(?)  + 0(-£^)E^(r)  . (7.27) 


The  diffracted  field  in  this  case  is  calculated  as  before,  and  the  only 
modification  necessary  is  the  replacement  of  x defined  in  (7.22)  by  the 
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generalized  definition 


2 < " 
— sin  — 
ra  m 


7T  7T  + ii» 

cos  — - cos  z- 

m m 


(for  wedge) 


(7.28) 


where  the  angle  has  a magnitude  equal  to  the  angle  of  the  continuous 

A A (J  A ^ 

rotation  about  t that  brings  o onto  a without  crossing  both  and  T,.  The 
angle  ,,  has  a magnitude  equal  to  the  continuous  rotation  about  t that 
brings  oj  , on  The  signs  of  are  determined  by  e^,r  as  before. 

In  general,  when  both  faces  of  the  wedge  are  illuminated  by  the  incident 
field,  two  values  for  each  of  the  angles  ^ and  are  obtained  which  lead 
to  the  same  value  for 

Using  Equations  (7.6)  and  (7.13)  together  with  the  generalization 
introduced  in  (7.27)  and  (7.28),  the  total  field  solution  according  to  GTD 
for  the  typical  wedge  diffraction  problem  sketched  in  Figure  28  is  obtained. 
Here,  once  again,  we  write  down  the  final  version  of  the  expression  for  the 
diffracted  field  excited  by  the  incident  field  given  in  (7.25)  as 
used  later  in  this  work. 


-*d  **  d 

Ea(r)  - g(kaa) 


'l  + ad/R, 


+ 0(k“3/2) 


(7.29) 


where 


i - r 

X + X 


(’.30) 


d n d ,N  d i i 

g(kc  ),  R^ , 3 , , e^(0),  and  e (0)  have  the  same  definition  as  in  (7.14), 

(7.15),  (7.19),  (7.20),  (7. 251  and  (7.26).  The  coefficients  x*  and  \r  are 


given  by  the  expression  (7.28),  which  is  the  generalized  version  of  (7.22) 


In  the  present  work,  we  apply  the  formula  (7,29)  to  the  special 


case  of  a curved  wedge  formed  by  a circular  cylindrical  surface  intersected 


by  a plane 


7.2.3  Diffraction  by  a cylindrical-planar  wedge 


The  structure  under  consideration  has  been  sketched  in  Fig.  30 


The  incident  field  is  the  creeping  wave  originating  from  Q,  which  meets  the 


The  curve  of  the  wedge  is  a circle  of  radius 


a lying  in  the  xy  plane.  The  direction  of  incidence  at  P _ (diffraction  point) 


o lies  in  the  tangent  plane  to  the  cylindrical  surface  at  P . making  angle 


8 - with  t (tangent  to  the  wedge  at  P 7) . The  polar  unit  vector  corresponding 


to  the  direction  of  incidence,  3 


also  lies  in  the  plane  formed  by  a and  t 


The  azimuthal  unit  vector  a 


The  unit  vector  along  the  diffracted  ray  a is  determined  by  the  direction 
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3p  - Jp  x Od  (7.34) 

Now,  we  consider  different  quantities  in  expression  (7.29)  and  evaluate 
them  for  this  special  case. 

Incident  field  (el,(0),  e1  (0)):  The  incident  field  reaches  the  wedge  at 
od  on 

grazing  angle.  Consequently,  the  only  nonzero  component  of  the  incident  elec- 
tric field  is  oriented  normal  to  the  surface  of  the  cylinder.  In  other  words, 

ep5(P)  - 0 (7.35) 

The  normal  component  of  the  electric  field  (the  phase  factor  eliminated), 

'i  -Iks^P  ) 

namely,  ep^(P^,)e  w2  at  P^,,  is  determined  by  evaluating  the  expres- 

sions (b.5a,b,c)  at  point  P These  formulas  give  the  surface  magnetic 

W 4 

field  at  any  point  on  the  body  of  the  cylinder.  Now,  if  the  surface  field  is 
assumed  to  be  locally  a plane  wave,  the  simple  relationship  between  the 
components  of  electric  and  magnetic  vectors  of  a plane  wave  can  be  utilized 
to  find  the  desired  expression  for  the  normal  component  of  the  electric  field 


4<PvP'"JkS  lP“2> 


H ,(P  ,)  . 

g w2  ' i 

2 *p 


(7.36) 


where  H.,  is  given  by  (6.5)  and  2 is  the  characteristic  impedance  of  free  space. 

P O 

It  is  noted  that  because  of  grazing  incidence,  the  incident  and  reflected  fields 
coincide  so  that  one  half  of  the  total  surface  field  should  be  considered  as 
Che  incident  field  at  the  wedge. 

Computation  of  R^  and  Divergence  Factor:  The  principal  radius  of 

n i ' 

curvature  R,  of  the  diffracted  wavefront  at  P , in  the  plane  (o  and  t)  can  be 
l w2 

determined  by  either  of  the  procedures  explained  in  subsection  (7.1.1). 
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According  co  the  first  procedure,  can  be  found  by  formula  (7.19). 


In  this  formula,  8 is  replaced  by  is  the  curvature  of  the  wedge  which 


is  equal  co  — in  the  present  problem  and 


o - (sin  0 cos  $ , sin  0 sin  $ , cos  0 ) » k 
O O 0 0 o o 


np  - - cos  $u2x  - sin  <J)ij7y  . 


w2 


(7.37) 

(7.38) 


Now  in  order  to  find  the  derivative  of  8 with  respect  to  the  arclength  along 


the  wedge,  d8/dn|  , we  note  that 
'8-6. 


w2 


can  8 « £2/a4>  * ^/n 


(7.39) 


and 


d6 


-e. 


2 , dn 

cos  8 


n2 


(7.40) 


d8 

dr,| 


sin  8 


w2 


(7.41) 


wJ 


w2 


where  °w2  is  the  length  of  the  geodesic  QP^.  Substituting  the  above  expressions 
into  (7.19)  leads  to 


_L_ 

R. 


' “r 


w2  a sin"  g 


(7.42) 


w2 


The  second  procedure  gives  the  same  result  as  in  (7.42); 
however,  it  is  more  complicated  than  the  first  one.  Here,  only  the  main  steps 
of  this  method  are  explained  and  the  details  are  not  discussed.  In  order 
to  apply  this  procedure,  the  first  step  is  to  construct  the  incident  wavefront 
and  to  determine  its  local  geometrical  properties  in  the  proximity  of  the 


. till 
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diffraction  point  P^*  Having  determined  R^  and  R^,  R*  and  R^  can  be  found 
from  (7.21)  and  (7.20). 

The  incident  wavefront  can  be  defined  as  the  geometrical  locus  of  the 
end  points  of  the  rays  emanating  from  the  source  located  at  Q and  having 
equal  lengths,  say,  S^.  The  ray  with  length  Sq,  in  general,  is  composed  of 
cwo  parts.  One  part  is  a geodesic  along  the  surface  of  length  o and  the 
other  part  is  a straight  line  which  can  be  considered  as  the  tangential 
continuation  of  the  geodesic  part  in  the  free  space  with  the  length  d. 

Therefore, 

j + d » S . (7.43) 

o 

When  S is  a constant,  each  ray  defines  a specific  point  S on  the  wavefront 
o 

Si  - const.  (Fig.  31).  Therefore,  one  way  to  parameterize  the  incident 
wavefront  is  to  use  the  (a, 3)  coordinate  system  of  the  cylindrical  surface 
in  the  following  way:  Each  value  of  (a, 3)  defines  a point  on  the  cylinder 
and  also  determines  a ray  joining  Q to  this  point.  This  ray  can  be  continued 
tangentially  in  the  free  space  at  a length  equal  to  d * - j to  reach  the 

point  S which  is,  indeed,  on  the  desired  wavefront.  With  this  parameter- 
ization system,  the  equation  of  the  surface  can  be  written  as 


1 

a 

cos(f 

COS 

2) 

- (So 

- a) 

• cos  3 * 

sin 

,a 

(—  cos 
a 

3)  1 

| 

X(J,3)  - 

a 

sin(— 

a 

cos 

3) 

+ (S 

0 

- a) 

• cos  3 • 

cos 

,J 

(—  COS 

a 

1 

3) 

• 

i 

i 

S 

0 

• sin 

3 

j 

The  principal  curvatures  and  directions  are  now  determined  by  means  of 
expressions  (9.6)  and  (9.7)  in  [43]  or  the  formulas  given  in  (44J.  The 
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coefficients  of  the  first  fundamental  form  I are  listed  below: 


, 3X  3X  <S0  ' ar  4 
E - to  ’ to 2 c°*  6 


■5V  5?  **  0)“  i 

F • H • If ~2 5l"  4 • ';°*  4 


2 .2 

iv  , 0 (S  - 0)  , 2 

C ■ ~ • t4  ■ S*  +• • sin*  S * cos  8 . 

38  38  o 2 


The  coefficients  of  the  second  fundamental  form  II  are 


?■» 

3 X ' o 

e - * N - — 


S - 0 


cos'*  8 


30* 


1 

f .-as. . a - - 

3o38 


o(S^  - a) 


• sin  8 * cos  8 


a 

'y 


- 0) 


- sin*  8 • cos*  8 


where  N is  the  normal  to  the  wavefront. 


3JC  , « 
30  38 


3X  x 2X 
30  ' 38 


(7.45a) 


(7.45b) 


(7.45c) 


(7.46a) 


(7.46b) 


(7.46c) 


(7 . 46d) 


The  "mean"  and  "Gaussian"  curvatures  can  be  derived  through  the  coefficients 
of  the  first  and  the  second  fundamental  forms  via 


E*  - 2(T  .+  «C  . 111  - JL  + 


1 


2 (EG  - F ) 


2S  2(S  - 0) 

o o 


(7 . 47a) 


•»  ~ r - 


'G  EC  - F*  12  ‘ 


(7.47b) 
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Therefore,  Che  principal  curvatures  are 


1 S 


*2 " • 


(7.48) 


In  order  co  find  the  corresponding  principal  direction,  we  note  that 

♦ ♦ 

the  curvature  in  the  direction  u • u^  + u,  -jg  is  given  by 

2 , , 2 

eu  + 2fu  u,  + gu? 

<(t) 1 1 • 

Eu“  + 2Fu^a,  + Gu“ 

Using  (7.49)  together  with  (7.48)  leads  to  the  following  results  for 

A 

principal  directions  u^  and  u.,: 

u ■ (cos  (^  cos  8)  , sin  (^  cos  8)  . 0) 

i A 4 


(7.49) 


(7.50a) 


u,  - (sin  8 sin  (—  cos  8) 


- sin  g cos  cos  g) 

*3 


cos  8)  (7. ‘'Ob) 


For  those  points  of  the  wavefront  that  are  on  the  surface,  we  have 


(surface  is  a caustic) 


(u^.N)  is  a tangent  plane  to  the  cylinder 

u,  is  orthogonal  to  the  cylindrical  surface. 

Based  on  these  relationships,  one  can  conclude  that  the  radius  of  curvature 
of  the  incident  wavefront  in  the  direction  of  the  tangent  to  the  edge  is 


R*  • — * S - a . 
o o 


(7.51) 


In  deriving  (7.51),  Equation  (7.21)  has  been  utilized. 

Substitution  of  (7.51)  into  (7.20)  and  the  exploitation  of  the  fact  that 
' i A 

J • n„  ■ 0 lead  us  to  a result  exactly  the  same  as  what  has  been  given  in 
1 

(7.42). 


Diffraction  Coefficients  x*'r»  The  planar  wedge  whose  faces  are 
tangent  to  the  original  curved  wedge  and  approximates  it  locally  is  a 
right  angel  wedge,  sketched  in  Fig.  32.  Referring  to  (7.28)  and  Fig.  32, 
m is  equal  to  3/2  and 


-2/3  n 

x = 1/2  + cos  (2t/3  + 2ip/3) 


(7.52) 


where  X “ X*  ” and  ^ = ip*  • 'pr . The  magnitude  of  angle  ip  can  be  expressed 
in  terms  of  the  observation  angles  and  the  location  of  diffraction  points 
in  the  following  way: 


COS  (Ip  I 


[gd  - • t)t]  • [a*  - (5*  • t)t] 


sini 2  S 


(7.53) 


w2 


or 


cos  jipj 


sin  9 cos  S » sin  (0  - <t»  0)  +•  sin  3 » • cos  9 - cos  S 

o w2  o w2  w2  o w2 


■ 2 „ 
sm  3 


w2 


(7.54) 


Using  the  relation  cos  gw?  = a • t = sin  So  • sin  (<p  - -p^)  * 

Equation  (7.54)  may  be  simplified  as 


cos  'P ! 


cos 

sin 


(7.55) 


i r 

The  sign  of  tp  is  determined  by  e ’ . Therefore,  in  order  to  define  ip 
completely,  the  position  of  the  observation  point  with  respect  to  the 
shadow  boundary  (S.B.)  should  be  examined.  Fig.  30  shows  that  the  plane 
(t,c  ),  which  is  tangent  to  the  surface  of  the  cylinder  at  is  the  local 

shadow  boundary  for  the  incident  surface  ray.  Therefore,  the  observation 


directions  (9  ,i  ) which  satisfv 
o o 
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sin  0o  cos  «»o  - 4>w2)  < 0 (7.56a) 

or 

l*o  ' *w2l  > 11/2  (7.56b) 

belong  to  the  shadow  region. 

The  singular  nature  of  x causes  some  difficulty  when  the  observation 
point  approaches  the  shadow  boundary.  As  a matter  of  fact,  when  •■h  tends 
to  zero,  x becomes  infinitely  large.  This  very  well-known  deficiency  of 
Keller's  GTD  diffraction  coefficient  has  been  the  subject  of  investigation 
by  many  authors.  Several  generalizations  and  modifications  have  been 
introduced  to  make  GTD  uniformly  valid  over  the  transition  region  (around 
the  shadow  boundary).  Among  the  "Uniform"  theories,  we  can  refer  to  UTD 
[101],  [25]  and  UAT  [100],  [98].  Another  approach  to  circumvent  this 
difficulty  is  to  employ  STD  (Spectral  Theory  of  Diffraction),  developed  by 
Mittra  et  al.  [82],  [83].  This  approach  attempts  to  correct  the  root  of 
this  difficulty  by  closely  examining  the  earlier  stages  of  the  derivation 
of  GTD  formulas  rather  than  attempting  to  modify  the  final  results  post-fact. 

In  the  familiar  example  of  diffraction  by  the  half-plane  considered  by  them, 
it  was  shown  that  the  singular  nature  of  the  diffracted  field  comes  from  the 
fact  that  the  spectrum  of  the  induced  current  on  the  half-plane  contains  a 
pole  in  the  direction  of  incidence.  When  the  observation  point  is  far  from 
the  shadow  boundary,  the  simple  saddle-point  technique  applied  to  the  integral 
of  the  spectrum  of  the  induced  current  leads  one  to  Keller's  GTD.  However,  when 
the  shadow  boundary  is  approached  in  the  observation  space,  the  pole  becomes 
very  close  to  the  saddle  point  of  the  integrand,  the  simple  method  used  in 
the  previous  case  fails,  and  more  refined  methods  are  required  to  carry  out  the 
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current  integration  (e.g.,  introduction  of  Fresnel  integrals,  ...).  It 

was  also  pointed  out  that  the  existence  of  the  pole  in  the  spectrum  of  the 

current  can  be  considered  as  the  contribution  of  the  physical  optics  currents 

of  semi- inf inite  extent.  Therefore,  the  scattered  field  was  interpreted  as 

the  superposition  of  the  spectrum  of  a physical  optics  current  and  a 
-1/2  1 

higher-order  (k  ) component,  called  J , without  the  pole  singularity. 

In  the  later  application  of  STD  to  the  strip,  the  diffracted  field  was 
considered  proportional  to  the  spectrum  of  the  induced  current  J,  composed 
of  three  terms, 

J * J + J.  + J.,  (7.5/) 

where  J*5  is  the  physical  optics  current  on  the  strip  and  J^r,  are  the  0(k  ^“) 
contribution  of  the  edges  1 and  2 truncated  over  the  strip. 


J 


tr 

1.2 


(7.58) 


1 -1/2 

i here  J are  the  0(k  ) currents  Induced  on  the  semi- inf inite  half-planes 

erected  at  edges  1 and  2,  and  J.  , are  portions  of  these  currents  outside  the 

strip.  Final  results  for  the  field  were  the  same  as  those  predicted  by  UAT. 

Generalization  of  this  procedure  to  the  curved  surfaces  required  some  special 

remarks  on  the  definition  of  J*3.  When  the  case  of  the  curved  wedge  is  dealt 

with,  J*3  cannot  be  simply  replaced  by  the  physical  optics  current  on  the 

-1/2 

fictitious  half-planes  tangent  to  the  wedge  faces.  The  0(k  ) current 

t r 

components  J 0,  more  or  less  represent  localized  distribution  of  current  in 
proximity  of  the  wedge,  therefore,  the  local  properties  oi  the  surface  should 
be  sufficient  to  provide  us  with  a fairly  accurate  estimate  of  that.  Whereas 
jk  corresponds  to  the  global  distribution  of  the  Induced  current  on  the  whole 


the 
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body  which  can  deviate  substantially  from  what  one  would  obtain  based  upon 
the  local  geometry  of  the  surface  around  the  wedge.  Initial  results  indicate 
that  using  the  approximate  expression  for  the  surface  current  (instead 
of  physical  optics  current  on  the  tangent  half-plane)  for  J*3  will  hopefully 
improve  the  accuracy  of  the  existing  solution  for  a curved-surface  diffraction 


problem. 

In  applying  the  above  principles  to  our  problem,  J*3  has  been  obtained 

by  truncating  the  current  distribution  that  would  exist  on  an  infinite 

cylinder.  Thus,  the  radiated  field  due  to  is  derivable  by  subtracting 

the  contribution  of  the  currents  of  the  extension  of  the  finite  cylinder 

from  the  field  radiated  when  the  cylinder  is  infinite  (Chapter  8).  The 
- t r 

contributions  of  J ^ to  the  far  zone  scattered  field  are  proportional  to 
the  spectrum  (Fourier  Transform)  of  these  currents  which  are  computed  under 
the  same  assumptions  employed  in  GTD.  These  contributions  do  not  have  any 

W 1 Wfc 

singularity  and  behave  smoothly  over  the  whole  space.  The  fields  F and  F 

(diffracted  from  the  lower  and  the  upper  wedge,  Figure  19)  are  generated  by 
t r t r 

and  J,  and,  based  upon  these  principles,  can  be  approximately  deduced 

from  the  GTD  solution  upon  the  subtraction  of  the  pole  singularity  from 

Keller's  diffraction  coefficients  \.  In  conclusion,  the  point  is  emphasized 

once  more  that  using  GTD's  \ as  it  stands  provides  us  with  a scattered  field 

due  to  the  total  induced  current  on  an  infinite  half-plane,  whereas  in  our 

w 

formulation  (see  Fig.  19),  F is  just  the  localized  effect  of  the  wedges  on 
the  total  diffracted  field.  As  a matter  of  fact,  the  field  due  to  the  current 
induced  on  the  body  has  been  calculated  separately  and  represented  by  Fs. 

In  view  of  the  above,  the  following  definition  has  been  used  for  \ 
for  the  computation  of  FW 
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Xp’r  * x1,r(GTD)  - xP°  (physical  optics)  (7.59) 

i r 

where  x ’ (GTD)  has  the  expression  given  in  (7.52)  and 

A 

X (physical  optics)  = cotan  -%■  * — jjp—  + •••  (7.60) 

and  we  also  have 

A-1 

X(GTD)  - — -=-  + Aq  + + •••  . (7.61) 


In  summary,  the  final  expression  for  the  far-zone  diffracted  field  due 
to  diffraction  at  wedge  point  in  the  observation  direction  (90»<P0)  can 
be  obtained  by  substituting  (7.59),  (7.42),  (7.36)  and  (7.35)  back  into  (7.19) 
which  leads  to 


(r) 


-jkr  -jir/4  ZkR, 
e e i 

kr  ’ Sin  Sw2 


'^oV!w£l 

..  2 J 


• exp{jk[a  sin  0 • cos  (0  - tp  _)  - a sin  0 cos  $ + cos  0 ]}  (7.62) 

O O W J.  o o i.  o 

-*d 

Several  remarks  should  be  made  on  (7.62).  First  of  all,  E^.,  is  polarized  along 
olp,  because  the  tangential  component  of  the  incident  field  at  is  zero. 

The  observation  vector  r is  connecting  Q (the  source  point),  which  is  chosen  to 
be  the  phase  reference  center,  to  the  observation  point.  The  diffraction  factor 
Dp  is  given  by 


D 


h 

F 


(7.63) 


It  should  be  noted  that  ?^_(r)  is  the  diffracted  field  due  to  a point  on 

w2 

the  upper  wedge.  Diffracted  fields  generated  by  the  points  on  the  lower  wedge 
are  constructed  in  the  same  way  only  when  is  replaced  by  and  geometric 
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quantities  of  the  lower  wedge  points  replace  those  of  the  upper  ones. 
Formula  (7.62)  gives  the  t'ar-zone  scattered  field  in  terms  of  its  single  dp 
component.  The  unit  vector  dp  depends  upon  the  location  of  the  diffraction 
point  P which  can  also  be  expressed  in  terms  of  its  components  along  the 

transversal  unit  vectors  along  the  observation  direction; 


cos  (4>  - $ -)  . 

A Q O V 4.  A 

ap  " sin  I6w2i  0 


cos  3 


sin  (0 


4w2> 


sin  3 


(7.64) 


w2 


7.2.4  Blockage  of  diffracted  rays 

Some  of  the  wedge  diffracted  rays  may  not  be  able  to  reach  the  observation 
point  because  of  their  interception  by  the  body  of  the  cylinder.  The  condition 
of  the  blockage  by  the  object  can  be  determined  easily  by  looking  at  the 
geometry  of  the  structure.  Blockage  of  the  upper  wedge  diffracted  rays  occurs 

when  0 > •=  and  the  observation  point  is  in  the  shadow  region.  Lower  wedge 

o 2 

diffracted  ravs  are  blocked  when  6 <4  and  the  observation  point  is  again  in 

O 4 

the  shadow  region. 


7.2.3  Caustic  directions 

Equations  (7.19)  and  (7.62)  express  the  dependence  of  the  amplitude 
of  the  diffracted  field  upon  R^  as  it  propagates  along  the  diffracted  ray. 

R^  is  the  distance  between  one  of  the  caustic  lines  and  the  diffraction  point. 
This  parameter  determines  the  rate  of  divergence  of  the  diffracted  rays.  As 
it  increases,  the  pencil  of  the  diffracted  rays  diverges  more  slowly,  and  con- 
sequently, the  diffracted  field  in  that  direction  increases.  In  the  limiting 
case,  in  those  directions  where  Rj  becomes  infinitely  large,  the  diffracted 
field  according  to  GTD  blows  up.  However,  it  is  known  that  real  fields  remain 
finite  in  these  regions.  These  directions  are  called  "caustic  directions"  and 
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are  determined  by  Equation  (7.42).  From  the  above  explanation,  it  is 
obvious  that  GTD  cannot  be  applied  in  these  directions  unless  it  is  modified 
in  some  proper  manner. 

In  order  to  derive  the  desired  modification,  we  reexamine  the  integral 
representation  of  the  far  field.  The  expression  for  the  far-zone  radiated 
field  can  be  shown  to  have  the  following  general  form 

r r jkfasini?  cos (4>-<p  )+zcos0  ) 

F(0  ,$  ) - I ip(0  .4>  ) * f(>J>.z)  * e ° ° ds  (7.65) 

o o J J oo 

S 

(derivation  is  given  in  the  next  section),  where  S is  the  surface  of  the  cylinder 
and  f is  a typical  component  of  the  surface  current.  For  large  values  of  k, 
the  stationary  phase  method  for  multiple  integrals  can  be  applied  to  derive  an 
asymptotic  expansion  for  (7.65).  Because  of  the  finiteness  of  S,  in  general, 
we  have  two  classes  of  critical  points  for  the  phase  of  (7.65).  Class  1 consists 
of  those  stationary  points  which  are  in  the  interior  region  of  S,  and  their 
contribution  to  the  total  field  can  be  identified  with  the  surface  rays.  The 
other  class  contains  the  critical  points  which  locate  on  the  boundary  of  S 
or  the  wedge  of  the  cylinder.  The  latter  class  generates  the  wedge  diffracted 
rays.  Based  on  this  argument,  it  is  expected  that  the  wedge  diffracted  field 
can  be  represented  by  a one- dimensional  Integral  around  the  wedge  stationary 
phase  points.  Some  of  the  authors  ([95],  [105],  [106],  [107])  have  attempted 
to  identify  this  integral  with  the  integral  representation  of  a fringe  current 
flowing  along  the  rim  of  the  wedge.  Ryan  et  al.  [105]  applied  this  technique 
to  evaluate  the  field  at  the  axial  caustic  formed  by  an  axially  incident 


plane  wave  upon  the  edge  of  planar  apertures,  disks,  or  the  wedge  formed  bv 
a ring  slope  discontinuity . In  this  procedure,  the  GTD  diffraction  coefficients 
for  the  planar  wedges  and  edges  are  used  to  determine  an  equivalent  current  at 
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the  edge  of  a discontinuity.  These  equivalent  currents  are  then  integrated 
using  the  radiation  integral  to  determine  the  scattered  field.  As  a matter 
of  fact,  the  diffraction  coefficients  specify  the  diffraction  at  each  point 
on  the  wedge  and  the  radiation  integral  sums  these  contributions. 

The  major  drawback  of  this  method  is  the  way  one  defines  the  "equivalent 

edge  current".  The  function  determined  by  the  diffraction  coefficients  is  not 

really  a current  in  a physical  sense.  First  of  all,  the  expression  used  as 

"equivalent  edge  current"  is  not  only  a function  of  the  wedge  point  but  also 

depends  upon  the  observation  angles.  Second,  at  points  on  the  wedge, 

removed  from  >£  - 0,  the  diffracted  ravs  lie  on  a cone  and  therefore  do 

w 

not  contribute  to  the  radiated  field  in  all  directions  in  the  entire  space. 
Consequently,  we  should  not  expect  accurate  results  by  this  technique  in  the 
regions  far  from  the  caustic,  3nd  the  "equivalent  current"  representation  may 
be  valid  only  in  a close  neighborhood  of  the  stationary  point  which  contributes 
significantly  to  the  field  in  the  caustic  region. 

In  this  report  we  have  applied  a procedure  very  similar  to  that 
discussed  in  [105].  The  type  and  expression  of  the  equivalent  edge  current, 
which  should  be  placed  on  the  wedge  of  the  cylinder,  are  determined  bv  compar- 
ison between  the  GTD  formulas  and  the  asymptotic  expansion  of  the  integral 
representation  of  the  far  field  radiated  by  the  fringe  current.  Because  of 
the  normal  polarization  of  the  incident  electric  field  at  the  wedge  (Fig.  30) 
and  i^-polarization  of  the  diffracted  field,  the  equivalent  edge  current  is 
assumed  to  be  of  magnetic  type  K(f>).  The  radiated  field  can  be  expressed  in 


terms  of  this  function  as 
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, jli'a  cos  6 -jkr  -jk(asint)  cos<i>  -i’,cosd  ) 
„w2  J o . e o o 2 o 


4* 


kr 


• e 


r 

. 


- jk<J(4>)+jkaslnc1  eost,.f>-j'  ' 

K (^)e  ° ° sin  (i>  - J>)d4* 

o o 


( 7. 66a) 


, 2 -Jkr  -Jk(asin0  cos4>  -<Lcosd  ) 

„w2  - Ik  a e J o o2  o 

K ■ ■ -**  — 1 ■ • ■ ■ ■ ■ • p 

6 4 it  kr 


$ 

f -jkcr(^)+Jkaslneocos(^-to) 

K (^)e  1 cos(l>  - i^di1 

. o ° 

4>w2~A<> 


(7.66b) 
( 7 . bob ) 


where  K(^)  has  been  assumed  to  have  the  following  form: 


K($)  ■ K ('j)e 
o 


■jka(^) 


( 7 . 66c) 


2 '»!/'» 

j(f)  - [(ac)"”  + t“]  ■ geodetieal  distance  between  the  source  Q and 

an  arbitrary  point  on  the  wedge  defined  bv  4>. 

( 7 . obd) 

Equations  (7.ooa)  - (7.6od)  have  been  written  for  the  upper  wedge.  Similar 
expressions  can  be  derived  for  the  lower  wedge.  Here,  once  more,  we 
emphasize  that  we  would  prefer  to  consider  (7.o6)  as  an  approximate  integral 
representation  of  the  field  in  the  caustic  region  rather  than  interpreting  K(j) 
as  an  equivalent  edge  current  which  generates  the  diffracted  field.  The  angle 


> determines  the  location  of  the  stationary  phase  point  on  the  wedge 


corresponding  to  the  caustic  direction.  It  is  believed  that  a one-d imensional 
integral  representation  for  the  field  similar  to  ('.no)  can  be  derived 
by  treating  asymptotically  the  exact  surface  integral  representation  of  the 
far  field  near  the  wedge. 


. A .W  . 


J 
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As  discussed  earlier,  our  basic  criterion  for  determining  the  unknown 


function  K ($)  is  the  assumption  that  one  must  recover  the  GTD  results  bv 
o 


asymptotically  expanding  (7.66)  around  the  stationary  phase  points  far  from 


those  corresponding  to  caustic  directions.  To  this  end,  we  apply  the 


stationary  phase  method  to  (7.66).  Integrals  given  in  (7.66)  are  of  the 


following  general  form: 


IU  ,0  ) - h(4)ejk^U)d<f 

o o 


(7.67) 


where 


Q(4>)  - a sin  6^  cos  ($  - ,$o)  - + i~  . 


(7.6S) 


The  stationary  phase  point  is  given  by 


&•(?) 


- a sin  0 sin  ($-$)- 
o o 


(7.69) 


fit  + (a*) 


sin  0 sin  (*t>  - $) 

o o 


(7.70) 


ft.  + (ad>)' 


The  stationary  point  4>  is  the  root  of  (7.70),  which  is  the  same  as  (7.3). 


Therefore, 


0 ■ $ - 
s w2 


(7.71) 


In  other  words,  the  stationary  point  $ of  fl(*>)  coincides  with  the 

diffraction  point  corresponding  to  the  observation  direction  (0  ,*•  ) given 

o o 

bv  a , . At  j * j . we  have 
w-  * 


9 •> 
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where  0%  - o(^)  and  ^ - sin"1  (t/oj.  Using  (7.42),  (7.  72)  can  be  rewritten 
as 


2 2 

-a*  sin'  3 

s 


(7.73) 


Formula  (7.76)  is  used  to  derive  the  first-order  term  in  the  asvmptoti 
expansion  of  (7.o6K  The  final  results  are 


1Q9 


, -jkr  +jir/4  /k|R  | jk[asin0  cos(4>  -<$  )-asin0  cos<t>  +<Lcos0  J 
,w2  e . e . / 1 . o so  o o 2 o 


kr 


2/2n 


sin  8 


cos  (d  - <fi  ) 
o s 


sin  6 


• [-jk  • sin  S • K(4>  ) ] . 


(7.77b) 


Comparison  between  (7.77)  and  GTD  formulas  (7.62)  and  (7.64)  determines  K(4>): 


K(f)  - 


ZQ  • D“(4>)  Hg($) 


jk  sin  S(4>) 


(7.78) 


,(2). 


In  a caustic  region  where  R,  + ® or  Q (4>  ) ■*  0,  formula  (7.76)  is  not 

i s 

valid  any  longer,  and  the  third-order  derivative  of  the  phase  should  be  taken 
into  account.  The  expansion  (7.74),  in  this  case,  is  replaced  by 


q(2> / * \ 

“ ($_)  •>  ••  \v  / -i 

fl($)  = r.(4  ) + ^ — <*-*)  + — (i>  - $ ) 

s 2 s b s 


, r.(3)(*_) 


(7.79a) 


where 


(3)  f3a  sin“  3s 

' ($  ) * a cos  6 
s s 


- 1 


(7.79b) 


and  the  integral  in  (7.75)  takes  the  following  form, 

,(3) 


r?.i2)  . 

a , 2 

s jk! (4>-0  ) + 


<®  ) V 

-T-^-(^c) 
t)  s 


(7.80) 


,(2) 


At  the  caustic  direction  where  (0  ) vanishes,  a fairlv  accurate  estimate 

s 

of  (7.30)  can  be  obtained  by  first  extending  the  domain  of  integration  to 
infinity  and  then  expressing  the  resulting  infinite  integral  in  terms  of 
Gamma  functions. 


> 


(1/3) 


In  very  special  regions  (tor  example  at  6 ■ tt/2  and  $ c it)  , ) 

oo  s 

and  become  extremely  small  simultaneously  and,  hence,  higher-order 

caustics  will  occur.  These  situations  can  be  treated  by  including  higher- 
order  terms  in  the  expansion  of  the  phase  function 


where 


and  the  phase  function  in  the  integral  (7.75)  should  be  replaced  by  its  value 


given  in  (7.32b).  Again,  in  this  case,  in  those  directions  where  i! 


and  H ) vanish  simultaneously,  the  integral 


can  be  approximated  in  the  same  manner  discussed  in  the  previous  case.  In 


this  case,  the  final  result  is 


Ill 


A careful  study  of  the  expressions  and  reveals  that 

their  simultaneous  vanishing  is  not  possible  and,  hence,  higher-order  caustics 

cannot  occur.  Also,  it  is  obvious  that  closed-form  expressions  like  (7.76), 

I 

(7.81),  and  (7.84)  cannot  be  used  in  the  proximity  the  caustic  directions. 

One  possible  approach  is  to  use  the  integral  (7.83)  together  with  the 
expression  (7.82a)  and  derive  a transition  function  to  replace  the  factor 

e+jir/4 


in  the  formula  (7.62)  for  the  field  in  the  caustic  region.  This  replacement 
provides  us  with  a fairly  smooth  transition  into  the  caustic  region  from  the 
observation  directions  where  the  formula  (7.62)  is  applicable. 


a sin 


(7.85) 
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8.  TRUNCATION  EFFECT 


According  to  our  formulation  of  the  present  problem  given  in  Chapter 

A,  the  far-zone  raaiated  field  was  decomposed  into  the  upper-  and  lower- 

wl  w2 

wedge  contributions  F ’ and  that  generated  by  the  current  induced  on 

s 

the  surface  of  the  cylinder  F . In  Chapter  5 the  far  field  due  to  a point 
source  located  on  an  infinitely  long  cylinder  was  calculated  and  the  surface 
ray  representation  was  used  to  express  the  final  results.  Frequently,  the 
far-zone  field  formulas  for  the  infinite  cylinder  are  used  as  an  approximation 
for  the  case  of  the  finite  cylinder.  However,  this  technique  yields  inaccurate 
results,  especially  for  the  observation  points  near  the  axis  of  the  cylinder 
where  the  expressions  given  in  (5.1)  and  (5.2)  are  not  valid,  or  in  the 
cases  where  the  distance  of  the  source  from  one  of  the  ends  of  the  cylinder 
or  the  total  length  of  the  object  is  not  large  enough.  In  order  to  carry 
out  a more  accurate  approximate  solution  to  this  problem,  the  current 
distribution  on  the  finite  cylinder  will  be  taken  identical  to  the  distribution 
on  the  infinite  cylinder  [73].  Then  the  contributions  of  the  portions  of 
the  infinite  cylinder  above  z = l ^ and  below  z = (Fig.  19)  are  determined 

with  the  assumption  that  the  current  on  these  semi- inf inite  cylinders  remains 
unchanged  when  the  cylinder  is  truncated.  Subtraction  of  the  contributions 
of  the  two  semi- infinite  cylinders  from  the  infinite  cylinder  field  yields 
a better  estimate  for  the  far-zone  radiated  field  due  to  the  finite  cylinder. 

In  this  chapter,  we  derive  the  approximate  expressions  for  the  far-zone 
radiated  field  due  to  current  distributions  of  the  upper  and  lower  semi- 
infinite  cylinders. 


8.1  Far  Field  Radiated  bv  Upper  Semi-Infinite  Cylinder 


The  far  radiated  field  due  to  a current  distribution  confined  to  a 
region  of  space  with  a finite  extent  can  be  expressed  as 
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A(r)  |acos0Q  • TUp  [sin(*Q-i)>) J^]  - asinOg  TUP  [J^] 
A(r)  ja-TUp  [cos (<J>0*-$) 


(8.6a) 

(8.6b) 


where 


A(r) 


ziii- 

4ffY0 


e~jkr  -jkasinSQ- cos«)!q 


(8.6c) 


and  T^p  is  the  two-dimensional  Fourier- type  transformation  defined  as 


T P [h(*,*)] 


2tt  00  j [ka  sinQ  cos(<H>n)+kz  cos©  ] 

/ d<f)  / dz  h (4> , z ) e 
0 i7 

(8. 6d) 


It  is  obvious  that  using  the  exact  solutions  for  J and  J in  the  above 

z 

equations  makes  the  double  integrations  a formidable  task  to  perform, 
especially  in  the  high-frequency  range.  Besides  this  point,  we  should 
bear  in  mind  that  expressions  being  derived  in  this  chapter  are  only 
corrective  terms  to  our  final  solution.  Based  on  these  facts,  approximate 
expressions  (6.8)  were  chosen  to  represent  the  surface  magnetic  field  or 
surface  current  components  on  the  semi-infinite  cylinder.  The  expressions 
(6.8)  have  been  shown  to  be  fairly  accurate  estimates  of  the  surface  field 
when  the  observation  point  is  two  or  more  wavelengths  far  from  the  source. 
Therefore,  we  have  (Fig.  33) 

J (P)  * H (P)  - (M-?)  [H,  sin2S  + H cos23]  (8.7a) 

Z $ BO 

J (P)  --H  (P)  » (M-z)  [-H  cos2S  - H sin23] 

i z So 


(8. 7b) 
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where  H (P)  and  H,(P)  are  given  bv  (6.8a,b).  The  special  form  of  the  final 

o p 


expressions  for  J and  J are  such  that  we  can  represent  them  as 

Z 'P 


J.(P)  - JzO(P)  ' e 


-jko 


y>  - Vp)  ' * 


-jka 


(8.8a) 

(8.8b) 


where  j is  the  arclength  of  geodesic  QP  (Fig.  33).  Substitution  of  (8.3) 
into  (8.6)  leads  us  to  the  double  integrations  of  the  following  general  form: 


F (VV 


2/  * jk[asin0  cos(^-*  )+z  cose.] 

/ ad*  / dz  f (*,z)  e jkJe  0 0 0 

0 l 

(3.9) 


, 1/2 


2 2 U 

where  a - [z“+(a*)  ] . Here,  we  attempt  to  simplify  the  computation  of  F 

by  deriving  an  asymptotic  expansion  for  the  inner  integral  (with  respect 


to  z)  on  the  assumption  that  k is  large.  To  this  end.  we  rewrite  (8.9) 


as 
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F <W  " ^ ad<,> 


jka  sine  cos(*  -*) 


0 


(8.10) 


where 


■sU 


f (♦)  - / dz  f(*,z)  e 

l. 


, 2 1 

Jk[z  cos0q  - •'(a*)  +z“  ] 


(3.11) 


and  then  apply  the  stationary  phase  method  to  (8.11).  At  this  point  one 
might  raise  the  question  why  it  is  not  possible  to  apply  this  method  to  the 
original  double  integral  (3.9).  Theoretically , it  is  possible  to  evaluate 
(8.9)  by  the  two-dimensional  version  of  the  stationary  phase  method,  however. 


U-..L—  
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it  can  be  shown  [85]  that  the  stationary  point  of  the  phase  of  the  integrand 
lS.9)  is  of  second  order.  This  fact  makes  the  resulting  two-dimensional 
asvmptotic  expansion  very  slowly  convergent  and,  in  addition,  the  coefficients 
of  the  double  infinite  sum  become  highlv  complicated  when  order  of  the  terms 
increases. 

In  order  to  apply  the  stationary  phase  method  to  (8.11),  first  we 
rewrite  it  as 

fU(>)  - / f(f,z)  eJkqU’*)  (8.12a) 

<2 


where  q(z,$)  ■ z cos30  - /(a$)“+z“  (8.12b) 

The  stationary  point  of  q is  given  by  the  following  equation 

3*  * COS90  - ■ 0 <8-13> 

► vai) ‘+z~ 

from  which  the  following  expression  for  the  stationary  value  of  z as  a 
function  of  $ is  obtained; 

( 4> ) ■ a|$:  cotgt)0 


(8.14) 
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The  phase  function  q(z,$)  can  be  expanded  around  z 


q(2)(V*>  2 

q(z,$)  : q(zg,4>)  + 2 ( *~za ) + ... 


where 


(8.15a) 


q(zg,<fr)  - -a1o|stn0(. 


(8.15b) 


(2) , . -Sln3e0 
q (Zs’|)  a i ! 


< 0 (^0) 


(8.15c) 


The  case  in  which  * 0 is  of  no  interest  to  us  because  z (i)  -*■  0,  and 

s 

if  kt.,  >>  1,  then  the  stationary  point  is  located  far  outside  the  domain 

of  integration  and  thus  has  negligible  contributions  to  the  integral. 

The  most  significant  contributions  to  a typical  integral  of  the  form 

given  in  (8.12a),  when  k is  large,  come  from  the  end  point  z ■ C,  and  the 

stationary  point  z » zg(i).  When  these  two  points  are  far  apart,  these 

contributions  can  be  computed  separately  and  then  added  to  give  an  estimate 

for  (3.12a).  Of  course,  if  the  stationarv  point  z is  far  from  the  end 

s 

point  and  outside  the  domain  of  integration,  only  the  end  point  contributes 
to  the  integral.  A more  complicated  situation  can  occur  when  the  two  points 
approach  each  other.  In  this  case,  the  simple  formulas  for  calculating  the 
contributions  of  a single  isolated  end  point  or  stationary  phase  point  are 
not  applicable  any  longer.  The  details  of  the  procedure  for  treating  this 
case  have  been  explained  by  many  authors  (for  example,  see  page  421  of  [108]>. 
We  apply  the  final  formulas  to  (8. 12a) . The  first-order  term  in  the  asymptotic 
expansion  of  <3. 12a),  according  to  the  above-mentioned  procedure,  is 
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which  Is  the  well-known  formula  for  computing  tin*  contribution  of  a single 
isolated  stationery  phase  point.  The  factor  1.'-  In  front  of  this  formula 
accounts  for  the  semi-infinite  integration  Interval.  When  z is  ^ulte  far 
from  or  »k  s,(f'  " 1,  using  the  following  asymptotic  expansion  for  Q, 


+ »ft»  i-Rec) 


IS. 18a> 


x ' 0 


x < 0 


the  expression  (3.12a)  reduces  to 


7U  . 

f 1.')  . 


%—  Jkq(*  ,*)+J-v4 

-rr^ fis„')e  . S(«  -( ,) 

■ 


f ( 0 , • & ) • e 
Jkq’ («,,*) 


JMu  .,*>) 


(-)  , v 

, q (*  .*)  . 0 


The  important  limiting  cases  (8.17)  and  (3.18)  indicate  that  the  formula 

(3.16a’',  indeed,  satisfies  the  criteria  discussed  earlier.  In  the  special 
case  where  qu.i)  Is  given  bv  <v8.l2b),  one  has 
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s,U)  - sgr.U  -z  ($)) 
- 2 s 


i ♦ I Sin60-iocose0+  /(a$;  +* 


2.  .2 

9 


1/2 


(8.19a) 


<♦>  - /-iilii  . e-J"/4 

S * sin  0, 


(8.19b) 


h,(*) 


1^(0  I 


• e~J"/4 


l. 


ICOS0Q  - 


/ 2 2 
/(a*)N< 


(8.19c) 


Let  us  denote 


*2U)  “ S2U)  'k/2 


(S. 20) 


then 


Q(s2(*)»lt  e-1"^4)  - Q [ ( 1+J  ) ;,,($)] 


(3.21) 


Equation  (8.21)  can  also  be  written  as 


r-  (1+j);,  2 

Qfd+J)  «,(*)]  - J~  - i * e"X  dx 


(8.22) 


After  proper  changing  of  the  variable  and  some  algebraic  manipulations,  one 
can  express  (S.22)  in  terms  of  the  Fresnel  integrals  C and  S, 


Q[(i+j):2)  ■ , 2 


*S  - 


j ( 1-(1+J)[C(-^)  - jS(— =■)] 


C3.23) 


where 


» I' 


S‘W"al  ““  1 taM.  for  various 

Z8"  " "SU"""  °£  th"'  functions 

efficient  v.y  to  conput.  C a„0  s.  vPich  peen  di9cus„d  ^ ^ 

“ th'"  * - - — functions  f#  a„d  ^ 

reasonably  simple  expressions  : 


where  f (2) 


are  given  by  following  approximate 
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Substitution  of  (8.25)  in  (8.23)  yields 


2c,  -J2c,  2c,  -J2c! 


Q[(l+j)C,l  «/l  + JfQ  (~~)  e “'2  -g0P)  • e ~'2]} 


C,  > 0 


(8.26) 


When  c,  is  negative,  it  can  be  shown  that 


Q[-(l+j)  IcJ  J - ^ - Q[  (1+j  ) |c,  | ] 


(8.27) 


Combining  (3.26)  and  (3.27)  results  in 


Q[(1+1K2J  “ ^ ^ 9 (— ^ “ s8n(c0) 


2|c,!  -j2c;  2|c  I -J2c2 

[J  f0  (-^>  « 2 - So(~  e 2] 

C »tt 

» IT 


(8.28) 


q[(i+j):21 


**  e(-;2)  - w(c0)e 


(8.29a) 
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where 


The  above  expression  for  Q is  employed  in  (8.16a)  which  together  with  (8.19) 


leads  us  to  the  following  expression  for  fu  (4) 


The  function  f (o)  is  used  in  (8.10).  Then  the  j-integracion  is  carried  out 


When  the  arbitrary  function  f(o,z)  in  (8.9)  is 


numerically  to  obtain  F 


replaced  by  J n sin  ( 0> r.— 4> ) , J n cos^-o),  and  J (given  in  (8.3)  and  (8.7)) 


various  terms  of  the  expressions  (3.6)  for  the  components  of  the  radiated 


field  due  to  the  current  distribution  on  the  upper  semi-infinite  cylinder 


Far  Field  Radiated  bv  Lower  Semi- Inf inite  Cylinder 


semi- infinite  cylinder.  The  final  result  is  exactly  the  same  as  for  the 


previous  case  except  for  minor  changes  in  some  parameters.  As  a matter  of 


z- integration  is  carried  out  over  the  interval  (~» 
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The  total  contribution  due  to  those  parts  of  the  cylinder  that  have 
been  removed  is  the  sum  of  the  fields  generated  by  the  upper-  and  lower- 
semi-infinite  cylinders  as  discussed  above.  It  is  expected  that  subtraction 
of  the  radiated  field  generated  by  the  truncated  parts  of  the  structure  from 
the  total  field  will  improve  the  overall  accuracy  of  the  results. 


^ - - il.i 


APPROXIMATE  EXPRESSION  KOR  THE  El EU)  IN  PARAX1A1  REGION 


q. 

Squat  Ions  (5.11  and  (.  used  l»  Chapter  2 tor  computing  the 
contribution  of  surface  rava  arc  not  valid  In  tho  paraxial  region  0 . 

or  b ^ : t (this  point  la  discussed  in  more  detail  In  IRS)  or  in  Sec.  J..'), 

It  la  well-known  that  field  component#  In  this  region  do  not  have  any  atrauge 
behaviour  and  their  variation  in  polar  region  l#  more  or  lea#  smooth,  baaed 
on  this  fact,  which  la  supported  bv  theoretical  analysis  and  experimental 
data,  several  techniques  have  been  suggested  for  field  evaluation  In  paraxial 
regions.  In  some  cases  where  simple  field  expressions  are  available  outside 
the  paraxial  region,  Huvgen'a  principle  can  be  Invoked  to  find  the  field 
in  paraxial  directions  by  integrating  the  field  over  an  arbitrary  surface 
enclosing  tho  ob|ect.  When  the  asymptotic  formulas  based  on  certain  canonical 
problems  are  involved  ifor  example,  In  CTD  or  our  case).  It  mav  be  possible  to 
modify  the  existing  formulas  bv  studying  the  exact  solution  to  corresponding 
canonical  problems  for  the  special  directions  In  which  the  existing  formulas 
are  not  valid  [«s|. 

Our  problem  Is  to  find  the  far-sone  radiated  field  generated  bv  a 
surface  current  distribution  excited  bv  a point  source  on  a tlnite  cylinder 
in  the  paraxial  region.  Using  the  same  argument  given  in  the  previous  section, 
wo  assume  that  the  current  distribution  on  the  finite  cylinder  is  the  same 
as  that  of  an  Infinite  one  truncated  to  a finite  length.  The  ettects  ot  the 
wedges  are  accounted  for  when  the  wedge  diffracted  Melds  are  added  to  the 
total  field.  In  the  present  case,  as  shown  later.  It  turns  out  that  one  can 
compute  the  field  along  the  axis  (0  * 0,  n 1 bv  a verv  efficient  numerical 

procedure,  based  on  the  fact  that  the  field  behaves  smoothlv  in  the  point 
region,  one  can  utilise  the  field  value  along  the  axis  given  bv  the  method 
discussed  later  and  the  field  quantities  at  the  points  not  verv  close  to  the 
axial  region  obtained  bv  our  asymptotic  formulas  In  some  interpolat ion  scheme 
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the  field  In  the  paraxial  region.  We  will  now  outline  a very  efficient 
approach  for  calculating  the  field  along  the  axis  (6^  ■ 0,rt). 

The  far-zone  radiated  field  due  to  a current  distribution  on  a finite 
cylinder  is  given  by  equations  similar  to  (3.b).  The  only  difference  is 
that,  in  the  present  case,  the  z- integration  in  (8.6d)  is  carried  out  over 
the  range  [ — 4 ^ 

Referring  to  Fig.  3a  in  the  axial  region  6.,  * t^,  the  rectangular 

U ’T 

components  of  the  far  field  can  be  written  as 


TT 

A(r)  • / ad<J> 
-tt 


dz  JxU,z)e±Jkz 


(9. la) 


A(r)  • [ adi  / ~ dz  J ( i> , z ) 

, y 

-»  -l1 


s±Jkz 


19. lb) 


where  A(r)  is  defined  in  (8.6c),  and  and  are  rectangular  components  of 
the  surface  current  which  can  be  given  in  the  term  of  J , 


J "_J  sin^  , J * J cosit 

X -t  V it 


(9.2) 


Equations  (9.1)  have  been  obtained  by  evaluating  (8.1)  along  the  axis. 

It  follows  from  (9.1)  and  (9.2)  that  only  the  ^-component  of  the  surface 
current  contributes  to  the  field  along  the  axis. 

The  exact  modal  solution  for  the  surface  current  J ^ due  to  a 
magnetic  dipole  of  moment 
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Insertion  of  (9.4)  into  (9.1)  and  some  algebraic  manipulation  leads  us 


/ dk  V(k  +k) 
J z z 


k Hj2)(kta) 


(2) ' 

ktHl  (kt3) 


(9.5a) 


kt  Hi2)(kta) 


sax  - — | . S • / dk  v(k  +k)  • *- 

y 8rr2  r -i  2 2 H 2)  (k  a) 


(9.5b) 


where 


-jti  +Jti 

V(t)  - / 2 e-JtZ  dz  - * =S 


(9.5c) 


(minus  sign  is  for  0^  * 0 and  plus  sign  is  for  0^  ■ jt).  It  can  easily  be 


shown  that  the  above  integrals  are  very  slowly  convergent  in  many  cases  , 


the  fact  which  makes  their  accurate  evaluation  a tedious  task. 


In  this  work  we  have  employed  a technique  used  by  Duncan 


[109]  in  the  study  of  cylindrical  antenna  problems.  Here  the  method  is 


illustrated  by  discussing  its  application  to  (9.5a).  For  (9.5b),  only  the 


final  results  are  given. 


In  order  to  apply  this  technique  to  (9.5a),  let  us  rewrite  the  integral 


in  (9.5a)  in  the  following  form: 


where 


'•'Then  - k.  it  can  be  shown  that  F - -ka.  The  function  l' 


written  as  follows 


„ 


134 


+j  > J dk 

* 1 


cos(k  -k)i  - cos(k  — k) 2. , 

Z 2 2 1 

k -k 
z 


F (k  ) 
2 z 


(9.10) 


From  (9.9),  one  can  derive 


x 00  sin(k  -k)4»  - sin(k  -kH 

Re  C - Re  { -j  / dk  [ * i k]  F(k  ) 

0 z k2  z 


J(k  -kH  j(k  -k )i 
r ' 2 1 


+ ' dk 


- e 


k -k 
z 


• F(k  ) } 
z 


(9.11) 


The  first  term  on  the  right-hand  side  of  (9.11)  can  be  computed  by  (9.6b), 


* sin(k  -k)l  - sin(k  -k)i. 

Re  { -j  / dk  [ ^ _2_ * k 

0 z kz*k 


] F(kz)  } - 


k sin(k  -k) i.  - sin(k  -k )i, 

j dk  ? 2 z 1 

* -T 


k -k 
z 


F2(kz) 


(9.12) 


In  the  second  term  of  (9.11)  given  below 
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Re  { / dk 


J(k  -k)*2  j(kz-k)H1 

e - e 

k -k 
z 


F(k  ) 
z 


(9.13) 


the  k -integration  is  along  the  path  in  the  complex  kz~plane  (Fig.  35). 

It  can  be  shown  that  the  integrand  has  no  pole  in  the  first  quadrant  and 

also  the  integrand  exponentially  decays  as  |kz|  ■*  »,  as  long  as  Im  kz  > 0. 

Therefore,  if  the  integration  inside  the  curly  brackets  of  (9.13)  is  carried 

out  along  a cl'osed  contour  composed  of  C^,  a quarter  circle  at  infinity  in 

the  first  quadrant  going  from  Re  k to  Im  k , and  C„  (Fig.  35),  the  result 

z z j 

is  zero.  Consequently,  the  integral  (9.13)  can  be  evaluated  along  the  path 
C,,.  After  proper  change  of  variable,  (9.13)  can  be  rewritten  as 

. -Jk*2-U2 

*•  / dc  • 5 ’ 


Se  1 Ip<*2)  - W 


(9.14) 


where 


iF«> 


f dC  ■ Ifcfosk;  + J.kco,k;  - k.l,UH)|«-t;  . [-1FH;)  1 (9.15a) 

0 c +k" 


lp(0  - I“(0  (k  -*  -k} 


(9.15b) 


Combining  (9.15),  (9.13)  and  (9.12),  we  have 


j 
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U (k  ) 


y 2 


sin(k  -k)l  + sin(k  -k)l, 

S Z 1 

k -k 
z 


sin(k  +k)l,  + sln(k  +k)4, 
z 2 z 1 

k +k 
z 


+ j 


cos(k  -k)l  - cos (k  -k )l 
~ - z 1 

k -k 
z 


cos (k  +k) l - cos (k  +k)  l , 

2 •- Z 1 

k +k 
z 


almost  in  Che  same  manner  as  discussed  earlier.  We  have 


I 

v 


!k| 

/ dk 

o'  2 


-J(k  -k)l„ 


k -k 
z 


+ 


-j(kz+k)£ 

k +k 
z 


G,(k  ) dk 

- Z 2 


+ Im  [IgU2)  + VV  + IGU2)  + Xg<V  1 
+ JRe  fIG(i2)  ~ IG(il)  ' IGU2)  + 1 


where 


1^(0 


d;  Ucoskl-ksinkQ  + 1 (kcosk.£+;sinkc ) 


. 2 
k +C 


G(  j s)e 


-Cl 


IgU)  - lj(t)  (k  - -k> 

k^ 
k 

z 


G(k  ) 
z 


J + j G ,(k^ ) . k > 

. k < 


W 


(9.18c) 


(9.19a) 


(9.19b) 


(9.19c) 


(9.19d) 


I I 


' 
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G (k  ) - k 
1 z t 


Vkta)Ji(kta)  + Vkta)Yi(kta) 
Ji2  (kca)  + Yl2  (kta) 


(9. 19e) 


G2(kz) 


_2_ 

tra 


1 

J^2(kca)  + Y^2 (k^a) 


(9. 19f ) 


G3(kz)  * "lkt' 


KjTTk^TI) 


(9. 19g) 


Formulas  (9.17)  and  (9.19)  have  been  derived  from  0^  » 0.  The  same 
formulas  can  be  applied  to  0^  = t,  after  replacing  k by  -k. 

For  numerical  computation  the  expressions  derived  for  the  integrals 
I and  1^  are  much  more  convenient  than  the  original  ones.  The  new 
expressions  (9.13)  and  (9.19)  consist  of  two  parts.  The  first  part 
is  a finite  integral  with  an  oscillatory  integrand.  Because  of  finiteness 
of  the  integration  domain,  many  efficient  algorithms  exist  to  perform  the 
task.  The  second  part  is  a linear  combination  of  infinite  integrals  whose 
integrands  are  fastly  decaying  because  of  the  presence  of  the  factor  e s'", 
therefore,  their  numerical  evaluation  can  be  easily  accomplished.  It  is 
observed  that  from  computational  points  of  view,  (9.18)  and  (9.19)  are 
much  more  suitable  than  (9.6a)  and  (9.13a). 

Once  1^  and  1^  are  determined,  the  field  components  can  be  obtained 

from: 


8n2a 


(9.20a) 


I 


V 


(9.20b) 
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It  should  be  mentioned  that  the  above  formulas  have  been  derived  for  a finite 
cylinder  and  therefore  the  truncation  effect  has  already  been  taken  Into 
consideration.  The  only  corrective  term  which  should  be  added  to  (9.20)  Is 
the  wedge  diffracted  field  which  has  been  discussed  in  Chapter  7. 


141 


10.  NUMERICAL  RESULTS 

A computer  program  for  the  calculation  of  the  tar-zone  radiated 
field  has  been  developed  based  upon  the  expressions  derived  in  the  previous 
sections.  The  program  first  computes  the  main  constituents  of  the  total 
field  (surface  rays,  wedge  diffracted  field,  and  truncation  effect)  and 
then  sums  them  up  to  obtain  the  final  value  for  the  field.  Special  situations 
such  as  field  evaluation  in  paraxial  and  caustic  regions  are  being  handled 
properly  by  the  program  itself.  The  complete  listing  of  the  main  program 
with  all  of  its  subroutines  and  some  additional  clarifying  remarks  has 
been  given  in  Appendix  D. 

Extensive  numerical  results  have  been  obtained  by  this  program.  Only 
some  typical  results  for  various  cases  are  presented  here.  The  geometry 
chosen  for  the  numerical  analysis  is  characterized  by  the  following 
parameters : 

a - 2\ , t - b\,  i2  - 2A  (10.1) 

All  of  the  lengths  are  normalized  with  respect  to  the  wavelength  (\). 

The  numbers  given  in  (10.1)  indicate  that  the  sample  geometry  is  a circular 
cylinder  of  radius  2\  and  length  8\,  on  which  a point  source  is  placed  2\ 
below  the  upper  wedge.  The  E-plane  and  H-plane  far-radiated  field  patterns 
have  been  generated  for  different  polarizations  of  an  elemental  dipole  of  unit 
moment,  ;M  ■!.  As  far  as  the  polarization  of  the  source  is  concerned,  two 
cases,  the  circumferential  magnetic  dipole  (N-0)  and  the  axial  one  (M  «0) , 
have  been  studied.  It  is  abvious  that  any  other  polarization  of  the  source 
can  be  considered  as  a linear  combination  of  the  two  previous  cases.  In  each 
case,  the  graphs  (Figs.  )b  to  sit  illustrate  the  variation  of  the  amplitudes 
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Ej  versus  dQ (degrees)  when  the  source  is  a 
circumferential  magnetic  dipole  and 


a-2X 


for  different  values  of 
a)  ^0-0°,  b)  *0-3(T,  c)  *0 
d)  .i>0-90°  e)  t0-120\  f) 


(Fig.  19) 
-60°, 
-160°. 
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150 


gure  38.  j versus  ^(degrees)  when  the  source  is 

a circumferential  magnetic  dipole  and 

V6*’  l2m2X'  a"2x 


for  different  values  of 


(Fig.  19) 


53 


Figure  39. 


E 1 versus  $ (degrees)  when  the  source  is 
P ' 0 

a circumfc  rential  magnetic  dipole  and 


2A,  a=2.\ 


for  different  values  of  9 


0- 


(Fig.  19) 


a)  9 =30°,  b)  9 =50°,  c)  6 =70°,  d)  9 =80 
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F 

Figure  40.  | | versus  9Q  (degrees)  when  the  source  is 

an  axial  magnetic  dipole  and 
a=2A  for 

different  values  of  *Q.  (Flg-  19) 

a)  $Q=10°,  b)  $q=20° , c)  4>0=50° 

d)  $Q=70°,  e)  $Q=80°,  f)  d>0=90° , 
g)  $0=120o,  h)  $0=150°,  i)  $0=170°. 


|E*|  <<»> 


c90  (degrees) 

(a) 


Figure  -*1.  i [ versus  (degrees)  when  the  source  is 
an  axial  magnetic  dipole  and 
i,«2\,  a-2X 

for  different  values  of  $ . (Fig.  19) 
a)  V20°,  b)  , -40* . c)  f -70- . d)  * -1. 


168 


r 


0 50  ICC  I5C  2C0 

<p  (degrees) 


Figure  43.  |E  j versus  when  the  source  is 
an  axial  magnetic  dipole  and 


*1-6X,  i2- 2\,  a»2\ 

for  different  values  of  6^.  (Fig.  19) 
a)  t?0*30° , b)  *0-50°,  c)  90-120°. 


170 


171 


of  the  two  spherical  components  of  the  radiated  electric  field  with  respect 
to  <}> q (for  different  0Q  cuts)  and  0Q  (for  different  cuts).  The  amplitude 
of  the  various  components  of  the  field  have  been  given  in  dB". 


I E0 | in  dB  = 20  log10 ( I E0 | ) 


(10.2a) 


|Ej  in  dB  = 20  log1()(|Ej) 


(10.2b) 


Figs.  36  to  39  show  the  radiation  patterns  of  the  spherical  components 
of  the  electric  field  on  4>q  = constant  planes  and  9^  = constant  Cunes,  for 
the  case  where  the  elemental  magnetic  dipole  has  circumferential  polarization 
(M  = 0).  The  variations  of  the  dominant  component  | E_ | and  nondominant 
component  (E^j  have  been  shown.  It  is  observed  in  the  lit  region  (jbg|<  y) 
that  the  amplitude  of  the  nondominant  component  is  about  10  or  20  dB  below 
the  dominant  one.  But  as  the  shadow  region  is  approached,  the  difference 
tends  to  decrease.  Lack  of  smoothness  of  the  plots  is  mainly  due  to  the 
finiteness  of  the  number  of  points  in  each  curve.  Of  course,  the  fine  struc- 
ture of  the  curves  can  be  studied  by  increasing  the  number  of  the  points.  In 
all  the  cases  in  the  deep  lit  region  ( | $ | <45°  and  0^  ^ 90°),  the  basic  struc- 
ture of  the  plots  can  be  identified  with  those  of  the  y-directed  magnetic 
dipole  on  a flat  ground  plane.  In  this  region,  the  geometrical  optical  field 
is  a reasonable  approximation.  In  the  paraxial  and  shadow  region,  the  more 
accurate  results  deviate  substantially  from  geometrical  optical  predictions. 

In  these  regions,  wedge  diffracted  and  surface  diffracted  rays  contribute 
significantly.  These  latter  contributions  are  also  responsible  for  too  many 
fluctuations  observed  in  the  field  amplitude  specially  in  the  deep  shadow. 
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Special  narrow  peaks  near  the  ■ 180°-end  of  the  graphs  of  the  field  com- 
ponents versus  (Figs.  38  and  39)  are  mostly  due  to  the  formation  of  a 
caustic  region  for  the  wedge  diffracted  fields  along  these  directions  (the 
detailed  discussion  of  this  case  was  given  at  the  end  of  Chapter  7) . The 
formation  of  the  caustic  generates  a relatively  strong  field  in  the  corre- 
sponding directions  in  the  shadow  region.  The  strong  peaks  at  the  polar 
regions  in  Fig.  37  indicate  the  significant  deviation  from  the  G.O. 
(Geometrical  Optical)  and  even  the  exact  solution  for  the  infinite  cylinder 
when  paraxial  regions  are  approached.  Once  more  it  should  be  emphasized  that 
the  curves  show  a rough  estimate  of  the  variations.  The  other  point  which 
should  be  mentioned  is  that  for  the  field  amplitude  below  some  certain 
level  the  computational  error  can  exceed  the  value  of  the  results. 

The  precise  scudy  of  our  extensive  numerical  results  indicates  that 
for  the  specific  geometry  chosen  for  computations  the  truncation  effect 
can  amount  to  up  to  20*  of  the  final  result  in  some  cases.  Of  course, 
as  and  C.,  increase  this  effect  decreases. 

Figs.  -*0  and  43  illustrate  the  field  variation  with  respect  to  $ 

and  when  the  source  has  a z-polarization  (M  »0) . The  general  properties 

U 4* 

described  in  the  previous  case  can  also  be  observed  in  these  curves.  The 
behavior  of  the  dominant  component  of  the  field  E in  the  deep  lit  region 
is  almost  similar  to  what  can  be  predicted  by  G.O.  expressions.  In  this 
region,  E^  (Fig.  41)  is  much  weaker  than  the  dominant  polarization.  Again, 
when  we  approach  the  shadow  region,  the  patterns  deteriorate  and  significant 
deviation  from  G.O.  is  observed  (Figs.  40f,  g,  h,  i).  The  field  variations 
versus  are  depicted  in  Figs.  42  and  43.  In  the  deep  lit  region.  ^Fics. 

42c,  d,  e,  f ) , E patterns  are  quite  analogous  to  those  of  a magnetic  current 

$ 
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filament  on  an  infinite  cylinder,  and  the  finiteness  of  the  object  affects 
only  the  shadow  part  of  the  plots.  In  this  case  too,  like  in  previous  ones, 
a high-order  caustic  formation  is  responsible  for  a fairly  sharp  peak  towards 
the  <J>q=180°  end  of  the  plots. 
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11.  SUMMARY  AND  RECOMMENDATION  FOR  FUTURE  WORK 

In  Che  beginning  of  this  work  (Chapters  2 and  3)  we  have  examined  the 
problem  of  radiation  from  sources  in  the  presence  of  smooth,  convex, 
impenetrable  objects.  First,  a brief  survey  of  various  high  frequency 
techniques  was  presented.  A generalization  of  the  geometrical  theory  of 
diffraction  and  two  new  techniques  based  on  the  spectral  domain  approach 
and  an  asymptotic  evaluation  of  the  radiation  integral  for  the  surface 
current,  were  discussed.  It  was  shown  that  is  is  possible  to  extend  the 
range  of  applicability  and  to  improve  upon  the  accuracy  of  the 
existing  asymptotic  theories  by  solving  new  canonical  problems  in  which 
some  of  the  earlier  restricting  assumptions  underlying  GTD  and  other  high 
frequency  techniques  have  been  removed.  The  accuracy  of  the  new  approach  was 
evaluated  by  comparing  the  numerical  results  derived  from  the  new  formulas 
with  the  available  theoretical  and  experimental  data. 

In  the  second  part  of  the  present  work  (Chapters  5 through  10),  the 
special  case  of  source  radiation  in  the  presence  of  a finite  solid  cylinder 
was  studied.  For  the  problem  of  high  frequency  radiation  of  an  electro- 
magnetic point  source  on  the  surface  of  a conducting  circular  cylinder  with 
finite  length,  a technique  combining  the  main  features  of  different  asymptotic 
theories  has  been  presented.  The  STD  interpretation  has  been  employed  to 
derive  the  diffracted  field  and  the  total  radiated  field  has  been  obtained 
by  adding  the  contribution  of  the  "body"  currents  and  the  field  scattered  by 
the  wedges.  A modified  version  of  GTD  has  been  applied  to  wedge  diffraction. 

The  modification  consists  of  removing  the  singularity  of  Keller’s  diffraction 
coefficient  which,  according  to  STD,  can  be  associated  with  the  field 
radiated  by  a current  distribution  with  a semi- inf ini te  support  on  the 
faces  of  a straight  wedge.  The  incident  surface  field  at  the  wedge  has 

I 
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been  obtained  by  asymptotic  expressions  derived  by  Lee  and  Safavl-Naini  [28]. 
The  contributions  of  "body"  current  have  been  determined  by  approximate 
asymptotic  formulas  developed  by  Safavi-Naini  and  Mittra  [20].  The  effect 
of  finiteness  of  the  cylinder  has  also  been  studied  and  accounted  for 
in  the  final  expressions  for  the  total  field.  The  special  situations  where 
the  field  evaluation  in  the  first-  or  higher-order  caustics  are  involved 
have  been  treated  successfully  by  introducing  the  notion  of  "non-uniform 
equivalent  edge  current"  [46]  and  modifying  the  GTD  in  a proper  manner. 

Based  upon  this  new  formulation  of  the  problem,  a computer  program  has 
been  developed  which  can  achieve  field  computation  for  all  observation  angles 
taking  into  account  all  various  effects  or  specific  situations  whose  hand- 
ling requires  special  care. 

The  accuracy  of  the  procedure  has  been  partially  justified  by  the 
demonstrated  validity  of  the  different  theories  used  in  the  present  method 
when  applied  to  other  similar  problems.  Finally,  it  is  recommended  that 
a thorough  testing  of  the  accuracy  of  the  procedure  outlined  in  this  paper 
be  carried  out  in  order  to  fully  evaluate  its  usefulness.  This  can  be 
done  by  following  a procedure  suggested  in  [110]. 
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APPENDIX  A:  FOCK  FUNCTIONS 

In  studies  of  radio-wave  propagation  around  the  earth  by  Van  der 
Pol,  Bremmer,  Pryce,  Fock,  and  others,  and  also  the  later  studies  of 
diffraction  of  electromagnetic  waves  by  certain  bodies  of  revolutions 
([51],  [52],  [53],  [54],  [55],  [48],  [56],  [57],  [58],  [59],  [60]  and 
[50]),  a class  of  universal  functions  was  introduced  which  can  be  used 
to  predict  the  amplitude  and  the  phase  of  the  reflected  or  diffracted 
field  by  smooth  convex  surfaces  [17].  An  exhaustive  treatment  of  these 
functions  which,  in  general,  are  defined  as  Fourier  integrals  having 
combinations  of  Airy  integrals  in  their  integrants,  has  been  carried  out 
by  N.  A.  Logan  [61].  (See  also  Bowman,  et  al.  [1]  and  Logan  and  Yee 
[17]. 

Since  the  first  extensive  application  of  these  functions  to  diffrac- 
tion theory  was  done  by  Fock,  many  authors  named  them  after  him.  Here 
we  list  only  the  most  important  formulas  and  expressions  for  these 
functions  without  going  through  the  details  of  their  derivations.  We 
have  followed  Logan's  set  of  notations  for  these  functions  [61].  However, 
since  his  time  dependence  factor,  exp(-iwt),  is  different  from  one  we 
have  used  throughout  this  paper,  namely  exp(+ja)t),  our  expressions, 
listed  below,  are  conjugates  of  what  have  been  presented  in  [61]. 

First  we  start  with  general  definitions.  Fock's  most  general  form 
of  the  "Van  der  Pol-Bremmer  diffraction  formula"  is 


w,(  t - y^  • | v(t  - v<) 


n 
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where  w^(t),  w^(c),  u(t)  and  v(c)  are  Fock-type  Airy  functions,  defined  as 


u(t)  » >•  tt  Bi(t) 


, v(t)  = /n  Ai(t) 

* 


w1<t)  » u(t)  + jv(t)  , w .^ ( t ) » w^t) 

We  note  that  w^  and  w,  can  also  be  defined  as  in  Sec.  2.  y>  and 
y<  are  the  larger  and  smaller  of  the  two  numbers  v ^ and  yn.  V(x,y^,v,,q)  is 
proportional  to  the  attenuation  suffered  by  an  electromagnetic  wave  generated 
by  a source  located  at  reduced  height  y^  above  the  surface  of  a smooth  con- 
vex body,  when  it  reaches  the  observation  point  located  at  reduced  height 
y,  above  the  same  surface.  x is  the  reduced  distance  between  the  source  and 
the  observation  point  along  the  surface,  and  q is  dependent  upon  the  impedance 
of  the  surface.  Let  us  consider  some  useful  limiting  cases. 

When  » y7  ■ 0,  then  V(x,0,0,q)  is  denoted  by  V^,  where 


jw/4 


Jxt 


VQ(x.q)  - 2 


_ 4 , — r » e J ' w (t) 

e _ . * . f 5 . dc 

J*  J_ oo  w2(t)  “ <’w2(t) 


(A.  2) 


We  also  have 


jff/4  r-  , - e_jxC  w,(t) 

v00  . V„.0)  - — fi  j_  - v‘(d  ■ 


(A.  3) 


f 

u(x)  • lim  |-2jxq“Vri(x,q) 


When  y^  = 0 and  v-,  -*■  °°,  then  V •*  V^(x,q): 


V^(x,q) 


I ej  3,1/4  3/2  f * e W2(C) 


• x 


w^(c)  " qw,(t) 


J-a 


W,(t) 


dt 


(A. 4) 


U.  >> 


and  also 


g(x)  * V^(x,0)  * — 


. ? 

wi(t) 


dt 


y «T  * — OD 


(.A.bt 


■xsr.:^ 
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f(x)  * lim 

q-M» 


f 

qV1(x,q) 


_1_ 

v'tt 


f " e'jxt 
i_=o  w2(c) 


dt 


(A.  7) 


A.l.  Functions  f and  g 
n ®n 


Based  on  Equations  (A. 6)  and  (A. 7),  a class  of  functions  can  be 
defined: 


f(n)(x) 


(-.1)1  . 
Jn 


r w2(t)  dxn 


(A.  8) 


g(n)(x) 


IzJL)  . 

A 


r tn 


-jxt  n , . 

-T/-x  * ^ 

W2(t)  dxn 


(A. 9) 


where  T is  any  path  in  the  complex  t-plane  which  comes  from  in  a sector 

TT  TT  TT 

defined  by  - it  <_  argit)  < - y and  goes  to  +°°  in  the  sector  - -j  < arg(t)  < 

In  what  follows,  we  will  give  the  suitable  formulas  for  f (x)  and  g(x)  in 
different  ranges.  Tabulated  values  and  graphs  of  these  functions  can  be 
found  in  [57],  [54],  and  [61]. 

When  x is  very  large  and  negative,  the  following  asymptotic  expansions 
for  f(x)  and  g(x)  can  be  used  [61]. 


f(x)  . 'V  ! + -L  ♦ 1 JV5  J2f.  + J3i«U5+  ... 

4x  2x  64x  16x  “ 1024x 


(A. 10) 


g(x)  - 2e 


jx3/3 


1 - 


1 , j 469  , 5005  j 1122121 

' r • “ , i.  " 


4x 


6 */  9 ' A/  12 

x 64x  6**x 


1024x 


15 


(A. 11) 


The  above  formulas  are  valid  and  accurate  for  x <<  -1.  For  moderate 
values  of  x,  namely,  -1  £ x <_  1,  it  is  difficult  to  find  an  appropriate 
expression.  Although  there  are  some  analytical  techniques  like  "stationarv 
phase  method"  or  "Poisson  summation  formula"  which  mav  be  used  to  evaluate 
f^n^  and  for  these  values,  another  possible  way  which  is  probablv  easier 

and  more  efficient  is  to  interpolate  the  tabulated  values  of  these  functions 


in  this  range. 
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In  the  vicinity  of  zero  (|x|  = 0) , the  Taylor  expansion  can  be  used 
to  calculate  f and  g.  The  coefficients  are  given  by 


r(n) (0) 


j (.  5n7r /6— tt/ 3) 


3tt 


(2/ 3) (n-1/4)  *> 

• I A (n) 
m»0  m 


i , 1m 


2 m — 


4n  - 1 


g(n>(0) 


-j5mr/6  r- 
e • v tt  • 


4n  - 3 
Tfm  ‘ 6 


✓ 7T  • 

|) 

•+) 


■>_)  (2/ 3)  (n-3/4)  00 

f ■ l Bm<n)  • TT 

2'  m=0  m l3ir‘ 


11 

* M 


i> 


where  x(\,u)  is  the  generalized  "tau"  function: 


T(X.U)  - y , \ > i 

n=0  (n  + u) 


A0(n) 


A-  (n) 


Vn) 


B,(n) 


1 , A^n)  = Tg  (n  - 1) 


5[5„-  - - 32) 

l lb  J 


B:(n)  - -7 (n  - 3/2)/48 


(49n2  + 364n  + 39849/16) / (29  • 3 2) 


(A. 12) 


(A. 13) 


(A. 14) 


When  x is  large,  and  positive,  residue  series  can  be  used  to  compute 


,(n)  , (n) 

t 3nd  g 


f(n)(x)  = ej(2+7n)-r/6  y 

p-1 


< ,n  , -j57T/6. 

(r  ) exp(r  • x • e J ) 
? £_ 


Ai'<-  v 


i A . 15) 


di M 
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g(n>(x)  , ej7,Tn/6  ? (rP>""1  eXP(rp  * * • e~j5n/6) 

p-i  Ai(- 


CA. 16) 


where  Ai(-  rp)  - 0 and  Ai'(-r^)  = 0 for  p = 1,2,3,... 

A . 2 . Functions  u,  v , u ' , v 1 : 

Fock's  functions  u and  v defined  in  (A. 3)  and  (A. 4)  may  be  re-expressed  as. 


v<«  - i ± ( ^ .-«« de 

/tt  J r 2y  j 


(a. 17) 


■«)  - .13,/V/2  A-  f .-«*  dt 

K if  W2(C) 


(A. 18) 


where  integration  contours  V ^ and  T , are  sketched  in  Figure  44,  and  w^(t)  is 
the  derivative  of  w?(t). 

Zeroes  of  w^(t)  and  w!,(t):  They  are  given  by 


tn  = !rn!e-J t/3  , and  t^  = | r' 


n n 1 


(A. 19) 


respectively.  The  magnitudes  of  the  first  ten  zeroes  are  listed  in  Table  A.l. 


A. 2. 2 Residue  series 


representation.  For  real  positive 


V(C)  - J1'2  l J"« 


-j5t' 


(A. 20) 


u<£)  « ej7r/42/7  £ 3/ 2 l e~JCCn 


(A. 21) 


(£)  « j e“j7T/4/n  £_1/2  ? (1  - j2£t’)(t,r1e  ^ n 

n^l  n n 


(A. 22) 


»’(£)  *e^/43^£1/2  l (l  - j j£t]e'j'C"  . 

n=l  n 


(A. 23) 


6 GO 


Figure  44.  Contours  and  on  the  complex  t (or  z)  plane.  Ti,  for 
example,  goes  from  00  to  0 along  the  line  Arg  t = — 2rr/ 3 and 
from  0 to  00  along  the  real  axis. 
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TABLE  A.l 

The  modulas  of  the  zeroes  of  the  functions  w^Ct)  and  w^(t) 


TABLE  A. 2 

*%  ” ;1  - (Residue/stnall  arg.)|  x 100 

j Difference  at 


Mag.  (X)* 

Phase  (deg.) 

0.00 

0.00 

0.11 

0.01 

0.02 

0.08 

0.09 

0.15 

0.10 

0.90 

A. 2. 3.  Small  argument  asymptotic  expansion.  For  real  positive  £ and 


£ - 0, 


»«>  „ 1 - £ .J"'V'2  m ■ 10- v * . . . 

(A. 2 


u(0  ^ 


1.jj.i-/V/2,3:sit^e-i'/v'2.J. 


701  * 10“2C6  + ... 


V.(M  % ~J  3n/4  1/2  + ?!  2 63*5  e-J»/4  7/2  _ 

v 8 ^ 20  ^ 1024  e 4 


485  * io~:£5  + . . . 


-•«)  - $ * .“j3nMC1/2  + % £2  + ^ e^'V'2  _ 2>2 


21  * 10-1C5  + ... 


A. 2. 4.  Numerical  evaluation.  For  £ ^ £Q,  the  residue  series 
representation  with  the  first  ten  terms  in  the  summation  may  be  used. 
For  £ £Q,  the  small  argument  asymptotic  expansion  with  the  first 

five  terms  may  be  used.  It  can  be  shown  that  the  smoothest  crossover 
is  obtained  if  £q  - 0.6.  In  the  present  study,  we  set  £^  - 0.7,  where 
the  difference  in  the  two  representations  shown  in  Table  A. 2. 
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APPENDIX  B:  DERIVATION  OF  FORMULAS  (3.8)and  (3.9) 


Here,  we  consider  only  the  derivation  of  the  asymptotic  expansion  of 
$ for  a circumferential  magnetic  dipole.  In  this  case,  $ may  be  written 


as : 


where 


JweM 

£ 

(2tt)2 


;(kt>  - i 

n--“ 


dk 


- j k z 

e Z • S(kt) 


(2) 

. . Hv  '’(k  p) 
,-jn*  . n tM/ 

<2,V> 


Applying  the  Watson  transformation  to  (B.2), 


S(kt)  - 2 


1 . 


lli“><ke*»  e~J  V<^-7T) 


,.. ..  u(2)  ,.  . sin  vtt 

C+D  H (k  a) 

v t 


dv 


where  C and  D are  shown  in  Fig. 45  . Or, 


S(kt)  - j • 


i(2) , 


I — »— J6  cosv(*  - $)  H^‘,/(k  o) 

JtS  sinvTi  • H^“^(k  a) 


dv 


Substituting  the  expansion 


COS  V(T  - i)  ( y y "J  V C ^ +2  IT  1 ) 

sin  vrr  * e 


i-1  1-0 

where  4> L - 4>  and  - 2 n - $,  in  (B.4),  the  result  will  he: 


!(kt>  » l l 


-j^-'(kto)  -jv(*  +2*1) 

"77) * e • dv 


i-l  1-0-'— -J6H  (k  a) 

J v t 


(B.l) 


(B.2) 


(B.  3) 


(B.4) 


(B.  5) 


( B . 6) 


Each  term  of  the  above  expansion  is  associated  with  a "creeping  wave* 
travelling  in  a counterclockwise  (i  - 1)  or  clockwise  (i  - 2)  direction 
around  the  cylinder.  Following  the  ra”  concept,  each  creeping  wave 
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appears  to  be  travelling  along  a specific  surface  ray.  Now,  as  e * ® 
(far  rone)  for  each  fixed  v,  we  have  l°-l 


H(2)(k  ■>)  * f-rr—  • 
v t ^ wkto 


-J  lkto-vn/2-’»/4) 


(B.7) 


On  the  other  hand.  It  can  be  shown  that  the  significant  contribution 

to  Slk^)  comes  from  a small  neighborhood  of  kfa.  In  this  neighborhood , where 

1/1 

k^a  and  v are  large  and  close  to  each  other  (|k  a - v|  < |v|  ),  the  Hankel's 

asymptotic  expansion  (B.7)  is  not  valid  any  longer.  In  this  case,  it  is 
necessary  to  expand  Bessel's  functions  in  terms  of  Foek-type,  Airy  functions, 
w (t)  and  w,(t),  and  their  derivatives  lib] 


HU^  (x)  , — ^ /w  ,(t) 1 4 1 w,(t)  + t*  w (t)l  + ... 

m>  a J bOm*  l ' | 


H (x)  - v',(t)  -f  c 1 4 1 w', it)  + (o  - t ) w.,(t) 

m"n  “ bOm"  ^ 


IB. 3) 


v B.  9) 


where 


> 1/3 

x 1 


V w X 

t - (m  is  verv  large) 

m 


Inserting  (B.")  and  the  first-order  terms  of  iB.SI  and  lB.91  into  tB.b> 


and  IB.  1),  we  obtain 


(h  Jt/4 


2 JO  vV 

y y 1 dk 

i-l  1^0  2 


'il  . . . . m* 

• Vni^  • ~T73 

kt 


iB. lOi 


J 


- ilc  a 2) 


il  • V * V* 


*11  • 9{it  * -,'1 


a(:.  * 2t l 


■ 
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Introducing  a new  integration  variable  a: 


k - k sin  a 
z 


k£  • k cos  a 


and 

we  have: 


au  - tan"l{z/lP  + a(*t  + 2^1  - ir/2)] 


(B. 11) 
(B. 12) 
(B. 13) 


• kR^  cos(8li  - a) 


(B . 14) 


where 


Ra  -|z'  + (p  + a(^L  + 2 w 1 - ir/2)ij 


1/2 


Now  (B.10)  takes  the  following  form: 

9 

u)cM . 


♦ - • & • -3y  • .J*'*  • I l 


(2t)“  V 0 kJ/2 


l-l  1-0 


da  • e 


-jkR  jCOSia-B^) 


cos“5/6  a • f0Uu> 


(B. 15) 


y is  the  path  of  integration  in  the  complex  a-plane,  which  is  shown  in 
Fig. 46. 

Now  we  deform  the  path  of  integration  into  the  "steepest  descent  path," 
SDP,  passing  through  the  saddle  point  of  the  phase  of  the  integrand.  Per- 
forming the  "saddle-point  integration,"  we  can  derive  the  asymptotic 
expansion  of  <. B . 15)  for  large  k.R ^ j • The  first  order  term  is: 

1/3  ■>  .»  “JRB, 

12  j 


acM 


. ‘f  . J’r/2  |ka)  r r , , . 

> ' 7 e • Tj  * > l (cos  8,,) 

2ik“  1 ‘ ’ i-1  1-0 


-4/3  e 


lls 


ils 


t0Uils) 


(&.  lo) 


where  R. , _ and  are  the  values  of  these  parameters  at  the  stationary  point 


specif ied  by  a - £ , . 


J 
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Equation  (B.lb)  is  the  creeping-wave  representation  of  the  far  field. 

If  the  cylinder  is  large  (ka  >>  1)  and  | t \ is  not  very  close  to  i , then  onlv 
the  first  term  (i  • 0,  1 ■ 1)  has  the  most  important  contribution  to  the 
total  infinite  sum,  and  the  other  terms  are  not  significant.  Neglecting 
the  other  terms,  we  obtain  the  result  given  in  (3.8)  and  (3.9).  It 
should  be  emphasized  that  (3.8)  and  (3.9)  are  not  valid  when  |3|  is  close  to 
t/2  (paraxial  region),  because  in  this  case,  kta  is  very  small,  and  (B.7), 
(B.8)  and  (B.9)  no  longer  apply. 

The  other  formulas  can  be  derived  in  a similar  manner. 
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APPENDIX  C:  ASYMPTOTIC  EVALUATION  OF  THE  RADIATION  INTEGRAL 


Consider  the  following  double  integral 


jk*(x,y) 


where  g(x,y)  is  rather  slowly  varying,  and  $(x,y)  has  a stationary  point 
(x  »y  ) inside  domain  D.  The  objective  is  to  derive  an  asymptotic  expansion 
for  (C.l)  when  k is  large. 


Suppose  g and  } have  the  following  forms  around  (x  ,y  ) 


N.  Chako 


has  derived  the  following  asymptotic  series  for  U 


where 
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g (x.y)  - l g.  ,(x  - x )k(y  - v J1 

k,l-0  S 


P(x,y)  - V a (x  - x )m(y  - y )n 
*•  , mn  s s 

on-n>  1 


Q(x,y)  - l b (x  - x )m(y  - y )n 
. “ , mn  s s 

m+n>l 


a10  , bio] 


'NjO  ‘ 800  ’ A10  " 810  800(V0  + l)  5 + t 

a ;f01  . Vl 

A01  " S01  800^  5 (u0  L)  t j 


In  order  to  apply  this  procedure  to  the  integrals  of  the  type  (50)  for 


which 


j(x,y)  - • -(R  + a) 


(C.4) 


g(x,v)  » F(j , 3 ,P)  — 


(C.  5) 


When  F is  one  of  the  components  of  J(1  - RR) , one  should  first  determine 
the  stationary  point  of  (1,  wherein  its  first-order  derivatives  vanish.  The 
second  step  is  to  compute  the  various  order  derivatives  of  (1,  J,  R,  ...,at 
this  point,  and  then  insert  them  into  (C.3).  We  just  give  the  main  formulas 
needed  for  these  derivations. 

Suppose  the  surface  of  the  body,  x(o,3),  is  parametrized  by  a geodetical 

polar  coordinate  system.  As  discussed  previously , in  this  system,  J is  the 

— ^ 

arc  length  of  the  surface  geodesic  connecting  the  pole  Q to  x(a,S),  and  8 is 
the  angle  between  the  geodesic  and  some  fixed  reference  geodesic  at  Q 
( cig. 1' ) 

The  element  of  length  in  this  system  is  given  by 


da*  + C(c,8)  di' 


(C.6) 


192 


Let  us  denote  dx(u)/du  bv  x ; then  we  nave  the  following  set  of  relations 


onst 


const 


is  the  outward  unit  normal  to  the  surface.  Another  quantity  of  interest  is  the 


geodetical  curvature"  < given  by 


Using  the  above  relations,  we  can  derive  the  following  expressions 


which  hold  true  at  the  stationary  point 
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aa2  3a 


_G_  3G/3 a 32G/3a2  _ L^f 

•■>  ' ^ ^ 


(C. 16) 


(C. 17) 


33ft  m 33C/36  3G  3G  _L_  So  66  1 32G 

3g3  2R  36  ’ 3o  ’ 4G  2 3o38 

where  is  the  radius  of  curvature  of  the  geodesic. 

Equations  (C.12)  and  (C.13)  determine  the  location  of  the  stationary 
point.  At  this  point  R - x^,  which,  if  we  introduce  the  ray  concept, 
tells  us  that  the  surface  rays  leave  the  surface  at  the  "point  of  diffraction" 
tangentially.  Equation  (C. 14)  indicates  that  the  stationary  point  is  of 

second  order,  so  that  we  need  higher-order  derivatives  of  the  phase . LJ° , L^° 

8 8 

and  L are  coefficients  of  the  second  fundamental  form  of  the  surface 
evaluated  at  the  stationary  point.  They  are  defined  as 


. aa 


08 


.68 


acr 


X63  * X3 


'3  ’ “ a 6 3 ’ 

Using  the  relationships  given  above,  one  can  find  the  expansion  co- 
efficients g. , , a , b and  A in  (C.3).  Zeroth  and  first-order  terms  in 
Kl  mu  mn  p cj 

(C.3)  give  us  formulas  (3.39),  (3.40),  and  (3.41). 

A few  remarks  should  be  made  concerning  the  expansion  presented  in 

(C.3).  First  of  all,  (C.3)  is  a doubly  infinite  series;  therefore,  for  each 

fixed  power  of  k 3 a finite  number  of  terms  should  be  summed  up.  The 

coefficients  of  various  terms  in  these  finite  sums,  namely  A 's  , become 

pq 

very  complicated  when  p and  q are  greater  than  0 or  1.  Another  difficulty  with 
this  series  is  that  when  the  stationary  point  of  the  phase  is  of  an  order 


higher  than  1,  the  difference  between  the  order  of  the  successive  terms 
(when  they  are  ordered  according  to  the  descending  power  of  k)  becomes  very 
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small,  and  consequently  the  infinite  series  converges  very  slowly.  For 
instance,  in  our  problem  where  5-3  and  r - 2 (stationary  point  is  of 
second  order),  sometimes  the  difference  between  the  orders  of  successive 
terms  is  k-^"^,  which  indicates  the  weak  convergence  (in  an  asymptotic 
sense)  of  the  expansion  in  the  cases  where  the  frequency  is  not  very  large. 

i 


i 
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i IS  THE  O.B.  POINT 


Based  upon  the  formulas  derived  in  this  report,  a computer  program 
for  calculating  different  constituents  of  the  far  field  has  been  written. 
The  organization  of  the  program  is  shown  below  and  its  listing  is  given 
in  the  next  pages. 


APPENDIX  D: 

COMPUTER  PROGRAM  FOR  COMPUTATION  OF  THE  FAR  ZONE  FIELD  RADIATED 
BY  A TANGENTIAL  MAGNETIC  DIPOLE  ON  A FINITE  CYLINDER 


INPUTS : 


FINITE  CYL.  - 


INFINITE  CYL.  - UPPER  AND 


IN  PA RAX.  REGION? 


LOWER  SEMIINFINITE  CYL. 


Ai’PROX.  FORMULA 


OUTPUT 


FOR  PARAX.  REGION 


DIFFRACTED  FIELDS 
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The  symbols  used  in  this  program  have  been  listed  below. 

Input  parameters: 

"PHID"  - 4>o(  in  degrees) 

"TETAD"  - 0 (in  degrees) 

"MZ"  - M , "MPHI"  - Ma 
z <p 

"ZO"  - U - i2)/ 2 (in  wavelength) 

"A"  • a , the  radius  of  the  cylinder  (in  wavelength) 

"EWTET"  and  "EWFI"  - E*  and  E* 

"ESTETA"  and  "ESPHI"  - contribution  of  surface  rays 
"EDTETA"  and  "EDPHI"  - contribution  of  direct  rays  (G.O.) 

"TRUNEF"  and  "TRUNET"  - truncation  effect; the  radiated  field  due  to 

the  removed  semiinfinite  cylinders 

"PARAXX",  "PARAXY"  ■ x-  and  y-components  of  the  field  in  paraxial  region 
"ETETA"  and  "EPHI"  ■ 0^-and  ^-components  of  the  total  field. 

All  the  field  components  have  been  normalized  to  the  factor  e J /r. 


MAIN  PROGRAM 


eiaiu 


ft 
)» 


eaidd 

0150 

dlbd 

0170 


"u»’'nn  ,OoTt*UT ) 

0<*iPU:<  eTETA#e*,Ht#trETaU,ti*nlG,r  r£  TA1  ,tHhI  1 ,fc 

rta.iiiftsw* 

xfeAO9tf,TeTa01*»tii0,f'Z,^PHl#2U,C,  A 
90  7GH«A[ (7P7.2) 


rtTA2,£PMIi 


IF  THE  OBSERVATION  POINT  IS  IN  PARAXIAL  REGION,  COMPUTE  THE  FIELD  BY 
INTERPOLATION.  IF  IT  IS  NOT,  THEN  CALL  "FIELD"  SUBROUTINE  TO  COMPUTE 
THE  FIELD. 


33(5:3 


US'.J 
id  i 9J 

0ddl<' 


9«i  7'Um-aT  ItTCM) 

GO  1u  £00 


ll 

► 9 ,4, 5 V , 


,L,  A,trtr*F£f»'il) 
7n  tfcPMlj  «,F9,n) 


THIS  PAGE  IS  BEST  QUALITY  PftiOSflMi 

from  oopy  ruwnsHfi)  ro  ddc  — 
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THIS  PAGt 

ram  spot* 


1U 


UO.Q 


,-KAor; 


FIELD  COMPUTATION  IN  PARAXIAL  REGION  BY  PARABOLIC  INTERPOLATION 


03*30  100  CONTINUE 

00*40  PRINT  U<),  TfeTA0,?nIO,M2, 

03*5U  110  FORMAT  (5X,5HTfcTA*,F7 


MPmI,Z0,L,  A 

2,5X,4MPriI« 


00*6O*|SHMPMl*,F5t2,5X, 3HZ04.F6. 2,5a,2HL»,F6 
00*70  UT*0jBi).4IF(re  rA0,GT,9a,J  Tkl  400*193, 


F7,2,5x, 3mmZ»,F5,2,5X 


6'""r  f •Cii*»',nrn*i/"  r«Cf3Af  jnn4»ir3, 

3HZg«.F6f2,5A.2HL»,F6t2#bX,^rlA.(Fit2j 


00200  print*," 


•PARAXIAL  REGION* 


00*90  PRINT*,"  P I £LQ  COMPUTATION  BY  3-PUInT  INTERPOLATION1 

00300  TETA0i«A3INU,/ (5.*pI*A)  )*iau./Pl*i, 

00310  IA  UtT  AU.UT.90.)  ‘|TAOi«l(  “ “ 

0032a  ]kT A02«TtTACl*5t»lF (TE TaQ, 

0033)  IF C rtTAU.UT.90. J l£TA02«T( 

00340  PRINT*!^  _POIN?^NO.  lr 

i * M I i 


lOtf.-TETADl 

LST.90.)  7£TA0*«TtTA01-5, 
r£TAOl-S, 


[0,MZ,NPtii,ZO,L,  A,ETtTAI  ,tPMll) 

,MZ,MHiH,Z0rL,  A,£TtTA2,£Pnl2) 

■MZ,MphI, ZU.L, A,ETkTAO,tPHlU) 

• I N T feRPOL  A TEO  FltLOA" 
dPHit; 


00350  Call  f JELD CtEtaoi .pmio, 

00360  PwInT*  " PoInT  nO,  2" 

0037O  CALL  FItLUlTETA02,Pnin, 

QfljOO  PRINT*,"  point  no.  a* 

0034-1  CALL  FIELUITETAOO^hIo, 

00400  PRINT*," 

00410  PRINT*." 

004PO  t T 1 *C  A«S CETET A 
0043O  tr2*CA0s ttTETA, 

00440  tTO*CAdS  iETF.t  aD 

00450  0EPlVT»ItT2-eri. . , . _ 

00460  0EPIVF*ftF8-EFn  /(TfeTA02-lETA0iJ 

004 TU  CT liUEHIVT/ tT£fAO\-TETAOU)-t£T l-£TO) / ( TtT AO  1 -TE T AOO ) **2 

00400  LT2»0tWlYT-2.»rtTA0l*CTl 

004  90  CT3*£T0-Cn*feT  A0Q**2-Cf  2*  TtT  ADO 

■ 0£PlvF/tT£TA0l-TETAoOJ-(£Fl-EF0)/CTtTAUl-rETAO0J**2 
0£RI'/r-2,*TtTAGt*CF  l 
£?0-CF  ^»TtTA0O*«P-C^2«, 


StF2*CA68  (tPP!2 
StPU*CAo3  fF.PnlU 
/(TtTAC2-)£T AOi 


00500  CF 1 »OER l vF / 

0351D  CF  2*0  _ . 

0052*1  tf  3*£F0-CF  l»ftTA0O*«2-Cr2*T£T  aOU 
0O53J  ET*CT  1 *T£Tao*»2*CT2*TETaO*C  T 3 
0354  0 ef*cf  i*rerA0**2*CF2*TtTA0*rF3 
00550  PRInT  l50,£T,fcF 

00560  150  F0R4AT(5x,BritfcTETA)«,F9.4,5<,7H[£PHIJ«,F9t4) 

0O57O  CONTINUE 

00500  200  CUNTiNUt 

0054«  STUP 

00600  tNO 


SUBROUTINE  "FIELD" 

THIS  SUBROUTINE  COMPUTES  DIFFERENT  CONSTITUENTS  OF  THE  RADIATED  FIELDS  AS  A 
FUNCTION  OF  THE  OBSERVATION  POINT  AND  GEOMETRIC  CHARACTERISTICS  OF  THE 
STRUCTURE  AS  DISCUSSED  ABOVE.  THE  FINAL  OUTPUT  OF  THIS  SUBROUTINE  ARE 
"ETETA"  - EQ  AND  "EPHI"  - E^ 


WW5  LOMPCci ^^AL^TbiiSiA^ocPuJpfifcb.VSISciuJo^cS 

00o3O*, t*i.)G€.t«iTtTAJlt*PMlt*(iicJLT41t-><2#26eTA2tt3T 
0064O*,O0,Ul,tTtTA,tPNl*6*^t‘Tl£.'iA'l,ciftT*,6uPrtI.tV 
5*  ? 2 1 * A \ X r1 11 0 1 ? * P Z U F * H 1 » *■  u * * 2 . c J*  6 i \ , c .1 J - 1 2 , L j p N I 5 . C J i 
0^66-*CJ4>3,NF-ii,,-irN2,nfa6,n*Ti,N^T0,NrT3,!:JAcl,cJiL2.Ci 
0 //♦Cm»pp  i , CMA-ip  l,*Pu'tC,fM'oT#^ti-P,^itT4M,Kj\cr,rN 

i^.V;*C^:V  #C2F?L2.C2FPLl,C2PPL2,C2»'iLi,w2iNL2,L*20PLl  # 

i5*S^fer‘,UJ5eTV  12HL2,  III.  l2t  ,tJLi,cJL2,ci 

00  ? O'?*****  4 a f , r>  41  MZ  , 


CJZ1 ,C JZ2, 
CiiFu 

f ^ \J  N £,  r t 

C2mPl2, 
M4R44X, 


rH I S PAGE  IS  BEST  QUALITY  FMfUMAUi 

FROM  CXM*Y  MNUSHB  TO  ODO  — 
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0*71.4  MEAu  U,<Sl.Ol^Pri  (14)  ,31f>Ti(N)  , M8c T A , MZ , NPMI , 83 II , RSI2 

S)f§$  l31*Vul"*t (§i^ATt}*C^Ci»^*^^Al«.^?,a,SX,SNi.2«,F5.2(5K 

047604.5hmphib.ps  $,5*  3moa,F 6.2,55, 2nca,f6, 2, sx.iMAa.F*. 2 J 

00»7*  ttTAaTETAO*PJ/100.SPHlapnIOaPI/180, 

007  80  UU»MNf  TET  A)  *C03CPHll»V0»hlNifttx)  *3  IN (PHI) 

00  79  ) wiiCOSfTeTAjlAASBa.^PltilNiTtTAi «A 

00000  UGVO*UO*tfO/ (4,«0iyiuQNQaUU*'*O/  r J,  *PI)  SVUWUavOaRU/  C4,aPI) 
00414  WU01»-(l-VO«VO)/C4.*Pl)*»<U0ia»Cll-aU«^OJ/C4,«Pl) 

00820  LiaL/2,aZOSL2*L/|,«Z0fdtTAlaAr AN(U2/(PX*A)} 

00434  4tTA0«4TAN(ul/(PlaAJl 

00840  8lNT«SlN(TtTAJJC03T«CU5(ftTA) 


COMPUTATION  OF  WEDGE  DIFFRACTED  FIELDS: 

FIRST  STEP:  FIND  THE  WEDGE  DIFFRACTION  POINTS  BY  SUBROUTINE  "WDGPT" 
SECOND  STEP:  COMPUTE  THE  DIFFRACTED  FIELD  CORRESPONDING  TO  EACH  WEDGE 


DIFFRACTION  POINT. 


PmJnT*?1 
•irtAXaS* 

wOGPTi 


COnTM  18 JT I ON  OF  *£0Gfc  DIFFRACTED  RAYS’* 


ILL  W0GPTUL/2.-Z01  / A,  TtT  A,  PHi|,NMAX,OlFPTi , iROOTl ) 
flNT  112,  (oUPtl  U)  . J«l.  iWOOTl) 

4 FORMAT  ISX,  34rt0!FFRACTl0'9 


POINTS  JN  T8t  UPPER  *El)Utl,2X 


1973 
>984 
00990 


00450  PR|NT 
00844 
04470 
00840  CA1 
00894  PM] 

84900  1)0 

BRiaa.dFR.S) 

00920  Call  aOGPT  t («l/2.«Z(J)  / A , Tt  T A . Pr«I , N*AX  , 0IFPT2,  IR00T2) 

04930  PRINT  124,  (OIFPtS U) . IMu6T2) 

JJ440  120  FORrtAr (3X, 30MOIFFMACT IJn  POINTS  On  r*t  uJ»£M  *EDGcl,2X 
0344O**  4F9. 3) 

004oO  lrfOOT4|MOGTl*|R’JOr2 
fe*TtT«£*»F  l»0 
00  200  1 ■ 1 * ImuOT 

If  t i , gt . Iroo r i ) gq  to  150 

01OJ0  PHlCaOlFPrlll) •Pl/140ti««L/2.-lO 
01010  IF  (lA4a(P»lI»P«lO)  # A I #*I/2«)  « ANU  , ( TE  T A , G I , P I / 2 , ) ) GO  TO  173 
0 1 420  Go  TO  1*4 

41330  150  P«tOaUIf PT2(I-l*U0Tl)«P W140.*H».c/2.-ZJ 

Slip  .{Srifs^i3S?iii?iJt5),-*}iA5i2i'T£T*-LT-'‘,i-M 

iissi 

01104 

Sit  S 

01130 


Kj'SttJKitfWMIiSf 

8|TA«ATAN>(ri,  A*PnlO) 

IF  (slNC rETA) »CO3(PMl-P«I0j ,N£,0,)  GO  TO  174 
168  CONTINUE 

PKlNfa,"  OBSERVATION  POINT  13  On  JNAOOm  40UNOa«»" 


GO  TO  173 


A ■•RPM I • 3 1 N fbETAlaRZaCUSlScTA)  »8£T  AlabETA 

a,EJ,PI/2, ) ,0R t C8E1 A ,t i,0. )>  dt T A 1 a 4£T A *4 , 40440 1 
SlvlaA/C04(86TM  Jaa2lMH08EraA/IlNl4tTAlla«2 
a (P I / RriOS IG a #21  l «/3,J  aSIGM A 

GNAa  (4,.-0.S44j  *Ct  X*  l(J,,»l,la2,aPlaSlli*A)/SIt.NA 
"8ETA«T  f 1 .•14. , 1 .) / ( 2,*Pl*3lG"*) 1 #F JC A V (AS  I) 

15Ma)  aa2«(4tl  t . ) ajPuCav  (441 J / Uaal42l  • 

tl  * H • 4 l • ) Jt*0  18»  1)  a (RnG31»/*NwSc  n 
J*2.*Pl*<?No3i'»i«a(2,/J,)J*aSiGNA 
„./U  ,<Si'»'‘'*aSlN(  TET  AJ  aCuS(P*l«PM  J)  /l  AaSlN(0t  Tall  a«2) 
t 4«) SRhP«A45 (RhQ J 

41224  feP»UN«lt*IPl*8J*(PMl-PNXJ)  ,tl  , PI/2,)  iPSwNa*), 

4iaI>Ta  • P > u P r * • 


t150  r 

1160  CiAlFAa 

II74»-FOCNU(nS1) 
J 130»Pla2,aP.‘*u5l 
1 190»/  tS',4  r (2 


01200  *"Oai 
OUto-H  *MhP 


•*(2,/3, 
~ I *4"<0; 


012  3 4 irt£TAi«3l"« 
dltfUO  l6(AiS(3lN{c%tTAi 
4l25<)*)  aSiGMl  ,4,  SI  «idl 


i ) ,u  , A8S  (0  J j i Tc.  T A J J J 

. r f i ) j 


S8ETAtaA8SiCua(TtTA) 


THIS  PAGE  IS  BEST  QUALITT  PRAOTIO^m 
FROM  rOPy  FORMAS  HMD  TO  DDC 
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(31060  *3  [•►.HSU**  *4  3 (.ACC5  ICUSC  f t r A)  /*Ar.T  \1 *1 .0} ) *1  ,?lA 


)• t-i .)-*./ ta* 


31  H-3  lr  ( HJ4  fr'ol  J .LA, 0,0311 ) ut  Tc  l q 4 
UUh  Irv«03(C<ja(PSI/^.)j  ,Lt. . 0,  o 0 .100  l ) 

A 123  0 on»b,4£.i,ijj2/il,343,*LuSl<i.*Pl/J, 

U10J5*  U»Sl  /<?.  ) 

01  <240  ViH  s 1 , juU  10  !M 
3130/)  1 6 4 0*»b,  9282032/9,  SG9»1  , 

3131')  1M  CO  * I l 'iu£ 

*31320  u>*e.4A2««A*A*jI.i(McTA)«*?/<*Mu 

0133 «’  (."£uA3«4*CU5(JtTAJ*(3.*A*A*.jlN(:1':J*)*«2/31li.’,A**2 
01i40*-l,j 

0 1 5 5 0 l'^r  '>  4 4 ■ » 0 i£  G A 2 * A • A * a 1 •<  (,  5 1 T A)  *•<*•(  1 ,*3»*A**/SIGiiA**2 
31Jb<l»*15««A»A*C03tPtTAJ**3/Sl<i'lA**3)/5>l'j'iA 
0137.)  IF  IA43(U  'ti»A3J  ,LT,/,i>n0000t  J 4McGa3»0.0«)0vi0O1 
0134  ) IF  i Af*s  t j^tijAai  tu r , 3 .'H) 000, :u  U"F.ga««.3  .yotfvH/r 
01344  Tii3«AiS(u«*ttaA^0idi*3jii24*Ada(urtL1iA2,-s--><-  J 


A 2 / 0 ‘"l  6 ii  A 4 ) 


IF  THE  POINT  IS  IN  THE  CAUSTIC  REGION,  USE  EQUATION  (7.83)  BY  CALLING 
FUNCTION  "PHIDEL"  OTHERWISE- EMPLOY  GTD  FORMULA  TO  FIND  THE  DIFFRACTED  FIELD. 


i - 


I 


0 1 4 .i  a 

A l 4 1 J 
oi  4.;* 

014  j< 

0144  J 
01  ’4  SO 
0 1 4 D 

alar,) 

0 1 40  if 

0 1 4 4 o 

01  »>1* 
OlSio 

0152/ 
0153  J 
0154/ 

0}5?‘ 

Olio/' 

015  N. 
01  5 i. 
01590 
O 1 o J •’ 
0 1 8 1 <•* 

0 103-' 

0 1 O 4i 
01634 
0160  0 


IF c (T23,Lfc,l5.) . ANJ. (T?4,t£.  10.) ) GO  TO  1777 

£'|Oi.fc»uH*(i  ,«i  J • r (y»ir'/2,j  /2.)  *C£x0i  (*.*  *1 . ) *2,«i*I*  l 

♦ a*sin  (^tT  a)  • ecus  c^Hi-Phio)  -cus  c^’ii  n ♦M*cna  ( rSr  Aj  j 

♦ J*OH*tiAi.F4/Ul'«(A9S(d«;TAi)J*0.**sl) 

IF  ( h M J , (.  T , /,)  t 'itJGta  (0  , , l.J  * t * t !)  6 c 
iG  T J 1 7 I 3 

1 7 7 7 UKTi-MUt  **»Ht  4T«,  " CAUSTIC  0 T 4?C T I 0* » 

6 *F.i;iic«A*H(««*Ct*»»  ( ^ '1,  , l . J *2,  1 * l A*SI  \t  ^TcT*  ) * (C  J8 

♦ VJ«Iu»unIJ»r;j>»lP"<i)  {♦"■CuSlltl  4 ) J ]»'J>i*tl4uFA»r"iIO 

♦ CLiPi*  O^cGA^, Hi"  4*3/3, ,Pl«'-fifcU,14/12, )/(«*. *01) 

1/M  CU.'*Unol 

t/i  Tt  T 4»ot  .i6l)Gt  *Cno  {HhJ-PMIu)  / i 1 < (,\tiS  (St  T A)  ) 

C -Hi*:  a*  vt  Out*CuS  (.  r t 1 A)  *ai  '«  IH*  I -►>.-»  1 J J * C-l  ,00^000)  / 

♦ a IN l A ) S ( 01  T a 1 ) 

f*i-  1 1 130,1,*:  •T£TA,t*»>*l 

♦ *Y  NU«»  ia,SX,BH6T6T4«,F«i.3,F9.5.S*, 

6U  1 u l 7 4 

17a  '4cr.i,.Af»r.i,j  ».pKi.,r  177,1 

177  r 4<,  »t  a T 13*,  7hkat  >0 . , 1 c, 2X , 2«n  I S JlOCM.0  8T  Tit  UAJtCTj 
1 7 7 Cjv  i i -lut 

c*Te.T«F.  «!»  T*c./<T£TAsfe.-(KI«C(HFi4t/<04i 


2 < - C > *6  T l <bt 

041^1  l;<5.t’’7&T,£rtF'l 

1 o S*  0 4 <1  a f \ \ i a , 7 .It « r c.  T a » , F 9 . 5 # r 1 / , 5 , 5 X , 


AainT  2 jo 
210  ►‘04*1/, T ( 1 J x , / ) 


6 wFl  /i 0 3 1 a , F 9 , 5 , f 1 0 , 5 ) 


COMPUTATION  OF  THE  FIELD  GENERATED  BY  THE  "BODY"  CURRENT: 

IF  THE  OBSERVATION  POINT  IS  IN  THE  DEEP  LIT  REGION  OR  PARAXIAL  REGION, 
USE  SPECIAL  EXPRESSIONS  DERIVED  FOR  THESE  CASES,  OTHERWISE  PROCEED  TO 
COMPUTE  SURFACE  DIFFRACTED  RAYS  BY  FORMULAS  (5.1)  AND  (5.2). 
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.1.7,1,)  GO 

, tt  , l / 4 , ) 


TO 

, >jK  , (2*  I , Ut  , 


01570  IHM3.L 

0U60  1FU*M. 

•non  p-xivTa,* 

017  1*  ^ttfcT  * 

0 j 7 10  P><  |M  • # " 

4 1 7 < 3 S1G  «*l  «*•  (P"I-Pl  /*,  ) / COG  l J£T4)  Si<Hyj 

01  73J  A5ti«iPj*AaCCdv.ofc!An*4U,/Jt)*i2H 

Jl74*  tNi«'i.ieT4*i.'1#  t . } * *,S*Cc*.3u'»  • 1,1 


. 7 , «PI /«, 1)  GO  TO  450 
COif^lrtoriUN  0?  Subnet  9AY3  " 
J«tT Ati'PBl  *CU5  (T£  Ta)  ♦,i/*aI'<CTcT4J«oeTA«PI/rf.-TtTA 
Nl\!,,'  «<4f  NO,!" 

COG  l J£M)  JG<No:JIGaA/CuS  (4fcT4)  •«* 

, , , ...  1 1 1 *2  • «►***  H«aIli*A  t «S1M  c TtT4 

l790*l»42*««alN(reTA)«COa(P4;))(«u0(A3li)*ftP’iI«31<\i'3cTA)«CbXP(t 

01  T70*(Pl*4»ni!»l3)  ••  1 1 ,/  J - J «G\  (A  A t U / (4.*Pl*A) 

0178/  trtcl  4i«'*»»Hi4toaI,li:i4)*(fil*r<HOSlMji 
01  74*1»2,*P1*  i 8 1 U M 4 1 • Sl',(ltl  4)  * < 

4l:lf  - l M-  *40  tfcNl ,t*tTAl 


"111  / r 

• (i*i4(<M0siGj**ia./3,)  *cc*i»(  c o,,-i  .)• 
4)  *4i44*3iNC  reTA)  *tCo  ) 1 AP  4 l ) 


' 'i  l ” i ris  ><;nt  < • l 

010*0  *<»0  *0  4.*' AT  t 10X,  4ht  10.5,5X,7*Et4fcTAii,*9.5.M0,5) 

01430  2 5<*  04I*T«"  «•»  »0,2M 

01640  4liinA^«AAC314pl/^t-r»4n/C0Stt}cT4) 

01*50  A5ia«ioi44*{os(oci4))4*ii,/i,)4ij.4Hi/a,«HHn 

>0,5*CbX0(l«},*«lt)*a,4rl«(5IG 

, O * 1 \ 1 l ' • . . . J 1 ^ I _ f . _ f v # 


CN*«i4»€  T A«  J 1 J j ,0 


01*30  tN*«.4r‘irA«f>5lfi,)«0,5*ccxi>uo,.-ii)**.*ei«(srG4  4a«siNiTtT4 
0io7d4)4«2^A4$l*OftiAj*CoMP4l))i4GtU,\5l*)«MPwl4ai'*(dtTA)*CcXP(l 
0 1 3».  3 ♦ * , . »l,)**,*Mi,(alG''A0«aixiTc.f4)4«^44*dl'(ii<ifA)*COSiHnl)))* 
01©404liSi*H,lOSlvi)»*ll,/3.J*G1143l,»)/l4t*?l*4) 

019o0  taETAaA^^nldCoaUtTA)  ,CtXt»(  10.  .-I  - ) «a.*r,l  * C31GflAa*SlNCTtr4 

0»<>iOA)***!AA.5i'i(itr4),v:oo(ij4i)))*i*,i«'4hOsi4)**i,a./3.)»r0i\sia) 

oi^ao^/u'.A^iAAj 

0193O  P41M  OTI.feNa,*;  3K  T 40 

01  «40  a TO  FuS'iAt  (l0&,<intN4«(F9,S#»-  1 0 . 5 # 5X  , TnEhfc  T 4a  a , F 9 , 5 , F 1 0 , S ) 
— 04i\r  ait 

350  Ck'NTlMvJt 

tSTfcTA««(t(S»'.TAl4t!j6l4a)  >6  5^nl»-ltNl4d42) 


0T 
11940 


950 

?rJ 

a; 

01440 

3*000 

a*  oi  0 
0*J*0 
0*030 
0*043 
0*0  As) 

0*05.3 


P*I\T  J«.) 


350*  J4  -« 


o0.tSlar4^t5H*ii 

AlUO4#7*'i:iftTA4#f4,3,rl.l,5,5*,o^tv'>»,fiIa,r4,3i's10.5) 


‘4  1nT  *io 
r*H  J \ I « , * 

ertT4a^ilr.T*4txriir4cp 


!'.»•  [SlTc  CTL  l*ot* 
1 1 *ca*,<il*c  «r  X ♦») , 


0*1, \T  S8/#tTt.U,tP4l 

T6TA«,*9.5,M  3.5,5<,5"F.PrUa,F4,5.Fl0,5) 


3«0  ► U<4MA  T l 1 ox  # fertc 
* JO  C O',  T 1 \^C 
40  TO  51-1 


THIS  SECTION  COMPUTES  THE  GEOMETRICAL  OPTICAL  FIELD 


02070 

0*040 

0*34' 

0*100 

liw 


45/  CJ*TlN0t 

tu T 1 1 4 a ( 0 , # • 1 , ) • 4Pn I *CJ3 iPnl ) , . . 

tl<H-la(.»t,l%)«C-,4.-<l«SINi/iU*'tJ5(T6T4)-^2«ii'<Crr.  TAJ) 
fcrfT4atoTc,TA4c«idf»t',,,,l*c0','il*t<>M 
POl<T,,"  OiPfcCT  4Af  CONT*lBu TIOm" 

vcic«<  PkInT  4(*0.eoTtl  4,f)PHl 

0*13/  4 50F04s4r  C10X#  7HtUTtT4A,F4,3,f’  1 / . 5 . 5 X . b«t  0 
0*140  PKlnT*."  iMflsllS  CtLi^OcP 

0*150  H 4 I * T 350»cTtT*#*,'^4i 

■»  T « L 1 W • 1 V t 1>  i .k 


OiPfeCT  4Af  CONT-ilriuTlO'j' 

b«tOHhl«#F4,5#Fl0,5) 


0*150  510  CuNTlNUt 


COMPUTATION  OF  TRUNCATION  EFFECT  USING  FORMULAS  (8.30)  AND  (8.31) 


0*170  PHI  xta, * T no nC 4 T 1 On  erFtCt" 

0*13/  '<P  T»  j.'.a3#< . f i(*«Pll  i (i.  1 , L*  ) ••*)  )M  . 

0*143  )P^iPa2,*Pl/ ('>  Ngs«t) 

j **oo  :wM.a«u*'<lu’ii  ii)«o,« 

O ? 2 l o t.#M  iai't*Pi(.'lf»l,l**,«Pi«Gl*wx.dl) 

0**2  0 rk’*'i  1 a *,  ‘ 

02  2 i v CjM^TtA^aPiai'  l«Gt:xPll>  ,#»5.)«2 


J*24. *(1  ,/3.‘i 
*2*50  JU  a O o IT4.Nh'ial#xPT><j.'; 


I/l*«)  /115*J.*(2,*P**4)4* 


JBISPASE  IS  HKM  QUALITY  nOnKUMWA 
IBOM  OOPY  KVnA.oiiaU  TO  UOO  ^ 
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olvKMi 

► 


^ 2 4 ft  c 
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32<*«. 

W2240 
323.)  * 
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Ar*Hy*))  DJ  VI0«2 
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'•'■r.  | »V«wr  *r-  I '^CTWf  J J " V * (I  i U • I V*  C *4  Jr  n 4 ^ 

CrtJ»rIa$,i4?(<i*A»4aaU»niPi/al'tT<«*j5*CtxP((u..»l.)*Pl/4t) 
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c,_.  ( c**»>.'*rK) *«2 *l2»l21  i 

' «c< 


] * ( Adi  KeJ«2$t»"  l ) ,ut  , ,35  J l ■£"$£ I 

ifl-Ll.td.O*"!)  ul«t  1 *3,2 


2l. 

223a 

«2592  Zti  * * • 5 « * t (H!)*sukt  i a p j»  i 


•'•5j*l(H!)*au*TlApJ»l»A«*aatPHlP)«;>JlNT*LlaCd$T* 
l (AAHnJp)  «*2*L\*L1)  J JaHWMl  ..LlAZSPru) 


324?<)Ot»4  I , 

<2 24  i »:  C"A(«r 

32422*At'5tCC$f*ci/a.i<T{l**P<4lri)**2*k.l**2)) 

3243  3 I ► l *c  s (L  1 ♦ 25“"  l)  ,cfc  . , <?*»)  Ca  aa2  1 ■(.'  *>/  1 
3 2 4 4 rt  aIl*'"A>«i«NT  l 14*P;aIP)*«2*Z3P'A(*«2) 

2 2*0  0 A ( «5  T 

224fc.'  l (A*Hiilr»)*»2*U2«l.cJ 

224  / ,»  Cjap.,»«A*»>'<i^/siitfi4i3sr'.ai«Za»3^I/Si(»flAsii;s,iS«*»coaas**2 

32  4*1-1  (lwaol»A*^*iIP/SlGlA}*il'tJl«»Ul/^lp'''Al^CSOl^«CJJ»dl«*2 

324  5.)  uuSi^«,\*P«Ir*/aX(Ji,4  4i»S(  xo2*l.2/al»i«A^ 

325.10  Wfu)TS«*/UOSaan25r.».oTl»*/CJ&iM  ■ •2»4rtGTcl*4/CUa32**2 

325  \ o 

325*.'  ’v3Iiai*l/!4HDTi««£)*«ii'/3a)*M'j'’M 

3 2^50  t'ale!*01l/<nUT2**p)«*fi,/j,i*aLn.', A^ 

21 25-4,*  »^;J6«Ci;P2*Aat/3,)*Ct<P(-,3 


)/ 


32$V  ♦CONal  / a j*  I (2.*PiOl'jr<Aa)  *ir  (alb,  Aa,tt,<i,  1 )*r  3S».lt  l-O.n’UlklOO 

j»  ^ ^ M t ^ M t I « i i » - r . » j i _ j .i  * - i 4 * • . 1 7 C . >•  . ~ - 3 . x . 
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.diki'jll-iJ,,  l,)*,il»kSU)* 
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1 •’ 
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02*7  c 
2264-' 

3 2 a *7  •) 
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02  f 1 ■■) 

02  1 20 
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2J 
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• C l a T A aj 
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r r»«i«£j\c2 
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► ...  J . r / « t 


32  > *31*  (Cj^^ 

2 2 » 3*)*Cr.  Af*  ( ( . . , . 
22^o24g2CI  Ar'/aWrtf  i<? 

22  7 7 04*CJALl*Ct*P( ( 

22^ai'**CMari«CJ2Hia 

22  7 P-’.  + C JPfi 
22t»224Ct  »►*  ( ( 
32«M»'*yZtlA-/5w4f  (2*1 

“ * “ * r v.  J m T • " 


?2<J2  .>  lr  (Mf»nl 

fltfrtj.’  r Z«C  j-»u2*Cc<2  ( ( 

•?  2 o •*  £ * k «C  i r i • c j z a • 

22  a >j« 

’2^  ' -UV  M,  1-1,1.  »■  V.V  - •*-(. 

^ 2 ? ■«  i/S  ,h  i 

*2'J4.,4^«a  ••►  ( «C  J 23*i:  Aar  I)  / i c , • i 

.*"(  •*•*  u-a(wrlr’)*'3i  \ T)*i.Ajr  l*„jZ3A.4tT*A^  j, ,1(^5 

J 2 "7  l J 7d'a  1. 1.  \ I i i j t 


3 a ) rnJ  To 


i)  r(«CJHi'kC»ru  '«*»U1*£«*'4»suai  J» 

) * £ ♦ k • C a S r v • ».  J l a • >itur/j.  * f l *; « • " I ) *C1i » P ( 

•ij«k.'*arA*CJ4a*o2r-l«r'/3-«Ti2,) 
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YHIS  PAUf  is  bSST  4UAl*ITt 

ooi'i  fvwush*  ro  i»0  — ^ 
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32923  *>»"UaA*»''li*fJi,«Ir*/OlrflU'aAaf**MaaI\(F*l.«HMlM) 

02930  F>»"i«}aFT»r'l2MGr-nlH/jlV 
32940 

a29v>  d.'v  c*.*\tIni'& 

22  9«v  f 'tu-t  I «Ct*F  (lu.  ,-l  , ) *2.  *»*i  *A*UvJ  • to,  • «t  • ) • l at1,  »•» I J • l COST* 
329 

329<*0  f «<uNtF«CtXi»  ( l4t.*liM?.aMa*ayO)al«),,-l.)a(oO,*PXjaft»«I2 
«f2°92  P4i«T  t'i'4  , Tntj.se.  1 , i ouNfcf 

3 3 30c  7 99  l-j»<>'AlCb<,7«TMU^T».F:J,4,f9t4,5x#7HT>?n.>1efta#Fd,4,r<»,4) 
SUBTRACT  THE  FIELDS  DUE  TO  THE  REMOVED  PARTS  OF  THE  CYLINDER  FROM  THE 
INFINITE  CYLINDER  SOLUTION  TO  FIND  THE  FINITE  CYLINDER  SOLUTION: 


ai? t j met a«c rt t ii-r^u.«tF 
al.Vc  ^isi»  « 

*3330  P*l  ,f  sAtf  ,tTe.TA,SP*i 

?j2U0  t»w  to  900 


rOMc  f-'IfeLO" 


THIS  SECTION  COMPUTES  THE  FIELD  IN  PARAXIAL  REGION  BY  USING  EXPRESSIONS 
(9.20a)  AND  (9.20b) 


*304c 
O3.'e0 
3 3 « 7 0 
a 3 o o 

330  9*’ 
J>  3 1 f *' 

2 3n  * 

03  1 2o 
*31 3o 

2314,* 

031V 
o 3 1 *>  a 
?3  17 
o 3 1 6 t 

/J5  Q.* 

A 3 2 0 *. 
< 3 1>  1 . ’ 
*322c 
*32  3*. 

2 324  0 
032*  <: 
?32o0 
0327a 
P3220 
*329*. 

0 3 3 »'  ** 

0331  o 

1 J3‘»«i 
C Ji  J ’ 
v 3 5 b « 

*3*  7 * 
i.35**. 
} 3 9 <■ 

3 3-<:* 
0 3 M • 
Oi«C  ■ 
03*3. 

.'344.: 
**3  4*1.: 

0 * *nij 


610  ?91>rat"  H4*a*UL 

a*  2,*2IaIF  (466(Tc7A«^I)  ,ur,rtl  4)  m-2,«P1 

^lN<scl«aI't(4*u.l)nCoa4Ll«tCS(4*ulJ>SiNaL2»3l''(A*L2) 

CUS»  l 2 «CoS (a  *L  2 ) 

121  a 1 p aO, 

*IaA*4wui.l,L2,»)*‘‘.’.MSC,*ZaAa3i*)/(NI-l) 

2 J lit  1 4 2 « l , > 1 

I y i w 4 • l . 

if  UlA7.t4fl),u*,(I*4tfc'J.  <I)  J 0IVI0*«2, 

*Z«  ( 1*2-1 ) 0*23  If  (*Z,c  J.oatM ) AZa*Z-.0301 

*»4«a4»<T  i(v#A»*4a<vzi 

C alw  "JEascL  ( » , ► *,  *■  \ ) 
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l • a a r ) 

*-.Ju2*Ci:4^U»'.#-l.i*(44«K)*!.2)/('i»A) 
tjLl«uE*f,UO,|-li)*(A:A|*)*Cl)/i’vZf'V) 
iH«In*Cl^4/4iVI0i<J*(*Jl.c»tJV.l)*7-2AZ 
1*1  ■ 12  1 ♦ (U*Z/U  l v Io*J  a (EjL2*e  JLl ) a*»2i\Z 
lit  LJ-sTlNCt 

COF.'il  l»C2f'*v.?»02F2n*C2fF(.a«C24-'Ll»C2ij'4L2»Cai»f»l.t«C2jf'L2«.'l 

2 c T 4 ix«-AuOi»  1 2.0  031  )/A*lM(Ll»C2>vJZcTAai»l/U0*«*A) 

OZt! A«i2.*u2ET4 


oat!  A2a2,*u2Er  A 

,.N2Er  A*2tfA,4»/v'ZfcTAiIM  ?*4cr  A,*»!  ,bCV-n  -*ZET  a«*>iO 
20  022  12*1,  \ 2 c 1 A 
Gl  viC2«l  .AIM  U.tG.  1)  'UUDIIZ, 

At  TA«  ( 14-n  *02t  t ASlc  (ZEt  A.Gr.itTAMX)  GO  TO  3J3 
c2H»CKM*2Er4*Ll)»tZu2*t.vM-2fcl**L2) 

IF  C A "IN  l (fcZLl,fe2L2)  .LC.O.JOl)  '’IC.TA«lt*L)Zt7*2 
Al»»*C)baLlA,'tT4*ai'**VLlJ>0l«2tTA.CJv»aLl«A**l'»\Ll 
*«»*•(-  s|'L2*2::r4.c>i  'M.2oo,:*4ti  A.uo4M.a*f'«ai 
f l<  »*  r J MI4*o  .1*1  it) 

C A u L ? t ■*  S £ L t » * ► •-  , f * ) 

A w 1 r.  f 4 ■ t ■:  1 .1  »f  t ( S ..  •>  r ( » ) » r •> ) ♦ 1 >,,l,)«2»a2M4/(Ola*aY*fw) 
ij  2 1 c.  T »»:•,.•»  1 irj.f  \ / K ^ ♦ 1 1 .)  •2./i,i>I«Aaf'0 

Vcc‘Ll*Ccfv*l.l*('t&r4/LlTijZ)*(AJ*sJ%,l4l*^l)*fJ2tlA*tc..l/' 
Uef  'i.«?«CrrvLO*i  ( C T A/..  1 » ) .‘t  1 « l >'c*  i-,  . I . ) J/t  I 4**  ’ 


*'l.  1 *Cr.f  1 ♦ (. 
'^iCer'LOM 


{ C • 4 / i,  | 


liv7iiHU«i;i;!w.W*  'AJ  • 1 a \ • ( %,i.J*-*li*AJ2ElA*e.2wl/i 
1.V*  -*.F»cOf  ^'-c-  v t - 1 4 / 1 » 1 . 2 ) • 1 *2»  * ■ . , i.  ) *o2  ) «F  f 2*  M*c2i.2/  ’ 
COG  1 1 » C 0 j 'l  1 ♦ * w 2c  7 ■»  t . 3 v i .».•«(■»  1 - ( ’ . , 1 , J • 4 t ) **.•  J 2L  T 4 •►.  2 . 1 / * 


Zc'.\/  i m.'I  1 , , t , j «*,>)  U«Cih2/  r 

C«*j'"Llai.*o!Jul*l.  2t!A/v.i''iJ*';ai^lA(.'.»l.J*Alj*GJ2r.  1 * •£  2L  1 / T 


* 1 c > » l 


1 . J • * «5  1 


.2^2/  » 


Cc?«>T'CeaCa**fua^l'.  / cTa/, 


ioz)  1 . t m2) 


THIS  PASS  IS  BEST  QUALITY  FRA  OIL 
FROM  OOPY  FVflUUSHMT  Tt'  DDC 


203 


034/u  222  cOnTInI'I:  ' " 

333  C 0 n T I n u t 

iilil2»?2.*i  1 iAAi4Aucc2F"L2*C2»'"Li-c2rPL2-i:2KPLi)*(.ij 

93S20A"tAL  [C27,PL2-C2ri'l'LlAL2F  PL2-*2FPLl ) 

03510  l2LlL2«-2,*J2lAAl*AU(C2liML2AC2Gr<Ll6C2UPL2AC2GPLl)*(a 

®35i0*»<LALlCiG'-l.a-CrfG'iLl»r.?uPv.£*CiG*,Ln 

® 35 PArfAXXald,,  l,)*nPnI*ilLlL2/(8,*PI*Pl4A) 

0354*  PaHAXTB— 42«J2L1L2/  (0.aPI*PI) 


• I 1 • ) * 


IF  THE  OBSERVATION  POINT  IS  VERY  CLOSE  TO  THE  AXIS  BUT  NOT  EXACTLY  ON  THE 
AXIS,  THE  SPHERICAL  COMPONENTS  OF  THE  PARAXIAL  FIELD  CAN  BE  FOUND  IN  TERMS 
OF  ITS  RECTANGULAR  COMPONENTS  AND  THE  TOTAL  FIELD  IN  THESE  DIRECTIONS  CAN 
BE  OBTAINED  BY  ADDING  WEDGE  DIFFRACTED  FIELD  TO  THESE  VALUES. 


*355®  C Ttl A»PaHaXX*C031PH1J  *COS ( TET A) aP AP AX Y *8 1 * (PHI i •CUilTfeTA} 
93569AAfct»  r fc  T 

2 3570  t*»hI».PAWAXX*3lN(PNl)4PANA*r«C0S(PHl)  aEaFI 


93540  PP I NT  8e6.PAPAAA.PAWAX Y 
23590  066  FoPmAT (5X,3MtX«,F8,4, 
9360*  PM1NT«." 

6it>  U’  PNl^TiV'i.f  TtTA,EPHj 
93600  CUMiNUfe 

23630  NfeTOK.N 
93642  tNL 


F9,«,5X#3Hfcr«lF8,a.F4j.4) 

TOTAL  r itLU" 


SUBROUTINE  "BESSEL" 


THIS  SUBROUTINE  COMPUTES  THE  EXPRESSIONS 

FM  - — 5 = , FN  ■ J (X)J'(X)  + Y (X)Y{(X) 

J’^(X)  + Y|  (X)  11  11 

AS  FUNCTIONS  OF  X,  TO  BE  USED  IN  EQUATIONS  (9.17)  AND  (9.19). 


0365?  Subroutine  dE$SbL(x,Fi,FN) 

93660  Pi*3. 1 415920536 
O3670  IF  CX.ti.w.)  X40.0091 
03000  IF  ( * ,GT  3 , ) GO  TC  20 

93690  XTax/3.03xT2*XT*XTtXTuaxT2*XT2SXToBXT4*XT2SXT8*xT4**(2) 
03790  XT10«XTA*XT2  ' # 

037  30  OB J 1 ■«, 5-3 . •,56249965a A T265.*,2 1093573*XT4« 7,» ,93954289a 
03740AXT069,  *,9044331 9aXT8-U ,*.90031 761 *XT 10 
03750  CLOGXbAlUGIX/2,) 

03  760  dU*  (2,  / (Pl*t  ,0n«£LCGX*dJ  l*(l  ,/X)  «(-,6300140*,221209i* 
O3770*X,2*2tled270  4*Ai4«i,3t64rt27*»r64,3l23'7  5l*xr8«,.’4  0J47o*ATlO 
037dO*4 .30278  73  *x  T U*  <T2) 


'1 7 4*  w3t  l »,/ri«09/u^.  »)*«5,;'i?n00U.10?00*  (2,  / (PI  *x**  1 ) ) «dj  !♦  (2. 7 (PI  • 
1360^*1. 0)J*cluuX«u 4 OlA(l,/ix*x))*(,e36ol4d*,2212J9i*xr2*31*2.l6 
9 3al0Ad2709*XTa»5,*it3io,lo2/*XToA7,*,3l2395l*xf's»9,*,0  4i)04*a* 
P3B2»;*X  ri/Al  ) 02278  73*  1 r l/-*xTi) 

9 3650  r**) HJl**2Aonvi*#2,*7',*uc4l*dJlAUdiri*fl<l 
03«4w  4ETj*>\ 

03852  24  IF  (x,&r,2,).)  oj  ro  32 

9366  0 *tl«P,*XAXc2*Xil**2'»xfe3axc2«xcl3Xb4*xfc34(XEl) 

038  72  i?l*l.Al5#/i2,*X62)*i4i  75,/(24,*xt4) 


mk  CAlit  IS  BEST  . 
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wtlUKN 


i*  Cy-Nr  iNUfc 

H E T (j  H n 

6.90 


FUNCTION  SUBPROGRAMS 
"FOCKV" , "FOCKU" , "DFOCKU",  "DFOCKV" 

THESE  FUNCTION  SUBPROGRAMS  COMPUTE: 


v(x),  u(x),  u’(x),  v’(x) 


03970 

03980 


COHPtfcX 
CunPLtx 
OIMtN 


function 
F v - t P 
AbSTFl 


FOCKV (X) 


03990  OIHtNSION  AbSTFtlO) 

S33rS.?a:0s05i:;<:aS)Sl?tf:S9«2::78?i3f29O'tABS1*'-9T*1 

CA  li  !>  J u.  r m i.s 


0U0f!(fi  rvay 
04030  Flai,l4lb92bb38 
04(040  lFU,Lfc.0,0)  i»u  Tu  20 
0405H  00  I*  vat  ,K( 

rPaAdJ>TP(N)*CtXP(«f0.,l.)<»Pl/3,) 

*v«F  y>CExPl-(0.,  1 ,W*fF)/rp 
>0  CONTINUE 

F0CKV*Fv«3Sr<TlPI)  aSUftr  U)  *CEXP(-(0,,  i . ) aPi/4.) 
« C 

60  CONTINUE 

04l404«'4,004l4laxa49 
P fc  T u 9 n 


04060 
090/0 
04000 
*4090 
0 4 1.12 

::  u: 


0415* 
04lo0 
041  li. 
0UJ80 


04190 
iUt 


0k 


tNU 

ccMPcta  Function  »-gc*u(x) 

COMPLt*  fu,tn 

DIMENSION  Ad S T (1.1J 

Art»r/2#33eii, 4-38795. 5. 50056.6, Tab 7 1.7, 94417 
09112-MO. 1 1,43632,  12.  928*8?  ' 

04220  F^ao, 

PI*3, 1415426S3N 
iF(*.Ut,0.6j  GO  To  00 
Ou  10  \ • 1 * 10 

T\aAdST  (n)*CE*P(-(0,.  l.J  4FI/3.) 

CW  *F  u*CE  “ _ " - * ‘ 


32265, 


0«  230 
A 4 2 4 r 
04050 
04260 
042  7 0 
04280 
04292 
043P0 
04iU 


1 N4A8ST  CNJ«Gfc*PC-(0,.  l.J  *PI 
FwaFu*cExF(-(0,, i ,)axaTN) 

JO  CONTINUE 

fOCiuaFU*2,a50NT(Pl)axa«(3, 


HE  7 uw.j 

04|U  2F  FOtKu«l. -3N4TCFI)aCExP(f0t,l,J. PI/4  )axa*l. 5/2. 4(0 
^•^^♦•<**3/12,»5,*8wPriPl)aCtXPU0,»»l,)*Pi)4,J*X*a4,S/b4, 
043i<+-^-37oi*X**o 

34342  HjiTjKM 

3«35*’  fcNO  

043o0  CU^cE*  Fjnc  riu'.  uFjCkv(A) 

C u m p c c a OF  V ^ TP 


/2,)*CfcXPU0,,  l,)API/4.) 

• CEXPU0.#  l,)»f  1/4,1  «XA*  1,5/2. 4(0,,  5, 


ulnfeNSIUN  AdSTP(l*} 

AdaTP/l,  U 8 N, 3.24*20,4.9201 3.6.1b 331. 7, 372 18,«, 4884< 
3 J 4 5 » 10 •52*60.  i l,n7i«,n,  12,3 <.4/4/ 


043  / 

04ja.)  . _ 

04*0/  uaTa 
0440c*, 4.33 

044  lc  JMa,: 

04420  *183.  la  15  /26536 
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04440  i.0  10  MB  1 , 1 o 
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94  459  T <»■  a(J 5 TM  (NJ  «Ct  ***(•(«.,  1.)  *^1/3.) 

0446  2 1>F  VaQPvMl  ,-10,,  l.  J *2.*XaTPj  aCCAPI«1<J,,  1,}  AX*TP)  / TP 
04470  10  tONTINUt 

*44  02  pcocitao*  V«C£AP  (•(*•»  I •)  •F>x/«.)*SU4T0»  l }*X  ••(•!, /£,)  / 2, 

04440  HcTuKn  ’ 1 

34500  20  jFgC*:(»3,*:>w4l  lPI)*CtXPn0,.-3,)aPl/4,)ax*a0,5/8,*l3.»7 
0451-C4. ) «x**2/20.*o3.«5  JHI  (PI ) »cexpt  (0,  ,-l  , ) *p!/4, )*x*a3,§/102&, 
0452O9-0.  J2405*  A*a5 
«453»‘  4fefgRN 
04540  t *u 

0455?  CU*PLt<  UNCTION  OPOC5U(X5 
94504'  COUPLE*  l>  U » Tn 
34570  10*  ABSUiO) 

345B3  OAT  A A0ST/2. 33811 , 4 , Cfl7  9 b , 5 . 50356 ,6 ,7ob  7 l , 7 , 4‘*u 1 7 , 9 , 92265 , 
34543413.  W4017,  n.aodSa,  n.93o0i,U.oa«iM?  * 

H4602  OHJ40 , 

04010  P I a3  , 1415420336 
940^0  lr  (*.1.6. 9, 6}  tiU  TO  20 

04030  00  10  

0404  0 TNSAB3T  (nJ  ttElT^l-  (3,  , l,)*PI/3,) 

34650  OfuaOPU* (!.-(«•»  1 ,)*(2,/3.)*x*Tn)  *CtXP (- (0. . 1 , ) * * *TN) 

04060  19  CONTINUE 

0«6  7 0 UFQCAuaOfu*tExP((0,,  1 ,5 •*!/».)  *3.*SuPT  (P 1 ) *X  **  ( 1 , /g  , ) 

34oB0  WETjHN 

04O4w  0 9 OF OCKUa. 75*60*  r ten  *CfcXP(  (3,.-,  75)  *Pl ) *X • *0 , $♦ (0 . » 5, ) • 

04  7 904X**2/4,445,*Sm4 r(Plj  *CExP(C3.,-lJ  *Pl/4,)*x **3.5/128. 

947104»9.222i*X**5 
94  720  4tfuPN 
94730  fe.NO 

SUBROUTINE  "WDGPT" 

THIS  SUBROUTINE  LOCATES  THE  DIFFRACTION  POINTS  ON  THE  WEDGES  BY  SOLVING  EQUATION 


(7. 3) FOR  "ROOTD"  - $w  AS  A FUNCTION  OF  "HA”  * *2/a,  "TETA"  - 0Q,  "PHI"  - $0. 

34 7 40  SubnOuT IN£  wOOPT(HA,TETA,PrH,NMAX.«OOrO,  IWOOTl 
04759  DIMENSION  PH1OO(«03J  ,Fg  nCT  (400)  .ftOOT  (20)  .400TU  (20) 

2 4/00  rCT  (P«l  J,  ‘iA.TtT  A,P**l)  «1 ,0J*P.-«1  J/SOBT  (nA**2*P"I0**2) 

0 4?  70*»SI'«  ITET  A J *SlO  (PNi-PHlU  J 

247*0  Pl*3. 1 al5 

04740  CtLF  I J«2.*PI/N«AX 

94009  NMAX1*N4AX*1 

04010  N**AX«3 

0 4 0c 9 OPn I S3 , 0200 1 

04OJ0  r 1«FCT  (-P1,HA,  TfcT A , Phi ) *F jnC T 1 1J  «F  1SX  1 a-P  I*  I*uC  T ■«) 

34040  510  00  200  IPrij0al,NPAXl 

94iJb0  P«Io»-PUUP"10-n  *L’Et,n  JIPunCT  UPHlO)«6Ct  (PHIC , h A , TE  T A , 

3 40o2*Pfit)  *PHlJl3(IPfH0j  4Pril3*l83t/Pl 
04370  JMAUNeUlPhlun  125(1 92*  105 

0 4009  102  IPOOfilWOuNl  S*UOT  U«UOTJ  *PhIOSHOOTu  (HOOT)  «Ph  10*180,/ 
0 4 3 40  *r 1 

24  409  PuIoaPHlO^OPHl  *x  1 »Pr*lO*f  1 aFO  T (X  1 , n A,  Tfe  T a,  PhJ  } 

9441.’  IF  in  J 200,102.220 

>44C*  105  IF(p1*FonCTIIPhIO)  .CT,0.)  JO  TJ  UJ 

«•  4 4 i <■  <la-'Hig>h  lah.jsCr  1 iP"l 0 J 
94440  wo  r \j  200 
04450  1 10  X2«pnc 

2<*4oo  JO  115  jai-NH.nA* 

74470  ia(<l«x2}/2.»FxaFcr  ix.ha.Tc!  a , Pni J 


9440J  IF l F X aF  1 ) 113,112.114 

P4  440  M2  I»JUTal400f*  lS4ggT  ( 14  JUT  J a<3P  JUTO  ( IWJCI ) ax«  100. 
.* 3 <1  ■) 0 P-«I-ax*/r^I>xlaPnlJSFiaHOrixl,"A,T£TA,Pr*lJ 


/■>! 


■'o'MO  ir(Fl)  2 >* 9 « 1 t c • 2 ' 0 
05022  Hi  * 2 a x l a X l a x » f [ a f a 
03ci«j  uu  To  113 
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S33?S  11?  SitWir- 

350*3  P*  I O«U*0«*l)/lFCT(!l3,nA,lErA,**Hi)-K))*( 

35(370  UPM  l mOttl-  l U / 3 * * ( *4  M A X - l J 

0508*3  IKQOT  ■ l«iJUT  ♦ 1 S 4 U U T ( I * J J f ) «Hm  1 Jirtuc  T 0 ( 1 «U0  T ) «Pm  1 0 « 1 t»  3 , /P  I 

33{i90  *i«H**iu»0(»HinFi»rcrc*i,riA,  rer*,Prij 

351H0  300  CONTINUE 

«5ita  *6Tu»n 

/5130  fc.iO 


FUNCTION  SUBPROGRAMS:  "FO",  "GO",  "Gl" 


THESE  SUBPROGRAMS  CALCULATE: 


fQ(x),  gQ(x),  gL(x) 

130  CCHPtEx  Function-  6U(X) 

4 In  CUhMLfi  Z li  / l 1 H'  ) . Z 'i  i i’ t 


C5130  COMPLEX  Function  Sy(X) 

a>5  I aib  COMFLti  Zli/U  H’)  * 2«a  C 1 pj  3 J ,108(100] 

05153  OATA  ZGA/tl.4990.*  |fl0$5j  . (1.4497.0.0(359;  . 11,9 
05l6U*o3),  ( 1 .9940,0.00671  , 1 1 . 9946 . A , W 0 7 3 J . 11.9995,0 
051  7 0*4.  ,0064)  , 11.9994.  , l 1 .9942.  0,0090)  , (1 .99 

051  63*  1 1.999, 0.01  IS  J . (. 999,. 3 1*1,  ( 1 .994,  P.rfU)  , h. 


,4947,0.03 
.0,0070) , (1.999 
9991,2.0106) , 


353  73*0  44,  , /4*o;  , (-.3*43,  .*38)  # 1-,  3*94,  ,33 Jo)  , (- • .'335 , , ai<?b)  , (« . 2 
35333430, ,336 ) / 

35390H70O 

05403  1<10  FC4MA  r (6(N9.u»F7 ,4J ) 

35410  IF '.<,07.4,5)  OU  TO  33 
3S4*w  IF  (a,lT.-4,5)  GO  Tu  43 
35433  00  13  1*1,4® 

05493  13  Ztik  aJ»2f»ACI) 

05453  00  03  I » 4 7 . 9 1 

35463  AU3(nsCQNJGUu0(l-4b)) 

05473  03  CONTINUE 
354&/  Zu3l90)«Zv,3l91) 

05443  U*l/,*XilAt*46*I< 

35533  IN  C*,lT,J,OJ  GO  TO  35 

355U  :»3«Z0  ' ( K1 1*1 10,  «<-!*)  * UuO  U A 1 ♦ l ) -Z63  1 1 X 1 J ) 

«J  5 **  * 2 ® t T 

35533  35  IVjtC«l*  tUm.tJ.l)  I n 0 I C * 3 

/ 5 5 4 v UkmZt/(I*n*Ci3,**-I*)«(Zjca*U-Zu3ClAt-lNillC))l«Ce*FU3 
0553*’*.»  l . ) • <**J/i,  ) 

355® « PcT  -,<*  • 

055  7/  3 / v*/«CEaH  ( 1 ,,••  A 7 4***  i-,0®6  'e5-l  1,  , ,51  ) ) / ,545735*CEx.‘»  I J.*4 
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